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MATHEMATIK UND MECHANIK 


INGENIEURWISSENSCHAFILICHE FORSCHUNGSARBEITEN 
Band 32 . August/September Heft 8/9. 


Jahrestagung der Gesellschaft für angewandte Mathematik 
und Mechanik‘ 


Die diesjährige Tagung der GaMM fand vom 4.—7. Juni in Braunschweig statt. Sie war 


von etwa 250 Mitgliedern aus dem In- und Ausland besucht. Nach einem zwanglosen Begrüßungs- 


abend am Mittwoch begann nach der Eröffnungssitzung am Donnerstag früh die Arbeitstagung. 


Auf ihr wurden mehr als 70 Vorträge gehalten. Den Inhalt des größten Teiles dieser Vorträge bringt 
das vorliegende Heft in Autorreferaten. Die Zahl der Vorträge war so groß, daß im allgemeinen drei 
Parallelsitzungen stattfinden mußten. Freitag früh fand die Hauptversammlung der GaMM statt. 
Nachdem der im verflossenen Jahr verstorbenen Mitglieder gedacht war, gab der Vorsitzende, Pro- 


fessor Grammel, den Jahresbericht. Die Zahl der Mitglieder ıst Dank der hervorragenden Werbe- 


tätigkeit von Professor Görtler im vergangenen Jahre von 340 auf 50 angewachsen, 120 davon £ 


wohnen im Ausland. Darauf wurde kurz über das auf der Ile de Porquerolles von der JUTAM 
veranstaltete Symposium über michtlineare Schwingungen und über die Massaufeier berichtet. So- 
dann wurde auf Antrag der Kassenprüfer dem Rechnungsführer Entlastung erteilt. Die Mitglieder- 
beiträge wurden auch für das nächste Jahr auf 5.— DM festgesetzt. Zum allgemeinen Bedauern 


lehnte Herr Grammel eineWiederwahl zum Vorsitzenden ab. Es wurden gewählt: Zum Vorsitzenden . 


Herr Walther, zum stellvertretenden Vorsitzenden Herr Meitler und für die aus dem Vorstands- 
rat ausscheidenden Herren L.Föppl und Meitler, die Herren @Grammel und Reh bock. 
Herr Walther sprach sodann Herrn Grammel den Dank dafür aus, daß er die GaMM wieder 


ins Leben gerufen und bisher erfolgreich geleitet habe. Als Orte für die nächsten, Tagungen wurden 


in Aussicht genommen für 1953 Aachen und für 1954 München. Am Schluß wurde dem örtlichen 


Vorsitzenden Herrn Rehbock, herzlicher Dank für die ausgezeichnete konkreitung. der Tagung 
"ausgesprochen. Mit einem Ausflug in den Elm und einem gemeinsamen Alec schloß am 
Samstagabend diese ausgezeichnet verlaufene Tagung. Auch an dieser Stelle sei nochmals vum Namen 
aller Teilnehmer Herrn Rehbock für die große Mühe gedankt, die er auf die Vorbereitung 
dieser schönen Tagung verwandt hat. . Willers. 
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A. Praktische Analysis 


Fehlerabschätzung für die zweite Randwertaufgabe der ebenen 
Potentialtheorie 


Von @. Bertram in Hannover 


Es wird ein Verfahren angegeben, nach dem man für ebene, konform auf den Einheitskreis 
abbildbare Bereiche mit stückweise glatter, doppelpunktfreier Berandung bei bekannter Ab- 
bildungsfunktion fertige Abschätzungsformeln bereitstellen kann, welche die Ermittlung von 
Fehlerschranken für eine harmonische Näherungslösung einer zweiten Randwertaufgabe der 
ebenen Potentialtheorie aus der vorgegebenen Quell-Senken-Belegung einerseits und der Näherung 
andererseits ermöglichen. Als typische Beispiele werden Kreis-, Halbkreis-, Quadrat- und Recht- 
ecksbereiche diskutiert. (Für die erste und bei geeigneten Voraussetzungen für die dritte Rand- 
wertaufgabe erledigt sich die entsprechende Frage unter Verwendung von Randmaximumsätzen.) 

Einzelheiten soll eine später in dieser Zeitschrift erscheinende Arbeit enthalten. 


Interpolation mit Hyperbeln zweiten Grades 
Von U. T. Bödewadt in Brunoy 


Vom Interpolationsfunktionen sind außer den Haupteigenschaften, stetig zu sein und die 
Sollwerte anzunehmen, je nach dem Anwendungsfall noch weitere Eigenschaften zu fordern. 
Rechnerische Einfachheit kommt erst in dritter Linie. Soll die Funktion z.B. positiv und (mög- 
lichst) monoton sein, dann können gebrochene rationale Funktionen (Hyperbeln) den ganzen 
(den Parabeln) vorzuziehen sein. Verf. untersucht die Kurven mit einem Pol erster und einem 
Pol zweiter Ordnung (Isothermen der Yan der- Waals-Gleichung). Sie haben fünf Parameter. 
Der Doppelpol wäre aus 5 Punkten durch eine Gleichung 18. Grades zu bestimmen; nimmt man 
seine Stelle versuchsweise an, so folgt die Stelle des Einfachpoles mit einer quadratischen Glei- 
chung aus 4 Punkten. Die Formeln für die Koeffizienten dieser Gleichung werden angegeben. 
Um der Stetigkeit willen ist darauf zu achten, daß beide Polstellen außerhalb des -Interpolations- 
gebietes fallen. 


Fehlerabschätzungen zum Iterationsverfahren bei linearen 
und nichtlinearen Randwertaufgaben. 


‚Von L. Oollatz in Hannover 


Eine in der Funktionsanalysis für das Iterationsverfahren aufgestellte Fehlerabschätzung 
(J. Weißinger, Z. angew. Math, Mech. 31 (1951), S. 245) liefert bei Verwendung passender 
Normbegriffe übersichtlich verschiedene schon bekannte und zugleich neue Abschätzungen. Bei 
nichtlinearen Gleichungssystemen für Unbekannte &,, ..., 2", (Komponenten eines Vektors &) 
erhält man bei der Norm ||«||—= Max |«,|das Zeilensummenkriterium, bei der Norm ||z|| = Summe, 
der |a,| das Spaltensummenkriterium und bei der gewöhnlichen Länge als Norm das Erhard- 
Schmidtsche Kriterium. Bei Randwertaufgaben ist die in der Funktionalanalysis häufig 
verwendete Norm ||f|| für eine Funktion f(2,, . ... , 2) (Wurzel aus dem über das Grundgebiet B 
erstrecktem Integral von |f|?) für numerische Zwecke wenig geeignet. Die Güte der Fehlerab- 
schätzung hängt sehr von der Wahl der Norm ab. Bei der Norm || {|| = Max | f| (Maximum in B) 
wird der Fehler in B durch eine Konstante abgeschätzt; bei ||/||—= Max |//Z], wo Z eine in B 
(evtl. ohne den Rand) positive Funktion ist, wird die Fehlerschranke proportional Z. Häufig 
erhält man mit Z als einer Eigenfunktion einer zugeordneten Eigenwertaufgabe eine besonders 
einfache Abschätzung; genauer fallen Abschätzungen aus, die die Greensche Funktion benutzen. 
Bei nichtlinearen Differentialgleichungen L[u]= F(&,..., an; u), wobei F bezüglich « 
einer Lipschitzbedingung genüge und L[u] ein linearer homogener Differentialausdruck 
für u(2, ..., 2) ist, ergeben sich unter gewissen Zusatzvoraussetzungen ebenfalls mit den 
obigen Normbegriffen Fehlerabschätzungen. Eine ausführlichere Darstellung soll demnächst in 
dieser Zeitschrift erscheinen. 
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Die Lösung algebraischer und transzendenter Gleichungen 
mit Hilfe von rekursiven Folgen 


Von H. Eltermann in Braunschweig 
In diesem Vortrag werden zwei Verfahren A und B wie folgt beschrieben. 


DasVerfähren A 
Gegeben sei eine algebraische Gleichung p-ten Grades: 


DU 0 Dre a ne (1). 
Dieser ordne man die lineare .Summengleichung: 
Tn-+p == a Int+p—1 + Os In+p—2 —. .. 4 4, %n OT Te Te (2) 
zu. Diese benutze man als Rekursionsformel zur Berechnung einer Folge von Zahlen: Ja Dan RG 
%n, ..., wobei man von p beliebigen Anfangswerten: &,, %,,...., ©, ausgehen darf. 


Das Verfahren stützt sich auf den folgenden aus der Theorie der linearen Differenzen- 
gleichungen bekannten Sachverhalt: 


Hat die algebraische Gl. (1) die verschiedenen Wurzeln: 0, 05, ...., die man sich nach 
fallenden Beträgen geordnet denken kann, so gilt für x, die Darstellung: 

REN a LEN 2 RE U = HS), 

Hierbei sind die C, an sich beliebige, aber durch die Anfangswerte &,, %,,..., ©, wohlbestimmte 


Polynome in n; und zwar ist der Grad des Polynoms C,, (n) gleich (k —1), wenn o, eine k-fache 
Wurzelist. Sind also speziell alle 0, einfache Wurzeln, so sind die C', Konstanten. — Das Letztere 
wird nun vorausgesetzt, wobei aber bemerkt wird, daß das Verfahren keineswegs an diese Voraus- 
setzung gebunden ist. — Im allgemeinen sind die ©, sämtlich von Null verschieden. Wählt man 


als Anfangswerte speziell: =, = -- = %-ı =, z, = 1, so ist dieses sogar stets der Fall. 
Es gibt nun verschiedene Stufen des Verfahrens A: Ay, As, ..., 4,,..., die resp. den 
Voraussetzungen: 
ER ER De 2) 
entsprechen. 
4,) Vor.:|a|>|eı:]. E 


Es gilt die asymptotische Darstellung: ®, = C, 0%. 
Daraus folgt: 


— 01 CR OENEROAEN EEE 


Nn>& v 
Also kann die Wurzel vom maximalen Betrag durch die Folge der Quotienten beliebig genau 
bestimmt werden. 
4,)Vor.: |&|>|o.|>|ee+.|., 
Es gilt die asymptotische Darstellung: ©,= 0, 07 +0,&. 
Ferner gelten die asymptotischen Darstellungen: 


In In+i 
Br In+1 In+2 57 Ci 0, 2 br) (01 0,)" 
und 
‚ %n In+2 
Fr In+ı In-t3 57 10, (0 — 1)” (OL le). rn en 6) 
‚Daraus folgt: 
ed | RR 
na mia et 


Es können also die symmetrischen Funktionen der beiden Wurzeln mit .maximalem Betrage 
beliebig genau bestimmt werden. 
| 4) Vor.: lalzle.| =... 2lal>]e-l- 


ilt di 2 4 ten- für 
Es gilt die asymptotische Darstellung: , I, +G.& +... + C, 0%. Ferner gelten. f 
ie ktonde rashige Determinanten asymptotische Darstellungen, die es gestatten mit Hilfe 


passender Quotienten die symmetrischen Funktionen der maximalen r-Wurzeln beliebig genau 


zu bestimmen. - 
5 16* 
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Zusatz: Das Verfahren ist in jeder Stufe auch dann anwendbar, wenn die in Frage 
kommenden Wurzeln mehrfache Wurzeln sind. Nur muß man dann jede der Wurzeln ihrer Viel- 
fachheit entsprechend zählen. Das Verfahren A, ist allerdings auch dann anwendbar, 


wenn die Wurzel vom maximalen Betrag mehrfach ist. 
Das Verfahren B. 

Dieses ist in gewissem Sinne als invers zum Verfahren A zu bezeichnen. Gegeben ist eine 
algebraische Gleichung: 


1=8 +, 4... ... 0000 (7). 
Man ordne ihr die lineare Summengleichung: 
hm atlma +: une. (8). 
zu und benutze sie wie beim Verfahren A zur Berechnung einerZahlenfolge: &, &,...„ Ins» +» + 
Analog zu (3) gilt die asymptotische Darstellung: ie 
= Got7+ dt Re en (9). 


Es gibt wiederum verschiedene Stufen B,, B,,..., B,,..., die den Voraussetzungen V;: 
alsle|Is---<|o|<|errl für 21 entsprechen. 
Im übrigen wird das Verfahren B analog zu A entwickelt. 


VerallgemeinerungdesVerfahrensB für transzendente 
Gleichungen 


Es gilt der folgende Satz: 

Satz: Es sei f(x) = 0 eine Gleichung. f(x)sei in einem Kreis mit dem Radius R bis auf 
Pole regulär und habe in diesem Kreis eine einzige Nullstelle og. f(0) sei von Null verschieden. 
Entsprechend der Taylor-Entwicklung der Funktion f(x) um den Nullpunkt gilt dann die 
formale Gleichung: 


1 -8 7 a SE a a er An (10). 


zu und berechne eine Zahlenfolge x,, wobei man von endlich vielen Anfangswerten 2, 2% -. - 
ausgeht. Dann gilt für x, die Darstellung: 


we oO LRFN 5er ee (12). 
Also ist dann das Verfahren B, anwendbar. 


Dieser Satz gilt in sinngemäßer Verallgemeinerung auch bei mehreren unter Umständen 
mehrfachen Nullstellen im Kreis mit dem Radius R. 


An dieser Verallgemeinerung des Verfahrens ist besonders interessant, daß die Gl.(10) an 
den betreffenden Nullstellen durchaus keinen Sinn mehr zu haben braucht, da die Ta ylor- 
reihe der Funktion f(x) dort nicht mehr zu konvergieren braucht. 


Verbesserung einer Iterationsfolge bei Gleichungssystemen 
Von Rudolf Ludwig in Braunschweig 


‚Ein Gleichungssystem, das algebraisch oder transzendent sein kann, habe als Vektor- 
gleichung geschrieben die Form eines Iterationsverfahrens: 


Les ALER) N. «.0f a RS Re (1) 


mit den Komponenten %i, nr = fi(@1,% +. +» %n,x). Für die Abweichung n; des k-ten Näherungs- 
vektors von einem Lösungsvektor £, 7, = I; — T, lautet die nach dem ersten Glied abgebrochene 


'Taylorentwicklung: 
NE = NE) gl 2 ee ee (2), 
wobei %: (E) = (fr,(E), +» fz„(E)) = (03), oder | 
BE Oi ee RT (3), 
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dabei ist in dieser Näherungsgleichung £ t durch f ersetzt. Zur Bestimmung des unbekannten Vek- 
tors £ müssen außer ‚dem Ausgangsvektor noch (n -+ 1)iterierte Vektoren bekannt sein. Die u-te 
Komponente dieser Gleichung lautet: 


Dun — Eu = Qu1ılm, ı— ı1i—4)+': Be 1—) - et ERET Bz (4). 
Faßt man nun diese u-ten Komponenten zu einer neuen (rn + 1)-dimensionalen Velibrabienee 
zusammen durch Einführung folgender Vektoren: 


%0 21 — %%o \ dy.o [ 1 \ 
a1) nf Mm L®| u) aa a 
In Tr,n +1 %r, e] \n) 1 r 


so lautet diese: 
Yılıt tan N)aut: a GB ch ER (6). 
Diese Vektorgleichung kann man auch auffassen als ein homogenes Gleichungssystem für die 


(n-+ 1) Unbekannten a1, @2, ...; yu— 1... @un, 1. Damit das System (6) lösbar ist, muß 
die Determinante verschwinden: 


det (31; 32» +: +3 3 HH) =, AN ERS ER (7). 


Zur Ermittlung der Komponenten von £ führt eine geometrische Deutung der Determinanten 
schnell zum Ziel: 

Die n vom Ursprung ausgehenden Vektoren 3, liegen mit dem Differenzenvektor d, in einer 
durch den Ursprung gehenden Hyperebene &,. Da einerseits aber die n Vektoren 3, jeweils mit 
einem der n Vektoren d, in &, liegen, muß andererseits 3, in der durch die d, ausgespannten 
Hyperebene g,liegen. Daher gilt für :9=Aıdı +4... +4,d„. DerVektor 3, liegt außerdem 
in-der durch e und 9, bestimmten zweidimensionalen Ebene g; : 1; = 051 9x + 0x2 e. Das Schnitt- 
gebilde von &, und g,ist aber eine Gerade, für die gelten muß 1,=t;. Dieser Vektorgleichung 
entsprechen aber (n+1)homogene Gleichungen mit den (+2) Unbekannten A,,..., A, 0k1» Ok2- 
Daraus ergibt sich für die gesuchte Komponente %, von: 


% = — 0x. 1/0r2 —= — det (d,; AR ECHh, d„,e)/det(d,, ...) o dx) Ba en (8). 
Durch Umformung der Determinanten ergibt sich schließlich: 
mtr deı 
= det (602.50, 80) . (k) dent" + dee 9 
en be) Eon a 9) 


Wendet man (9) auf den Fall» — 1 an, so erhält man die von Steffensen (1933) !) angegebene 


Formel: 3 
Im A .. (10). 


Wie hier nicht näher erläutert werden kann, entspricht der nach (9) aus x, hervorgegangene 
Vektor £ etwa einer einmaligen Anwendung des Newtonschen Verfahrens. Weiterhin kann 


FE en PER EEE a 0 | ba | 51 |o |5@ | .% 
dy, 1 dyatdıs dgat ds Agst ds Es € 2 dy2 | bye v1) v „2 v3) 
das | 1 da dpa dga || Cs | Cgg | a3 | das 
das | 1 0 0 | 0 er e* 0 | Ole 
—8, > —s(1) 0) Be 0, I |o| o, Io 


di,s= %i,k — Mit, 

b,e= — (die + c1,Rb,ı+...+0—1,% bi,r—1):Ck,k, 

Gk= dir +cı,sbi,ı+. + Gh kbii—, 
di CH 

Wera RI VOR 


Erklärungen: - 


c1,k= di,k, ER Skat &- 


3/Y'r; ZA; Willer 5% ‚Methoden der praktischen Analysis. 2. Aufl. (1950), $ 27,7. S. 258. 
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man nun f als Ausgangsvektor einer neuen Iterationsfolge wählen und dazu den verbesserten 
Vektor L berechnen. Der Vorteil dieser Methode ist, daß man den verbesserten Vektor nur durch 
elementare Rechenoperationen erhält ohne Berechnung neuer Funktionswerte oder Ableitungen. 
Das Ausgangsverfahren kann unter Umständen sogar divergent sein und doch ist die Folge 
to» &  . . konvergent. 

Die Determinanten von (9) unterscheiden sich nur in der letzten Spalte;; ihre Ausrechnung 
ist am einfachsten nach dem Matrizenverfahren von Zurmühl?). Die Tabelle auf Seite 233 
enthält das Rechenschema mit gewissen Summenkontrollen für 3 Unbekannte, das sich 
ohne weiteres auf beliebig viele Unbekannte erweitern läßt. i 


Bericht über ein in Bearbeitung befindliches Tabellenwerk 
für elliptische und verwandte Funktionen 


Von M. Schuler in Göttingen und H. Gebelein in Landshut 


% 
° 
! 
j 
Entsprechend der großen und stets wachsenden Bedeutung der elliptischen Funktionen 
und solcher, die sich darauf zurückführenlassen, gibt es für sie bereits eine erhebliche Menge von ” 
Tabellenwerken (zusammengestellt in dem umfassenden Bericht von Alan Fletcher, 
‚Guide to Tables of Elliptic Functions‘ in,,Mathematical Tables and other Aids to Computation“ j 
Bd. III, 1948). Es handelt sich dabei stets um die Wiedergabe von Funktionen zweier Veränder- w ' 
licher, so daß man sich auf ziemlich weite Intervallschritte beschränken muß, um den Umfang I 
in erträglichen Grenzen zu halten. Diese beiden Veränderlichen sind die Amplitudengröße ® 
und der Modul © der elliptischen Integrale erster und zweiter Gattung, oder sie hängen eng damit ö 
zusammen wie z.B. die Größen k=sin® und k?=sin?©, die ebenfalls häufig verwendet ! 
werden. Alle diese letzten Endes auf Legendre fußenden Tabellen sind insofern unrationell, - 
als die nahezu trivialen Ergebnisse für mäßig große © einen unverhältnismäßig breiten Raum 


einnehmen, die Umgebung der entscheidenden Singularität © — W° aber zu kurz kommt, indem \ 
dort die Differenzen außerordentlich groß werden und Zwischenwerte nur durch sehr mühsame j 
Interpolation einigermaßen genau erhalten werden können. 


Der Praktiker aber braucht Tafeln, mit denen durch lineare oder allenfalls quadratische 
Interpolation bereits eine gleichmäßige und hohe Genauigkeit erzielt werden kann. Um dieser 
Forderung zu genügen, wird bei dem neuen Tabellenwerk ein anderer Weg beschritten, der enger 
an die Theorie von Jacobi anknüpft und bei dem an Stelle des Legendreschen Moduls © 
von dem’ Jaco bischen Parameter qg Gebrauch gemacht wird. Tabuliert werden in erster Linie 
die beiden von den Thetafunktionen hergeleiteten Hifsfunktionen 


d,(x) — 2 g!/*sin & 


2 d;(x) — 1 
G(g%; 2) ee >) pi nz und H(o, 2) —— Bee 2 
mit 2=cos2% für die Bereiche 2= — 1(0,05) +1 (d.h.von 2=—1 bis z= +1 mit 


Intervallschritten 2= 0,05) und g*= 0(0,001) 0,100 bzw. q? = 0(0,002) 0,174. Außerdem 
sind Tafeln für die Logarithmen der Jaco bischen Funktionen, nämlich für 
sn u en u f nu 

lgdn u,1g ven und Ig a mit - = 27 =cos x 
für z— — 1(0,05) +1 und qg= 0(0,01) 0,56 in Bearbeitung. Der Bereich des Moduls erstreckt 
sich auf diese Weise bis © = 89°56’. Diese Funktionen sind dadurch gekennzeichnet, daß sie 
sich, obwohl sie transzendent sind, in den interessierenden Bereichen durchweg beinahe linear 
oder quadratisch verhalten. Dies wird dadurch erzielt, daß in Bereichen mit geringer Veränderung 
der Funktionswerte die Tabellenwerte locker verteilt sind, während ihre Verteilung bei Annähe- 
rung an die Singularität © = 90° immer dichter wird. Bei diesen Tafeln genügt für beide Ver- 
änderliche überall lineare oder quadratische Interpolation, um Zwischenwerte mit voller Genauig- 
keit zu erhalten. 

. Um diese neuen Tabellen gut anwenden zu können, ist es unbedingt erforderlich, den Zu- 
sammenhang zwischen dem Parameter q und dem Modul © bequem numerisch zu beherrschen. 
Auch dies geschieht durch Tabulierung eines fast linearen Zusammenhangs, den man gewinnt, 
wenn man Ig(1 — q) über — Igcos © als unabhängiger Veränderlicher aufträgt. Diese Variable 
—log cos © ist überdies auch geeignet, die vollständigen elliptischen Integrale erster und zweiter 
Gattung, K und E besser wiederzugeben, als dies in Abhängigkeit von © geschieht, wovon erst 
Kaplan Gebrauch gemacht hat. Daher wird vorgeschlagen, zur Vervollständigung des Ta- 


?) R. Zurmühl, Matrizen (1950), $. 251-257, 


Math. Mech. 31 (1951). 8. 356370, 
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bellenwerks auch noch die Größen K und K/E in Abhängigkeit von — lg cos © zu tabulieren. 
Somit sollen dem Werk Tabellen der drei Funktionen 

erh —lg(1—g), Kund K/Efür den Bereich —Ig cos 9 = 0(0,01)3 
beigefügt werden. 

Man erkennt sofort, daß alles, was man bei den geläufigen Aufgaben mit elliptischen und 
verwandten Funktionen braucht, ebenso wie aus den Thetafunktionen auch aus den Werten 
dieser Tafeln auf elementare Weise erhalten werden kann. Ein weiterer Vorzug dieser neuen 
Tafeln besteht darin, daß wegen der Einfachheit der Differenzen auch die Ableitungen der in 
Rede stehenden Funktionen mit hoher Genauigkeit entnommen werden können. Bemerkenswert 
ist auch die Unabhängigkeit dieser neuen Tabellen vom Winkelmaß, in welchem der Modul und 
das Argument angegeben sind, da nur trigonometrische Funktionen dieser Winkel vorkommen. 
Die Genauigkeit dieser Tafeln wird 7—10 Dezimalen betragen, und daher werden sie sich voll 
ausnützen lassen, wenn man sie zusammen mit einer Logarithmentafel und trigonometrischen 
Tafeln derselben Genauigkeit verwendet. - 2 

Das neue Tabellenwerk ist eine rationelle Ergänzung der Logarithmentafeln, wodurch 
diese so erweitert werden, daß eine Rechengrundlage für elliptische und verwandte Funktionen 
entsteht. _ 


. Über die Fehlerabschätzung beim Extrapolationsverfahren von Adams. 


} Von P. Matthieu in Zürich-Berlin 

"Das Problem der Fehlerabschätzung beim Extrapolationsverfahren von Adams wurde in 
letzter Zeit von zahlreichen Autoren bearbeitet. In zwei Arbeiten des Referenten, von denen die 
erste vor einigen Monaten in der Z. angew. Math. Mech. erschienen ist !), während die zweite 
nächstens folgen wird, wird dieses Problem auf einem neuen Wege in Angriff genommen. Dieser 
hat den Vorteil, daß er sich im Gegensatz zu den bisherigen Methoden ohne Unterschied auf die 
sämtlichen Varianten des Verfahrens von Adams anwenden läßt, sowie auch auf die Inte- 
gration von Gleichungen beliebig hoher Ordnung. Im übrigen läßt sich das Verfahren im An-. 
schluß an die erfolgte Integration einer Gleichung praktisch äußerst einfach durchführen. Auch 
ist die Fehlerschranke eine ziemlich enge. Sie beträgt im allgemeinen etwa das Anderthalbfache 
des wirklichen Fehlers. \ 

Der genannte Weg beruht im wesentlichen auf zwei Prinzipien. Das erste betrifft folgendes: 
Es sei z. B. ein Anfangswertproblem einer Gleichung 1. Ordnung in der Form 

LEERE y-fay),  Ya)=y.. re. + 

gegeben. Ferner seien y*(x) mit y*(x,) = Y, eine irgendwie erhaltene Näherungslösung und 
ö(z) = y(z) — y*(x) mit ö(x,) = 0 der zugehörige Fehler. Setzt man dann y(«) = y*(&)+ö(8) 
in (1) ein und entwickelt nach Potenzen der kleinen Größe ö(«), so ergibt sich, wenn man nur 


das Gied erster Ordnung berücksichtigt, für ö(x) die Differentialgleichung 1. Ordnung - 


9-3%.(e y)+f® y*) Er -Ö(&)) BE ME (2). 


Dabei ist noch im einzelnen zu untersuchen, unter welchen Bedingungen die Näherung (2) zu- 
= lässig ist. Aus der ‚‚Fehlergleichung“ (2), in welcher alle Koeffizienten bekannte Funktionen sind, 
° Jäßt sich dann der Fehler ausrechnen oder abschätzen. Das genau gleiche Verfahren, das im 


wesentlichen in einer Störungsrechnung 1. Ordnung besteht, läßt sich auch auf alle anderen Varııı 
anten des Verfahrens von Adams sowie auf Gleichungen höherer Ordnung anwenden, wobei im 


‘ letzteren Falle die Fehlergleichung von höherer Ordnung wird. 


Der wesentliche Koeffizient in (2) ist die ‚„‚infinitesimale Funktion“ f(x, y*) — y*. Sie stellt 


+ sich als Differenz von zwei annähernd gleichen Funktionen dar und bereitet bei der Abschätzung 
die meisten Schwierigkeiten. Das zweite Prinzip, das hier nicht hergeleitet werden kann, zeigt 
E nun, wie die ‚„‚infinitesimale Funktion“ in ganz einfacher Weise aus dem Differenzenschema ab- 


‚geleitet werden kann. Man kann also für diese Funktion mit Leichtigkeit eine Majorante bilden 
und diese in (2) einsetzen. Dann stellt die Lösung der so modifizierten Gleichung auch eine 


2 Majorante für den Fehler dar. Ganz entsprechende Verhältnisse gelten auch für Gleichungen 
_ höherer Ordnung. Man hat also auf diese Weise eine ebenso durchsichtige wie allgemeine Methode 


der Fehlerabschätzung. 


Da die Fehlerschranken recht eng sind, lassen sich ferner interessante Vergleiche ziehen 


über die Arbeitsweise der verschiedenen Varianten. Insbesondere führt eine nähere Untersuchung 
zu. der interessanten Tatsache, daß das Verfahren von Adams bei mittelbarer Extrapolation 
stets rund 13mal genauer ist, als bei unmittelbarer. 


3) Matthieu, P.: Über die Fohlerabschätzung beim Extrapolationsverfahren von Adams. Z. angew. 
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Eine neue Methode zur Bestimmung empirischer Abhängigkeiten 
Von R. Wartmann in Darmstadt 
Will man die empirische Abhängigkeit einer Größe von mehreren Variablen rechnerisch 
erfassen, so steht man vor der Aufgabe, durch eine von Beobachtungsdaten gegebene mehrdimen- 
sionale Punktwolke eine Ausgleichsfläche höherer Ordnung zu legen. Begnügt man sich mit 
einer Ausgleichse bene, so macht man für die empirische Abhängigkeit der Größe x, von den 


n 
Variablen &,, -.-, &, den linearen Ansatz X, = a +2 b;&;; die b; werden dabei in bekannter 


i=1 

Weise nach der Methode der kleinsten Quadrate aus » linearen Gleichungen mit » Unbekannten 
ermittelt. Ist dieser lineare Ansatz unbefriedigend, so kann eine Erweiterung in der Weise vor- 
genommen werden, daß man Quadrate «} und Produkte x; %; usw. in den Ansatz mit hinein- 
nimmt, jedes derartige Erweiterungsglied als eine neue selbständige Variable auffaßt und den 
zugehörigen b-Koeffizienten aus dem entsprechend erweiterten linearen Gleichungssystem berech- 
net. Der zur Aufstellung des Gleichungssystems notwendige Arbeitsaufwand, die Berechnung 
der Größen Q;; = S (x; —3;) (8; —%;), schwillt jedoch dabei so gewaltig an, daß dieser Weg 
bei größerem n (etwa von 10 Variablen an) praktisch ungangbar wird. Mit der Diskriminanz- 
methode wird ein Weg gezeigt, wie man durch Einschalten einer Hilfsgröße trotz im wesentlichen 
linearer Rechnung nichtlineare Abhängigkeiten erfassen kann. j 

Steckt sich die gewöhnliche Korrelationsrechnung das Ziel, die Größe x, unmittelbar aus 
den Variablen x möglichst genau zu berechnen, also $ (x, — X,)* zu einem Minimum zu machen, 
so erstrebt die Diskriminanzmethode zunächst lediglich, aus den » Variablen z; eine Kennzahl 
X, zu berechnen: Man teilt den gegebenen Wertevorrat in einzelne Gruppen mit angenähert 
gleichem x, und legt die Parameter von X, derart fest, daß für Elemente derselben Gruppe der 
Wert von X, möglichst gleich ausfällt, während X, von Gruppe zu Gruppe möglichst verschieden 
sein soll. Denn kann man aufgrund der Variablen x die Kennzahl X, angeben und somit die 
entsprechende Gruppe, so ist dadurch mittelbar auch x, angenähert bekannt. 


i n 
Macht man für X,, wiederum als Einfachstes, einen linearen Ansatz: X, = 3 b;z,, so 
i=1 
ergeben sich die 5b, aus der Bedingung, daß das Verhältnis der ‚‚Binnenklassen-Varianz‘ von X, 
zur „Zwischenklassen-Varianz‘ ein Minimum sein soll, also 


2% = b:b;0:; 6 * Q;; = Binnenklassen-Abweichungsproduktsumme 
= ——_ mi 
06,3 b,5,9;; 

5 


af [2 
Führt man die Differentiation aus, so ergibt sich, daß die 5; gegeben sind durch den zum klein- 
sten Eigenwert gehörenden Eigenvektor des symmetrischen Matrizen-Eigenwertproblems 


(& —A8)r =0. (Alle Eigenwerte A, des Eigenwertproblems sind positiv.) 
Trägt man nun die Gruppenmittelwerte x, über den Kennzahl-Gruppenmittelwerten 


X, = b;&;in einem X,x,-Koordinatensystem auf, so ordnen sich diese Punkte mehr oder 
weniger längs einer Kurve an. Man kann entweder mittels der entsprechenden Ausgleichskurve 
graphisch zu jedem errechneten X, das zugehörige X, ablesen, oder man schließt eine eindimen- 
sionale Korrelationsrechnung X, = f(X,)an, wobei man sich mit dem Ansatz für f nach der 
Gestalt der gefundenen Mittelwertkurve richtet. Ist diese eine Gerade, so ist die Diskriminanz- 
methode mit der linearen Korrelationsrechnung gleichbedeutend; andernfalls aber stellt die 
gewonnene Ausgleichsfläche eine mehrdimensionale Zylinderfläche dar, und diese paßt sich der 
gegebenen Punktwolke naturgemäß weit besser an als eine Ebene. — Es sei außerdem darauf hin- 
gewiesen, daß die Gruppenvarianzen erkennen lassen, ob die Streubandbreite sich mit x, verändert. 

Der Arbeitsaufwand ist nicht wesentlich größer als bei der linearen Korrelationsrechnung, 
abgesehen von der etwas schwierigeren Ermittlung der 5,, da die Anzahl der Produktsummen 
nicht gestiegen ist. — Eine Erweiterung des Kennzahl-Ansatzes sowie sonstige Verbesserungen 
sind im gleichen Maße wie bei jeder anderen Korrelationsrechnung möglich. 


Ein ausführlicher Bericht über den Vortrag: 


L. Cremer: Die notwendigen und hinreichenden Stabilitätskriterien bei Voraus- 
SR y setzung positiver Koeffizienten der charakteristischen Gleichung 
wird in einem der nächsten Hefte dieser Zeitschrift erscheinen. 


Referate der beiden Vorträge: 


Th. Kaluza jr.: Mittlere Aproximation mit linear unabhängigen Funktionen. 
A. Walther: er en Arbeit des Darmstädter Institutes für praktische Mathe- 
mati 


Q; ; = Zwischenklassen-Abweichungsproduktsumme. 


sind nicht eingegangen. 
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Einige Potenzreihen mit Werten der Riemannschen Zetafunktion 
in den Koeffizienten 


Von Otto Emersleben in Berlin 


I. 
Aus der Partialbruchzerlegung der ctg-Funktion erhält man für hal <uls 
I—maetgne 2 Scaler y)i)..2..2.2.222..2..2.M), 
kahl ©“ 


so daß als Koeffizienten die Werte der Riemannschen Zetafunktion positiven geradzahligen 
Arguments auftreten, also Werte, die sich unmittelbar aus der für R(s) >1 — absolut — kon- 


vergenten Reihe &(s) = In ergeben. Analytische Fortsetzung dieser Funktion ‚über die 
n=1 
Konvergenzabszisse hinaus ergibt £ (0) = — 5 und daher für |2|<1: 


ee Ne 
k=0 


r Die Koeffizienten dieser Potenzreihe haben hier also dasselbe Bildungsgesetz für alle vor- 
kommenden — geradzahligen — Potenzen von 2. Bei Schreibweise mit Bernoullischen 
Zahlen in den Koeffizienten ?)ist der Aufbau weniger einfach. Es bedarf auch der Klarstellung, 

welche der verschiedenen Definitionen (Zählweisen) dieser Zahlen gemeint ist. 


Mitunter interessieren Sonderfälle von (2), z.B. für2=—: 


2 
= (Ak 62% 1 
5 ea =R, An = Ent REN 


2 Fe 
oder für z = E 


2 2(o ea | an: 


Diese geschlossene Auswertung erspart numerische Berechnung. — 
Durch Differentiation von n ctgrz gemäß (2) erhält man: 


ey on Nager oder a \-2 Her-nen 3); 
k=0 k=0 


sin?72 sinn2 
allgemein ergibt sich durch n-fache Differentiation (n natürliche Zahl) ein geschlossener Aus- 
druck für _ \ 
::- Der Neenze. 
: k=1 % 


Ba: II. 

Einige der vorstehend und nachfolgend berechneten bzw. durch geschlossene Ausdrücke 
dargestellten Reihen erhält man z.B. bei Berechnung der Gitterenergie von auf einer Kreis- 
linie äquidistant verteilten Ladungen (gleichnamiger oder alternierender Ladungen) °). 

1) Siehe auch P.Stäckel, Eine vergessene Abhandlung Leonhard Eulers über die Summe der rezi- 
proken Quadrate der natürlichen Zahlen. Bibliotheca Mathematica 8; (1907—1908) S.37—54. — L. Euler, 
Opera omnia (1) 14, 8. 156—176, Leipzig 1925, insbesondere 8. 175: 

1 TC SR > 
re eote uU Dar Wr, 
2u2 2% Ex 2 1 ä 


wo 83, &(2r) (r natürliche Zahl) bedeutet. 2 
5 3) Nachteäglicher Zusatz: Verel, auch A. 0. Gelfond, Differenzenrechnung (A. O. Teushonn, 


Meuncnenne koHewnsIx pasHocrei), Moskau und Leningrad 1952, S. 324: 
REURmR 04 San] 3 2]. 
SINTTL p=1 | hi: 2 
©) Wiez,B.v.Mangoldt-Knopp, Einführung in die höhere Mathematik II, 9. Aufl., Nachdruck, 


ipzige 1951, S. 574. I : 
Di Ya, '0.Emersleben, Das elektrostatische Selbstpotential äquidistanter Ladungen auf einer 


Kreislinie. Wird demnächst erscheinen. 


be Rn 
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Dort treten u. a. Reihen auf wie: 


= [2k—1\2(2%k) 
Si, ee (4?) 
ne Den 


in der Konvergenzhalbebene R(s) > 0 der rechtsstehenden Reihe definiert ist und in der Halb- 
ebene R(s) <O durch analytische Fortsetzung. 


Für alle s+1 gilt die Funktionalgleichung 


Es) = (1 Ri) C () Wale Be 
also { 
HOT m Ne a a RN 
- 0) = (1 —2)2(0)=-5 (Sb) 
Daher erhält man wiederum aus der bekannten Partialbruchzerlegung für |2|<1: 
RS TEE 
sin zz BR a, 
1 
Hieraus ergibt sich z. B. für z = 5: 
Scah)_® ee BE 5 
Pe also Ai 2 933 1 8 ( a) 


Um eine Darstellung der gesuchten Reihe (4?) zu erhalten, differentiiere ich die aus (6) erhält- 
liche Gleichung für (r/sin x 2) (2% —1) mal. Es ist, wenn n eine natürliche Zahl bedeutet: 


d"(n/sin setz n! = ee 
ee) = I RR pet I ...(2k—n)z2tn-ı . . (4a). 
3, 
1 : 
leitungen sind in 2 Er sämtlich vorhanden und stetig (im Innern des Konvergenzkreises!); 


Die Funktion (n/sinzz) hat in z= = ein Minimum und ist zu2= — symmetrisch. Ihre Ab- 


daher verschwinden die Ableitungen ungerader Ordnung in 2= 3° Mitn=2x—l1l 


erhält man somit aus (4a) in z = x 


Rt 1 S ec N 2 K—2 x 
Kerr &(2k)z A TEN) 
also für alle x: $ - 
= /2k—1 &(2k) 1 
Den ER 


Neue Separationsfälle der Hamilton-Jacobi’schen Differential- 
gleichung 


Von N. H. Forbat in Mons (Belgien) 


Wir verallgemeinern die Methode von Levi-Civita 6) und gewinnen unmittelbar die 
notwendigen und hinreichenden Bedingungen dafür, daß ein rheonomes mechanisches System, 
dessen Kräfte sich von einem Potential V(q, : '-- q., £) herleiten, durch Separation lösbar sei, 


“) Die Entwicklung von (6) ist gleichwertig mit einer bekannten, Bernoullische Zablen in den 
Koeffizienten enthaltenden Potenzreihe beiK.Knopp, Theorie und Anwendung der unendlichen Reihen, 
4. Aufl., Berlin und Heidelberg 1947, S. 209 Gl. (b). Gl. (6) hat jedoch einen einfacheren Aufbau und läßt ohne 
weiteres erkennen, daß der Einheitskreis der Konvergenzkreis ist. 

5) Weitere Potenzreihen mit Werten der Riemannschen Zetafunktion — und zwar nicht notwendig 

anzzahligen Arguments — finden sich auch bei 0. Emersleben, Das Selbstpotential einer endlichen 
ihe neutraler äquidistanter Punktepaare, Math. Nachr. 6 (1951), 8. 155—170, insbesondere $. 166 und 


0.Emersleben, Das elektrostatische Feld einer Raumladung, Ann, Physik (4) 82 (1927), 8. 713— 774, 
insbesondere S$. 733, Gl. (23r). 


.%)Levi-Civita, Mathem, Annalen 59 (1904), S. 383. 


EN 
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d.h. daß di ' : i 
der Kork ie entsprechende Hamilton-Jacobische Gleichung eine vollständige Lösung 


Wa tt Ans I) = W (t) +W, (4)-+W; (9:) ee ie (9) 


besitze, wobei alle W überdies von n willkürlichen Konstanten 6, 'c„ abhängen. Die Separa- 


bilitätsbedingungen für die Gleichung 


oW oW 
are ne) = o(n = sn) 


sind 


®H oH oH Do oh OH oh «ol ®H oH 0H 


Bqu OR Op Op 00m ap du Sprcn 0m 6 Omen ir WM 
Kt=1-.::n; &-#]) 
und 
Ba OH See oB 
& Bells erGn 


og, ot Op Opndt 


- die identisch in den g, p und terfüllt sein müssen. Es kommen also zu den bereits von Levi- 


Ci vita gefundenen Bedingungen (I) noch unsere, für den rheonomen Fall charakteristische, 
Bedingungen (II) hinzu. 


Man kann dann den Stäckelschen Separationsfall wie folgt verallgemeinern. Es seien 


(tl), wi), Gi=h:n) 


‘(r + 1)? willkürliche Funktionen einer Veränderlichen die nur der Bedingung 


f® pi) 2a gÜ) 
Ai; #0 
pm pw) Sr pin) 


unterworfen sind. Man bezeichne mit Z,(= 4), L,::: Z„ die Minoren der Elemente der ersten 
Zeilein A. Die Hamilton-Jacobische Differentialgleichung 


ER wa TER 
ec taz ut ER ER (11) 
: ist durch Separation lösbar und 
4 w-2|y8 Cr gta [(Faf0+ dt 
ir k=1 i= 


ist eine vollständige Lösung. Wir nennen dies den verallgemeinerten Stäckelschen Fall. 


Im weiteren Verlauf wenden wir die Separabilitätsbedingungen (I) und (II) auf das Problem 


. der Relativbewegung eines freien Massenpunktes im Falle eine Potentials V(x, y, 2) der einge- 


prägten Kraft an. Wenn das ‚System O,,, sich bewegt, sind die Bewegungsgleichungen des 
Massenpunktes in bezug auf O,,, wie man weiß, sehr schwer lösbar, auch dann, wenn V%9,2) 
eine einfache Form hat (z. B. V= mgz). Nun erfordern unsere Bedingungen (I) und (II) wie man 
leicht nachrechnen kann, in diesem Falle, daß 

a) Die Separation möglich sein muß, wenn man das System O,,, als fix betrachtet. 

b) Das System O,,, nur eine Translationsbewegung ausführen darf. 

c) Die Geschwindigkeitskomponente des Ursprungs O bezüglich der &-Achse (bzw. y-Achse, 


z-Achse) konstant sein muß, wenn V von @ (bzw. y, 2) abhängt. 
All diese Bedingungen sind natürlich meistens nicht erfüllt, und dies erklärt die Schwierig- 


keiten der Integration der relativen Bewegungsgleichungen. 
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Konforme Abbildung von Kreisbogenpolygonen; ein Beitrag zum 
Parameterproblem 


Von F. Stallmann in Gießen 
Das Problem der konformen Abbildung von Kreisbogenpolygonen führt auf eine Differential- 
gleichung 2. Ordnung der Fuchh sschen Klasse, die wir in der folgenden Form schreiben wollen: 
n—l n—4 
an EB 5%) var sur 
v=1 »=0 : (1). 


n—1 
8 (x ee e,) 
v-1 


Es ist bekannt, daß der Quotient zweier linear unabhängiger Lösungen dieser Differentialgleichung 
die obere Halbebene auf ein Kreisbogenpolygon abbildet, und daß dabei die n singulären Stellen 
es &9 =, &n—1, oo in die n Ecken des Polygons übergehen, die dort Innenwinkel der Größe 6,77 
einschließen. In dieser Weise läßt sich die Abbildungsfunktion für jedes beliebige Kreisbogen- 
polygon gewinnen, wenn nur die — übrigens reellen — Parameter ö,, e, und b, entsprechend ge- 
wählt werden. Aber bei einer konkreten Abbildungsaufgabe wissen wir natürlich @ priors nicht, 
welche Größe die Parameter e, und b, haben müssen, und diese Schwierigkeit ist eben als das 
Parameterproblem bekannt. Die Schwierigkeit ist hier deshalb auch noch besonders groß, weil 
man mit den üblichen numerischen Integrationsmethoden keinerlei Einblick in die Abhängigkeit 
der Lösungen von den Parametern erhält; man kann so nicht einmal grobe Näherungen für die 
Parameterwerte angeben, die sich dann etwa nach Art des von Herrn Matthie u angegebenen 
Verfahrens !) verbessern ließen. Einen Einblick in die Abhängigkeit der Lösungen von den 
Parametern gestattet dagegen die Methode der asymptotischen Integration. Die Grundgedanken 
dieses Verfahrens sollen hier kurz angegeben werden. 
Wir schreiben die Differentialgleichung in der Form 
Y 


" 


% = 1(2) p8ol2) . Er (1a), 
wo 0?9(x) mit der Größe 
N. n—4 
Sr DAR I (© — c,) 
er) =— en 
4 II («— e,) II (©—e,) 
v1 


identifiziert werden soll, und üben hierauf die Liou ville - Transformation aus: 


&= [Yola)dz, uE) = Ypla) Yle). 


Wir erhalten dann eine Differentialgleichung der Form: 


" 


u 
Fri —ı 0°? + o(E) we El aa u ie Tr a ZDF (2). 
p(£) ist überall regulär analytisch mit Ausnahme der Stellen, die den Punkten = e,und — 6; 
entsprechen, wo es Pole 2. Ordnung besitzt. Sei &, eine solche Stelle, so kann unsere Differential- 
gleichung (2) in der Umgebung dieser Stelle durch eine Hilfs-Differentialgleichung der folgenden 
Form angenähert werden: B 


nn Bam, 
Prim be BE er (3) 
Die Lösungen dieser Differentialgleichung lauten R 


ve)=JE— Zu (ielE— 8); | 

wenn man unter Z,, eine beliebige Zylinderfunktion m-ter Ordnung versteht. Es läßt sich nun 
beweisen ?), daß — unter gewissen Einschränkungen — für große lel die Lösungen «(£) von (2) 
und v(£&) von (3) asymptotisch einander gleich werden. 

Nun ist &(@) ein spezielles Schwarz-Christoffelsches Integral, es bildet also die obere 
Halbebene aufein gewisses Geradenpolygon ab. Diese speziellen Polygone lassensich derartin Recht- 
ecke aufteilen, daß jedes dieser Rechtecke höchstens eine singuläre Stelle von (2) enthält. Für 
Lösungen u(£) von (2) erhält man so für jedes Rechteck einheitliche und jeweils im ganzen Recht- 

ı ” s * ” per “ 
Fe All in = er h i ° rar N en | en . der Schwarz-Christoffelschen Formel“, 

°?) R.E.Langer: „Onthe asymptotie solution“, Transactions Am, Math. Soc. 33 (1931), 8. 23—64, 


” 
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eck gültige Abschätzungen mit Hilfe von Zylinderfunktionen, und auch die Relationen, die 
zwischen den in den einzelnen Rechtecken abgeschätzten Lösungen bestehen müssen, lassen sich 
so bestimmen. Der uns hier interessierende Quotient zweier Lösungen von (2) wird sich also in 
jedem Rechteck näherungsweise wie der Quotient zweier Zylinderfunktionen verhalten. Nun 
hat das Gebiet, auf das unser Rechteck durch den Quotienten zweier Zylinderfunktionen ab- 


e, 62 [4 e3 5 


Bild ı1 


gebildet wird, näherungsweise die Gestalt eines Kreisbogenvierecks, wie man leicht aus dem Ver- 
lauf der Besiselfunktionen im Komplexen feststellt. Das gesuchte Kreisbogenpolygon wird also 
in derselben Weise, wie das Geradenpolygon aus Rechtecken zusammengesetzt ist, aus Gebieten 
zusammengesetzt sein, die annähernd die Gestalt gewisser Kreisbogenvierecke haben. Die 
spezielle Gestalt dieser Kreisbogenvierecke und die Art ihrer Zusammensetzung lassen sich an 
Hand der Abschätzungen leicht bestimmen. Das beigefügte Bild 1 zeigt die Verhältnisse im 
Fall eines Kreisbogenfünfecks. (Eine ausführliche. Darstellung des Verfahrens findet man bei: 
Stallmann ,,‚Zum Parameterproblem der konformen Abbildung von Kreisbogenpolygonen‘‘, 
Mitteilungen d. Math. Inst. d. J.-L.-Hochschule Gießen, Heft 43 (1952).) 


Über die auf dem Mantel eines Kreiszylinders ausgebreiteten 
Potentiale der einfachen und doppelten Belegung *) 


Von H. Schubert in Rostock 


Nach Herleitung Fourierscher Integraldarstellungen für die Potentiale der auf dem 
Mantel eines Kreiszylinders ausgebreiteten einfachen und doppelten Belegung wird unter der 
Voraussetzung von Rotationssymmetrie folgendes Randwertproblem in eine Fredholmsche 
Integralgleichung 2. Art übergeführt: Gesucht ist das Moment. einer Doppelschicht, die nur auf 
einem endlichen Stück des Zylindermantels ausgebreitet ist, das von zwei zur Achse senkrechten 
Ebenen begrenzt wird, und ein Potential erzeugt, das bei Annäherung von innen her auf der be- 
legten Mantelfläche stetige Rand werte annimmt. 


*) Eine ausführliche Darstellung erscheint in den „Mathematischen Nachrichten“, 


Ein ausführlicher Bericht über den Vortrag: 

H. Unger: Zur Praxis der Biorthonormierung von Eigen- und Hauptvektoren 
wird in dieser Zeitschrift erscheinen. 
Ein Berieht über den Vortrag: 

J. Heinhold: Zur Konstruktion involutorischer Kerne 


ist im Archiv der Mathematik, III (1952) S. 15—23 erschienen. 
"Referate der folgenden Vorträge: 


H.Bodemüller: Über aktuelle mathematische Probleme der Geodäsie, 
P.Matthieu: Über praktische Anwendung der Schwarz-Christoffelschen Formeln, 
A. Weigand: Zwei Aufgaben der angewandten Potentialtheorie 


sind nicht eingegangen. & 
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C. Mathematische Maschinen und Instrumente 


Grundgedanken und Entwicklungsstand des Darmstädter 
Rechenautomaten 
Von H.-J. Dreyer in Darmstadt 


Im Institut für Praktische Mathematik der Technischen Hochschule Darmstadt sind Ent- 
wicklungsarbeiten für den Aufbau eines elektronischen Rechenautomaten im Gang. In seinen 
Grundgedanken folgt dieser Rechenautomat der Linie, die H.H. Aikenin Harvard-Compu- 
tation Laboratory, Cambridge/Mass., bei seinem im wesentlichen fertiggestellten Automaten 


- Mark IV eingeschlagen hat. So werden Rechen- und Kommandowerk weitgehend aus Netzwerken 


von Selen- und Germanium-Gleichrichtern und Widerständen aufgebaut; im Speicherwerk tritt 
neben den Haupt-Trommelspeicher auf Magnetophonprinzip ein Pufferspeicher mit magnetostati- 
schen Spulenketten, deren Eisenkerne rechteckige Magnetisierungsschleife haben. Die Zahlen- 
stellung ist dezimal mit 01-Verschlüsselung der Dezimalziffern in 4 Parallelkanälen vorgesehen; 
die Dezimalziffern werden nacheinander (serienmäßig) verarbeitet. Im Gegensatz zu Mark IV 
soll die Kommastellung gleitend sein, um beim Bearbeiten physikalischer und technischer Pro- 
bleme langwierige Größenordnungsüberlegungen zu vermeiden. Jede Zahl z wird also in der Form 
z = + a10# geschrieben, davon wird 1 <a< 10 mit 12 Dezimalen, also 48 Werten in Ol-Ver- 
schlüsselung dargestellt, 0 <b <63 mit 6 Dualen, dazu kommen die beiden Vorzeichen. Für 
mathematische, insbesondere zahlentheoretische Aufgaben kann bei fester Kommastellung die 
Stellenzahl beliebig erhöht werden. 

Die bisherigen Entwicklungsarbeiten betrafen in erster Linie Speicherwerk und Rechen- 
werk. Für den Haupt-Trommelspeicher wurde untersucht, wie viel Rechteckimpulse sich je mm 
Länge auf üblichem Magnetophonband geeigneter Sorte aufbringen lassen. Im Versuchsbetrieb 
konnten bis zu 20 Impulse/mm ohne gegenseitige Störung aufgebracht und kurvengetreu abge- 
nommen werden. Die Frequenzbegrenzung liegt bisher im Schreib- und Lesekopf, und zwar bei 
200 kHz. Auf dem Umfang der Versuchstrommel von 40 cm Durchmesser lassen sich etwa 10 000 
Doppelimpulse je Spur unterbringen Wie eng die Impulsdichte im Speicher an diese Grenzwerte 
endgültig angenähert werden kann, müssen Betriebsversuche zeigen. Von beträchtlichem Ein- 
fluß ist die Sorgfalt bei der Trommelherstellung. Im Versuchsstück war die Oberfläche auf +2 u 
genau rund. 

Die Gleichrichter-Widerstands-Netzwerke wurden bisher für das Addierwerk aufgebaut. 
Die Frequenzgrenze liegt bei 10 kHz. Durch bessere Typen von Selengleichrichtern oder durch 
teilweise Verwendung von Germaniumdioden soll sie auf mindestens 25 KHz hinaufgesetzt werden. 
— Das Beschaffen geeigneter Spulenkerne für den Pufferspeicher war schwierig. Die Arbeits- 
frequenz liegt in der gleichen Größe wie bei den Gleichrichter-Netzwerken. So ergibt sich eine 
Additionszeit von etwa 0,6 ms. — Zur Eingabe von Zahlen und Befehlen ist ein lichtelektrischer 
Lochstreifenabtaster mit elektrostatischer Antriebs- und Bremsvorrichtung für den Lochstreifen 
in Entwicklung. 

Soweit nach dem bisherigen Stand der Arbeiten Schätzungen möglich sind, kann mit einer 


betriebsfähigen Laboratoriumsausführung des Rechenautomaten im Verlauf des Jahres 1953 ge- 
rechnet werden. 


Eine neue Integrieranlage 
Von Th. Erismann in Schaffhausen/Schweiz 


Der Konstruktion und dem Bau der Anlage, die ich Ihnen heute beschreiben will, ging eine 


Reihe technischer und wirtschaftlicher Studien voraus, deren Resultate sich kurz wie folgt zu- 
sammenfassen lassen. 


1. Eine mechanische Integrieranlage kann beim heutigen Stand der Technik nur dann Aus- 
sicht auf Erfolg haben, wenn sie die immanenten Vorteile der mechanischen Integration in vollem 
Maße ausnützt. D. h. sie muß billiger und bei der Verarbeitung empirisch ermittelter Funktionen 
handlicher sein als eine programmgesteuerte digitale Rechenmaschine, gleichzeitig aber die ein- 


fach und schnell arbeitenden elektrischen Integrieranlagen an G igkei ielseitiskei 
ARE RT g gen an Genauigkeit und Vielseitigkeit 


2. Der Bau großer Integrieranlagen von 10 Integratoren aufwärts dürfte auch in Zukunft 


den Ausnahmefall darstellen, da der europäische Bedarf a Ä 
Dies Ange at p n solchen Anlagen naturgemäß be- 


aßte Ausführung wurde also von Anfang an derart entworfen, daß 
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sie auch als kleine Anlage mit einem Minimum an Aufwand rationell betrieben werden kann. 
Allerdings wurden als Parallelentwicklung auch Projekte für Großanlagen bis ins Detail ausge- 
arbeitet, wobei vor allem der schwierigen Verkopplungsfrage sowie einem möglichst geringen 
Raumbedarf besondere Aufmerksamkeit gewidmet wurde. 

Auf Grund dieser Überlegungen wurde die Konstruktion aller Elemente in einheitlicher 
Weise nach einem Baukastensystem ausgerichtet, wobei die vollständig neue_Entwicklung aller 
Einzelteile eine gute Abstimmung ihrer gegenseitigen Eigenschaften ermöglichte. 

Im Interesse der Vereinfachung wurde versucht, möglichst 


wenige, dafür aber möglichst vielseitig verwendbare Elemente Jafegrator 

zu schaffen. Bekanntlich benötigt eine Integrieranlage zur 

Lösung aller in Frage kommenden Aufgaben folgende Rechen- % 
elemente (Bild 1): 


x 
1. Integratoren zur Lieferung von Integralen der Form 


Summentrieb 
[ y- da. 


2. Summentriebe zur Lieferung von Summen der Form J+h- ty 


64°% +06, °'%, wobei c, und c, einstellbare Koeffizienten sind \ 
(davon mindestens einer stufenlos einstellbar). 


3. Funktionstische zum Eingeben bzw. zur Registrierung Funktionstisch 


von Funktionen f(x). 
X x 
Bei der vorliegenden Anlage wurde nun der Schritt gewagt, ee a 
ein Element zu konstruieren, das sowohl als Integrator wie auch rin), N 
als Summentrieb verwendbar sein sollte. Dieses Getriebe, das Ne 
wir als Nebenschlußintegrator bezeichnen, besitzt 3 Eingänge TEE YKeCy Zu JeCONSR 
%, yund zund einen Ausgang, der eine Funktion von der Form ,_ den eyern0 
G' ji y-d&-+-c,:z liefert. Voneinem Summentrieb unterscheidet se 
‚sich ein solches Getriebe funktionell einzig dadurch, daß das It ne 
_ ursprünglich konstante Übersetzungsverhältnis c, nunmehr zu Ba 
einem auch im Betrieb variablen Verhältnis c, - ywird. M.a. W. 
der Ersatz des y-Antriebes durch eine Einstell-Handkurbel verwandelt den Nebenschlußinte- 
grator in einem Summentrieb. Die Anwendungsmöglichkeiten dieses Getriebes gehen nun aber 
über die einfache Summenbildung bei konstantem %, die gewöhnliche Integration bei z= (0 und 
die gelegentliche Kombination beider Rechnungen im Rahmen eines bestimmten Kopplungs- 
planes beträchtlich hinaus. Schließt man nämlich auch den z-Antrieb (den sogenannten Neben- 
schlußantrieb) an den Geber an, welcher den #-Antrieb steuert, so erhält man als Resultat 
SIoo y+ 1: de. SE 
Nun ist c, durch Zahnradübersetzungen einstellbar. Es ist also möglich, bei Integration 
_ über die durch den maschinell bedingten Arbeitsbereich von y gegebenen Grenzen hinauszugehen, 
indem man beim Überschreiten dieser Grenzen geeignete Koeffizienten c, hinzufügt oder abzicht. 
Da die mechanisch bedingten Integrationsfehler sich natürlich nur auf den ersten Term des 
resultierenden Binomes erstrecken, kann dureh Wahl geeigneter Verhältnisse grundsätzlich eine 
beliebige Erhöhung der Integrationsgenauigkeit erzielt werden. Bewegt sich der Integrand stets 


= in der Nähe eines gewissen Mittelwertes, so erhält man diese verbesserte Genauigkeit ohne weitere 


Umstände. - Durchläuft er aber einen wesentlichen Teil des gesamten möglichen Bereiches, so muß 
jeweils beim Erreichen der Grenzen für y die Übersetzung c, umgeschaltet werden. Damit ist die 
. Möglichkeit gegeben, neben der üblichen Integrationsgenauigkeit von der Größenordnung 1%; 

gelegentlich, allerdings mit einem gewissen Zeitverlust, beträchtlich genauer zu arbeiten. Im 
vorliegenden Falle kommen 5 verschiedene Übersetzungen von c, in Betracht, die sichere Ge- 
nauigkeiten bis zur Größenordnung 0,2%, liefern. 


In der konstruktiven Gestaltung weicht der Nebenschlußintegrator von den bisherigen 


Kugelintegratoren ebenfalls wesentlich ab (Bild 2). Es ist eine unangenehme Eigenschaft der 
_ meisten Getriebe mit Reibungsschluß, daß sie beim Übersetzungsverhältnis 0 schlechte Genauig- 


keit und erhöhten Verschleiß aufweisen. Beim Nebenschlußintegrator konnte man nun, da ein 


_ Differentialgetriebe ohnehin erforderlich war, den y-Bereich von vorneherein zwischen 2 positive 
Werte verlegen. Eine durch O gehende Integration ist dann bei Verwendung eines gewissen nega- 
tiven c, möglich. Um trotzdem die Kugeldrehachse in einem genügenden Bereich schwenken zu 
_ können, wurde zwischen den Berührungspunkten der beiden koaxial angeordneten konischen 
- Rollen für An- und Abtrieb ein größerer als der sonst übliche rechte "Winkel vorgesehen. Das für 
den Nulldurchgang verwendete c, wurde derart gewählt, daß beim Übersetzungsverhältnis 0 die 


er - Kugeldrehachse in der Symmetrieebene des Systemes liegt. In dieser Stellung rollt die Kugel auf 


beiden Rollen ideal ab, d. h. die Tangentialebenen an die Kugelin den beiden Berührungspunkten 


 ' ae): x E> va 2 u 6% 
by 5 Kay E ae & 7 
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gehen durch den Schnittpunkt von Kugeldrehachse und Rollenachse. Die Steuerung und An- 
pressung der Kugel erfolgt nicht durch die üblichen getrennten Anpreß- und Steuerrollen, sondern 
durch ein ebenfalls konisches und koaxiales um eine vertikale Achse laufendes und ideal auf der 
Kugel abrollendes Rollenpaar, das von einer geeigneten Tangentenmechanik geschwenkt wird. 
Man kann sich leicht überlegen, daß von den 

4 vorhandenen Rollen nur eine angepreßt wer- 
den muß, um das ganzeSystem unter Spannung | AR 
zubringen. i a 


Iengentenmechanik 


Bild 2a 


Zudem hat die Anordnung noch einen weiteren Vorteil. Bekanntlich besitzt die Kugel 
3 Drehungs-Freiheitsgrade. Jede auf der Kugel laufende Rolle vermindert diese Zahl um einen 
Grad. Die Kugeln der bisherigen Integratoren mit ihren 5 bis 6 unabhängig voneinander ge- 
lagerten Rollen waren also grundsätzlich in ihrer Bewegung überbestimmt und nur durch äußerst 
sorgfältige Herstellung konnte der beim erzwungenen Schlupf entstehende Verschleiß vermieden 
werden. Die koaxialen Rollenpaare der neuen Anordnung sind in dieser Hinsicht wesentlich 
leichter zu behandeln. Schließlich reguliert sich die Anpressung der beiden Arbeitsrollen auto- 
matisch beim Schwenken der Steuerrollen. Bei kleineren Übersetzungsverhältnissen wird die 
Abtriebsrolle, bei großen Übersetzungsverhältnissen die Antriebsrolle stärker angepreßt. 

Bei der Konstruktion der Funktionstische wurde keine der bisher üblichen Methoden für 
die Papierbefestigung übernommen. Statt der Aufspannung auf ein Reißbrett oder auf eine 
Walze fie] die Wahl auf ein perforiertes Papierband beliebiger Länge, das von einer geeigneten 

Walze vor- und rückwärts bewegt werden kann 

(Bild 3). Damit hat man es in der Hand, die 
auszuwertenden Diagramme stark in die Länge 
zu ziehen, wodurch die bekannten Fehler beim 
Abtasten steiler Kurven vermieden werden. 
Ferner ist es möglich, die ebenfalls auf per- 
forierte Papierstreifen aufgenommenen Meß- 
resultate verschiedener Meßinstrumente durch 
Abpausen oder bei entsprechender Papier- 
breite sogar unmittelbar zu verarbeiten. Wie 
bei der Nebenschlußintegration ist auch hier 
eine bedeutende Erhöhung der Genauigkeit da- 
durch möglich, daß die Diagramme nötigen- 
falls beim Erreichen eines der beiden Papier- 
ER as ii SIR ränder geknickt und somit große Ordinaten- 
ERREGER 0. höhen durch zickzackförmige Diagramme 
dargestellt werden können. Auch hier muß 

allerdings ein Eingreifen von Hand beim Erreichen der Endlagen in Kauf genommen werden. 

Die photoelektrische Abtastvorrichtung arbeitet mit moduliertem Licht, was den Tages- 
lichteinfluß ausschaltet und eine Wechselstromverstärkung ermöglicht. Eine geeignete Wahl der 
Zelle und der elektronischen Elemente gestattet die Verwendung verschiedener Registrier- 


ı 


N VE NEE NN an 3 2.55 


tinten, was bei der Verwertung empirischer Meßresultate von Vorteil ist. Es besteht sogar die 


Möglichkeit, das Papierband aus einem Resultattisch in einen Abtasttisch zu führen und auf diese 
Weise gewisse Verzögerungsprobleme zu behandeln. i 
Für ganz kleine Anlagen ist auch eine Ausführung ohne automatische Abtastung mit Hand- 
kurbel und Fahrlupe vorgesehen. 
Einen der wesentlichsten Bestandteile jeder Integrieranlage bilden die elektrischen Fern- 
übertragungen zwischen den Rechenelementen. Auch hier wurde ein Schritt in Neuland unter- 
nommen. An Stelle der bisher üblichen Sels yn- Systeme, die keine Drehmomentverstärkung 


’ 
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ergeben, sowie des induktiven oder kapazitiven Power-Selsyn, das einen beträ i 
Aufwand an Elektronik bedingt, wurde eine äußerst einfache Schritt-F ernübertragung ne 
Ein kleiner kollektorartiger Geber erteilt den 7 Adern der Übertragungsleitung bei der Drehung 
> eine Reihe rythmischer Stromimpulse. Diese werden von einem ganz ähnlich aufgebauten Emp- 
fänger derart verarbeitet, daß an seinen Ausgangsbürsten bei Gleichlauf keine Spannung, bei 
Vor- oder Nachlauf vorzeichenmäßig richtige Spannungen entstehen (Bild 4). Diese steuern das 
Feld des Servomotors, der den Geber zum Gleichlauf zwingt. Die Schrittzahl pro Umdrehung be- 
trägt dank geeigneter Anordnung der Kollektoren 28. Der dauernd unter Strom stehende Anker 
| des Motors wird durch einen Serienwiderstand vor übermäßiger Erhitzung bei langsamem Lauf ge- 
Be schützt. ‚Die relativ kleine Reibung des Gebers, welche beispielsweise für einen Gonella- 
Scheibenintegrator ein unüberwinliches Hindernis wäre, wird vom Kugelgetriebe mit Leichtigkeit 
aufgebracht. Die Antriebe sind mit einem hydraulischen Schwingungsdämpfer zur Vermeidung‘ 
von Pendelungen, einem Zahnradstufengetriebe für die Wahl geeigneter Maßstabfaktoren sowie mit 
je einem Zähler zu einer kompakten Einheit zusammengefaßt. Ein Geber ist in der Lage, eine be- 
i trächtliche Anzahl von Empfängern zu steuern, da nur die kleinen Feldströme geschaltet wer- 
E;; den. Dabei werden verzweigte Kabel bzw. spezielle Zweigstellen benützt. 


Bild 4 


Dank der von Anfang an auf kleine Abmessungen gezüchteten Bauart aller Teile war es mög- 
| lich, die gesamte bisher ausgeführte Versuchsanlage auf 2 Brettern von verhältnismäßig beschei- 
_ dener Größe zu montieren (Bild 5). Dies gestattet eine direkte Verkopplung durch Steckschnüre 
von Element zu Element. Daß eine solche nur für kleine Anlagen in Frage kommt, versteht sich 
Er, von selbst. { 
0 Die neuentwickelten Elemente haben sich durchwegs nach Überwindung der üblichen Kin- 
_ derkrankheiten sehr gut bewährt. Über die Genauigkeit der Integratoren wurde bereits ge- 
_  sprochen. Die Abnützung wurde in zahlreichen Dauerversuchen zum Teil bei stark übersetzten 

_ Drehzahlen untersucht, wobei nicht die geringsten Verschleißspuren an den durch Friktion 
gekoppelten Teilen festgestellt werden konnten. Die Funktionstische erreichen eine Abtastge- 
nauigkeit von wesentlich unter +0,1 mm senkrecht zum Rand der abgetasteten Kurve. Die 
Fernübertragungen, die bei dieser Anlage die Grenze für die Schnelligkeit der Arbeit festlegen, 
gestatten bei voller Genauigkeit den üblichen Zirkeltest in weniger: als 5 Minuten. 

RR Damit darf das gesteckte Ziel als erreicht bezeichnet werden, nämlich die betriebsreife Ent- 
wicklung einer Integrieranlage, deren Leistungen zwar nur in wenigen Punkten über das bisher 
AR Erreichte hinausgehen, die aber mit einem äußerst bescheidenen Aufwand gebaut und betrieben 
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0. Analogierechengeräte mit linearen Potentiometern 

eu Fe RN A“ N Von 4. >. ‚Speiser in Zürich | 
Unter Verwendung von linearen Potentiometern (d.h. von variablen ohmschen Wider- 
n, deren Widerstandswert proportional einem Drehwinkel ist) und von festen ohmschen 

Xechengeräte bauen, deren Genauigkeit je nach der Güte der verwende- 
dnung von 1% liegt. Die Funktionswerte werden durch Span- 
argestellt. Die elementaren Operationen, wie Erzeugung einer 
plikation, Division, werden auf bekannte Art ausgeführt 
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(vgl. etwa W.MeyerzurCapellen, Mathematische Instrumente, 3. Auflage, S. 137). Die 
Ablesung der Resultate erfolgt dabei stets so, daß durch Drehen an einem Potentiometer der Aus- 
schlag eines Messinstrumentes auf Null gebracht wird. Die Resultate werden also nicht auto- 
matisch angezeigt; es handelt sich somit um das elektrische Analogon einer „offenen Getriebe- 
kette“. Die Genauigkeit ist dabei wesentlich besser, als wenn am Meßinstrument abgelesen 
würde, und außerdem ist die Funktion umkehrbar (eine Schaltung, die multipliziert, kann ohne 
Änderung auch zum Dividieren verwendet, werden). 
Eine wichtige Eigenschaft der Geräte ist nun die, daß sie zur Erzeugung von analytisch 
gegebenen oder empirischen Funktionen verwendet werden können. Dazu kann man von folgen- 
der Tatsache Gebrauch machen: Gegeben ein elektrischer Zweipol, der beliebig viele Potentiome- 
ter, die auf einer gemeinsamen Achse sitzen, und beliebig viele ohmsche Widerstände enthält. 
Dann ist der zwischen den Klemmen des Zweipols gemessene Widerstand eine rationale 
Funktion des Drehwinkels der Achse. (Im gegebenen Bereich ist die Funktion nicht negativ, 
enthält keine Pole,und Nullstellen können nur am Rande des Bereichs vorkommen.) Der Verlauf 
dieser Funktion ist bestimmt durch die Konfiguration des Zweipol-Netzwerkes und durch die 
Widerstandswerte seiner Schaltelemente. Durch geeignete Wahl dieser Eigenschaften läßt sich 
eine große Klasse von gegebenen Funktionen — gleichgültig, ob diese eine einfache analytische 
Darstellung besitzen oder nicht — approximieren; die Darstellung ist nicht exakt, der Fehler 
kann aber durch Erhöhung des Aufwandes beliebig klein gemacht werden (abgesehen natürlich 
von der Toleranz der verwendeten Teile). Die Hauptschwierigkeit besteht darin, für die Her- 
stellung einer gegebenen Funktion eine möglichst einfache Anordnung zu finden, und anschließend 
die richtigen Widerstandswerte zu ermitteln. Zu diesem Zweck sind für gewisse Fälle besondere 
Nomogramme erstellt worden, welche die Anpassung auf graphischem Wege ermöglichen. 
Besonders nützlich erweist sich die Ausdehnung dieses Verfahrens auf Funktionen von zwei 
unabhängigen Veränderlichen, da die Darstellung von solchen Beziehungen in Analogiegeräten 
bis jetzt noch nicht befriedigend gelöst ist. Auch hier entstehen rationale Funktionen. Die Viel- 
falt der verwendbaren Anordnungen ist hier sehr groß; es erweist sich als nötig, die einfachsten 
davon zunächst herauszusuchen und in ihrem Verhalten zu katalogisieren. Besonders nützlich 
ist hier die Geometrie der Gewebe, daihre Anwendung es in vielen Fällen gestattet, eine gegebene 


“Funktion von zwei Variablen in eine einfache Beziehung zwischen Funktionen von nur einer 


Variablen überzuführen. 

Auf diese Weise konstruierte Rechengeräte können mit geringen Kosten kurzzeitig ent- 
worfen und gebaut werden undsind sehr robust und betriebssicher ; sie benötigen nur eine einfache, 
nicht stabilisierte Spannungsquelle, unter Umständen sogar nur eine Trockenbatterie. Es sind 
ausgesprochene Einzweck-Geräte, welche nur die mathematische Aufgabe, für die sie konstruiert 
wurden, lösen. Bisher wurden Geräte für die Auflösung gewisser nichtlinearer Gleichungssysteme, 
sowie für die Auswertung von Ausdrücken mit empirischen Funktionen gebaut; weitere befinden 
sich in Vorbereitung. 

Eine ausführlichere Beschreibung der entwickelten Verfahren wird in der ‚Zeitschrift 
für angewandte Mathematik und Physik“ erscheinen. 


Der Programmator 
Von K.Zuse in Neukirchen 


Moderne programmgesteuerte Rechenmaschinen werden durch Befehlsfolgen gesteuert, 
welche an Hand der auszurechnenden Formeln nach bestimmten Regeln gebildet werden müssen. 
Diese zum Teil sehr mühselige Arbeit soll durch den „Programmator‘ weitgehend erleichtert 
werden. Es wird zunächst ein Zusatzgerät für das programmgesteuerte Rechengerät Z 5 gebaut. 
Die Entwicklung fußt auf Vorarbeiten, die bereits im Jahre 1944 in Deutschland im Zusammen- 
hang mit dem Bau des Gerätes Z 4 in Angriff genommen wurden. 


Die Arbeit des Mathematikers geht bis zur expliziten Entwicklung der Formeln in Form 
von „‚Rechenplangleichungen“. Diese werden an der Tastatur des Programmators eingetastet, 
worauf das Gerät selbsttätig die Befehlsfolgen bildet, welche zur Durchführung der Rechnung 
in dem eigentlichen Rechengerät nötig sind. Die Kontrolle der Befehlsfolgen erfolgt durch 
Drucken auf einer Schreibmaschine.: Es ist auch geplant, die Rückbildung der eingetasteten 
Formeln aus dem Befehlsstreifen durch den Programmator vornehmen zu lassen. Es sind ferner 
Einrichtungen vorgesehen, um das Ändern, Ergänzen und Umschlüsseln (Variation der Speicher- 
zellen-Nummern) der Befehlsfolgen zu bewirken; 


; 
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Das programmgesteuerte Rechengerät Z5 


Von AH. Stucken, in Neukirchen - 


Eine Relais-Rechenmaschine wird beschrieben, die in Weiterentwicklung der in Zürich“ 
arbeitenden Z4 im Auftrag der Firma Leitz, Wetzlar durch die Zuse K. G.-gebaut wird. Das 
Gerät arbeitet im Dualsystem mit gleitendem Komma. Die Rechengenauigkeit beträgt 8 De- 
zimalstellen. Die Rechenzeiten betragen für Addition 0,1, Multiplikation 0,4, Division und 
Wurzelziehen 0,75sec. Zahlen werden über Lochstreifen oder Tastatur eingegeben, Resultate 
optisch angezeigt und gedruckt. Rechenbefehle für Haupt- und Unterprogramme erhält die 
Anlage gleichfalls über Lochstreifen. Errechnete Ergebnisse können durch „bedingte Befehle“ 
den Rechenablauf beeinflussen. Die Maschine ist zunächst im Speicherwerk und in der Programm- 

_ fertigung noch nicht voll ausgebaut. Ein Gerät zur automatischen Programmfertigung d.h. 
zur automatischen Umwandlung von Rechenvorschriften in Befehlslochstreifen wird als Er- 
weiterung entwickelt. 


Ein ausführlicher Bericht über die in den Vorträgen: | 
3 H.BillingundW.Hopmann: Die Göttinger elektronischen Rechenmaschinen 
ER FE eK GEung.a2,% 

A. Schlüter: Lösung von Differentialgleichungen mit der G1. 


behandelten Göttinger Maschinen wird in einem der nächsten Hefte dieser Zeitschrift er- 
‚scheinen. : 
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Ein Referat des Vortrages: u 
ni; Th. Fromme: Die Matrizen der Schaltalgebra 
ging nicht ein. FEW 
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D. Wahrscheinlichkeitsrechnung und Statistik 
Zur Verrechnung und Zusammenfassung von Versuchsergebnissen 
nn... mit der Varianzanalyse 
er, ? Von Hans Geidel in Rethmar (Hannover) 


+ pi = ) — r \ Yih 
Es %y:.-%p | Zi . Z; = Pk 
R un TR ee h “ = 
EN] 9 77° 
Bene Unp: Zn rS N 
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;chema erhält man mit der Varianzanalyse nach Abzug der Zeilen- und Spalten- 
als Quadrat des mittleren Fehlers der Mittelwerte bzw. des Mittelwertes 


z ah, + — — = — 
= re BET, CET 
ar —1)P—1) 


Zahlen 2;, zur V reinfachung der Rechenarbeit bei der Be- 
m, selbs nvariant bleibt. — Eine Erweiterung auf 
nzeln mit der Varianzanalyse verrechnet 
ınd die mit der Varianzanalyse berech- 
zwischen den einzelnen Versuchen werden 
e np Meßwerten «,, vorliegen, die sich ein- 
Mech. 29 (1949), 8. 188-190. R 
ar Die 3 17* 
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zeln in ein Schema aus je p Spalten und n Zeilen einordnen lassen, erhält man als Versuchsfehler 
bei Abzug der Spalten- und Zeilenstreuungen und der Streuungen zwischen den Versuchen 


n—1 5% ,, 9-1 
m? = — —— I m} + —— m*?. 
ung — np ung—i 
m, — Fehler der Einzelversuche »—=1,2,..., 9) 
m* — Fehler aus dem Schema der Mittelwerte der Einzelversuche. 


Über eine statistische Theorie zur Beurteilung des Gütegrades 
von Mischungen 


Von K. Stange in Karlsruhe 


Es wird die Frage beantwortet, wie groß das kleinste Teilvolumen einer Mischung (welche 
zwei Stoffe mit gegebenem Volumen- oder Gewichtsverhältnis enthält) sein muß, damit es 
innerhalb gewisser zulässiger Toleranzbereiche hinsichtlich seiner Bestandteile ‚‚homogen“ ist. 
Wesentlich ist die Unterscheidung der Begriffe ‚‚Volumenhäufigkeit“, ‚‚Gewichtshäufigkeit“ 
und ‚‚Teilchenhäufigkeit‘ für die einzelnen Bestandteile in der Mischung. Es wird vorausgesetzt, 
daß für die Teilchenhäufigkeit statistische Gesetzmäßigkeiten gelten, und gezeigt, daß die für 
Stichproben mit gleicher Teilchenzahl hergeleiteten Ergebnisse in erster Näherung auch für 
Stichproben gleichen Volumens oder gleichen Gewichts gelten, so daß der Übergang von den 
Teilchenhäufigkeiten zu den Volumenhäufigkeiten oder den Gewichtshäufigkeiten erlaubt ist. 

Die Ergebnisse werden auf Mischungen mit drei und mehr Komponenten erweitert. — 
Die Arbeit erscheint voraussichtlich in ausführlicher Fassung im Ing.-Archiv. 


Über asymptotische Entwicklungen der Wahrscheinlichkeitsrechnung ') 
Von H. Lenz in München 


> % +. ., %, seien zufällige Größen mit den Mittelwert 2, — 0, den Streuungen 2 —=o 


und den Verteilungsfunktionen F,(z,) = 9, (=) ; ferner sei = 0o?>kk—n und |a,|" <%, 

v »=1 n 

für einr>3. Dann ist die charakteristische Funktion e'*t der zufälligen Größe & = Er > Xp; 
o 


” * . . -1 
wie sich durch Taylorentwicklung (mit Integralrestglied) ergibt, bis auf Glieder der Größen: 


Er She n 
ordnung n!”/? gleich der char. Funktion e'v! der zufälligen Größe en a Yy; wobei die 
we 
y, zufällige Größen sind, deren Verteilungsfunktionen man aus den p, 5 erhält, indem man 
(0) 


%y, 


diese nach der Normalverteilung ® (=) und ihren Ableitungen entwickelt und nach dem 
y % 


r-ten Glied abbricht. Man erhält so die Formel 


oo 


eist — e-UR. > u (d)" +0 (mi—r12) . jejr 


m=0 
of ... on a 
1 n I In 
TER He,, (2) He,, (=) 
nm‘ 0, O7 


Mt +=m Kı 


mit 
AHy< r 


Hen(y) m _ 
m! 


en - - 
oder wegen H, (") Be u Su He 3 —=H (®) — 
% Oy 3 O0, er Oy 
THE Ton Aysı, für R<5 ee Bo FE LAEN 
u. ee 
wi an Pa ER e Kr, 
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!) Ein ausführlicheres Manuskript wurde der Redaktion 1. Mathematischen Annalen eingereicht. 
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Aus dieser-Formel gewinnt man leicht die Abschätzung 


He„(%) | < km . mIm/3]) — m/2 
m! 3 


und daraus asymptotische Entwicklungen für e'*t und (nach Cram&r) für F(x). Wenn die 
char. Funktionen ei%' der x, im Unendlichen mindestens wie t-° mit 6>0 gegen Null streben, 
kann man ohne ‚Verwendung der tiefergreifenden Methoden von Cram&r und Hsu mit 
Hilfe des Fourierschen Umkehrtheorems asymptotische Entwicklungen für F(x) und ihre 
dann existierenden Ableitungen gewinnen, nämlich 


d\r SE Le 
en) F(e) = a e-»l2 . > a Hem+p_1(8) +0 (ne+!—"12) 


mit beliebig kleinem &> 0. 


Moderne Forschungsfragen der Versicherungsmathematik 
Von @ünther Wünsche in München 


Die auf der Jahrestagung der GaMM über das vorbezeichnete Thema gemachten Aus- 
führungen berührten unter ausdrücklicher Betonung des Einflusses der neuesten Entwicklungs- 
tendenzen in der mathematischen Statistik auf die Behandlung versicherungsmathematischer 
Problemstellungen im wesentlichen die folgenden Gebiete und Forschungsbereiche: 


I. Einleitung: Neue Forschungen im Bereich der Versicherungsmathematik als Folge der 
Fortentwicklung mathematisch-statistischer Denkformen und Verfahrensweisen. 
1. Bemerkungen über die Beschränkung auf mathematische Fragen der Personenver- 


sicherung. 


2. Kurze Erörterungen des gegenwärtigen Standes der mathematisch-statistischen Fun- 
dierung und der diesbezüglichen Problematik innerhalb der Nichtlebensversicherungszweige. 


II. Fragen und Untersuchungen im Problemkreis der Entwicklung einer rationellen ver- 
sicherungsmathematischen Rechentechnik. 

1. Methoden zur Ermittlung des Bilanzdeckungskapitals ganzer Versicherungsbestände; 
Hinweise auf die Schwierigkeiten des Problems und den Näherungscharakter der Verfahren. 


2. Anwendungen von Stichprobenverfahren. 

a) Durchführung erster Versuche, insbesondere mit Sequential-Testmethoden; Erörterung 
der Grundzüge der Sequential-Analysis von A.Wald. 

b) Kennzeichnung der mit Hilfe von Sequential-Testverfahren zu behandelnden versiche- 
rungsmathematischen Probleme. 

c) Erfordernis eines ‚‚optimalen‘‘ Sequential-Testplanes. 

3. Weitere versicherungsmathematische Fragestellungen mit ausgeprägt statistischem 
Charakter. i 

4. Das ‚‚versicherungsmathematische Genauigkeitsproblem“ als Konsequenz aus der sta- 
tistischen Natur der Rechnungsgrundlagen. 

a) Wahrscheinlichkeitsstreifen und Näherungsfunktionen. 

b) Beispiele: Schätzung der Bilanzreserven im Zusammenhang mit der F-Methode von 
Jecklin ; Ausgleichung von Beobachtungsergebnissen zur Erstellung von Sterbetafeln. 


5. Nomographische Verfahren. 


III. Schluß: Zusammenfassende Hervorhebung des mathematisch-statistischen Blick- 
punktes bei der Behandlung versicherungsmathematischer Problemstellungen mit dem Ziele 
höchstmöglicher Rationalisierung der resultierenden praktischen Verfahren. 


Referate der Vorträge: Da 


H. Giesekus: Die Anwendung der statistischen Prüfverfahren auf geglättete Zu- 
fallsfolgen, 
U. Graf: Die technische Statistik in der Praxis 


gingen nicht ein. 
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Über Interpolation von Bewegungen 
Von D. T. Bödewadt in Brunoy 


Verlangt wird, von einem Punkte einer im Raum und Zeit gegebenen Bewegung einen 
glatten Übergang zu einer anderen im voraus festgelegten Bewegung zu finden. Bei widerstands- 
loser Bewegung werden Lösungen mit stückweise konstanter Beschleunigung gesucht. Für 
Translationsbewegungen ist die Aufgabe einfach [s. Weltraumfahrt, Frankfurt/Main, 
2, 109—113 (1951)). Übertragung auf Rotationsbewegungen ist ohne weiteres möglich bei 
Drehungen um eine feste Achse, das heißt bei nicht unterstütztem Körper: bei Drehungen um 
eine stabile freie Achse. Für den allgemeinen Fall beliebiger Drehungen ist die Kenntnis dr 
Bewegung unter dem Einfluß fester Drehkräfte erforderlich. Diese wurde für den symmetrischen 
Kreisel vom Verf. angegeben (Math. Zeitschr. 55 (1952), S.310—320), steht hingegen 
für den unsymmetrischen Kreisel noch aus. Es bleibt dann noch ein transzendentes Gleichung- u 
system zu lösen, welches Dauer und Beträge der Drehbeschleunigungen festlegt. d. Ehe, 


Die Bewegungsverhältnisse im Zweistandrollengetriebe 
Von K. Hain in Braunschweig ; 


(Institut für Landtechnische Grundlagenforschung der Forschungsanstalt für Landwirtschaft, 
Braunschweig-Volkenrode.) ) } j 


Als sechsgliedriges Zweistandgetriebe ist im Schrifttum das in Bild 1 dargestellte Gtriebe 
' mit den beiden Gestellpunkten A, und B, bekannt geworden. "Mit diesem Getriebe können : 
gegenüber dem einfachen Gelenkviereck wesentlich höhere Ansprüche, wie z.B. Raten ei 
weiterlaufendem Antriebsglied, erfüllt werden. Bringt man die binären Glieder A’Dund?AC 
zum Schnitt P, so entsteht ein nur für den augenblicklichen Bewegungszustand gültiges Ersatz 
Gelenkviereck B,CPA,!). Ersetzt man in Bild 1 das ternäre Glied BCD durch eine Rolleund 


x I : ee. 
r R =. =. Br: 
a5 D 4 < i 17,0 ee & 


Bild 1 Dr N it Dia a ae, 


ie binären Glieder A’D und AC durch ein die Rolle umschlingen. les Band, so ent 
standrollengetriebe nach Bild 2, dessen Übersetzungsverhältnis zwischen den | 

ell gelagerten Gliedern 4,44’ und B,B durch den Relativpol Q gekennze: 
/Qa, wobei der Punkt Q als Schnittpunkt der Geraden PB mit dem Gestell A 
1. Abwandlungen des Zweiständrollengetriebes entstehen nach Bild 3 
an Stelle der einfachen Rolle eine Differentialrolle oder eine Kurven 
Ä _ Das Zweistandrollengetriebe ist ein periodisches Getriebe, da 
' zwischen An- und Abtriebsglied sich stetig verändert. Es hat g: 
' nach Bild 1 den Vorteil, daß die Herstellung der Einzelteile und der 

stark erleichtert werden; denn es brauchen, um einen leichteren 


hnet- 
B, gef 


n, 


ee . I D % a 

,»)&K.Hain, Geschwindigkeitsermittlung mit Hilfe von Ersatzkurbeltriehe 
Beil age Geiriobetschnik 9 (1941), H.4,S. 177-2178. = = De vr 
_ &K.Hain, Untersuchungen an sechsgliedrigen Kurbeltrieben und d 
u: Präzision 49 (1041), H. 17718, 8.210228, nnd 
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reichen, nicht mehr 7 Lagerbolzen der Drehgelenke genau parallel zueinander zu stehen. Diese 
Getriebeform kann deshalb überall dort Verwendung finden, wo der Kraftschluß durch äußere 
Kräfte immer so vorhanden ist, daß das Band gespannt bleibt. Bandgetriebe sind z.B. bei 
Zahnradfräsmaschinen und bei Waagen als sehr genau arbeitende Getriebe bekannt, so daß das 


u 


Bild 3 Ne 


 Zweistandrollengetriebe Billigkeit in der Herstellung mit der Genauigkeit im Betrieb verbindet. 
Da auch bei Gelenkgetrieben sehr oft Gegenfedern zur Beseitigung des toten Ganges in den Ge- 
_ lenken verwendet werden, ist es ein folgerichtiger Schritt, wenn man diese Gelenke durch 
- ein Wälzband mit Rolle oder Kurve ersetzt. .- x 
Das Zweistandrollengetriebe kann z. B., wie im Vortrag anhand mehrerer Modelle gezeigt 
werden konnte, ohne weiteres als Umlauf- und Schwingrastgetriebe verwendet werden. 


SB | _ Einfachster synthetischer Aufbau der Schwingungsbewegung. 


> ur 


” 


NIAn 


einfachste 'Schwingungsbewegung eines Massenpunkts ist die Kreisbahn vom Halb- 
, die von der Masse m mit der Frequenz » durchlaufen wird. Die Größen m, s und v 
De jewählt werden. Die abgeleiteten Größen sind: 


et 
2 
8: 


el, den die Verbindungsgerade zwischen m und dem Kreismittelpunkt 0 
n Zeiten 7 bildet. Zu einer bestimmten Zeit T ist = vT. Die Geschwindig- 


Eee 


ıhn steht und deshalb keine Arbeit leistet, ist: Kt 
— re Fe j s : } 3 : \ i 


vr u FE PT z 


ee ne er er ar N RTH, { i i 
teekende Energie B: 77. 0... - 


senkrecht zur B 
3/2 ti Et 


ung auf dem Kreis 2xs zerlegen Be. 
n Z, und Z, der Kräfte Zundin 
tung, die durch die dynamische 


are. Ken NE sen 


Kraft Z—= R sinv T = Zy,sin 9, 
T zerlegen kann (Bewegung I). 
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Die Bewegung wird im Beharrungszustand unter der Voraussetzung betrachtet, daß keine 
NT 27 
Dämpfung auftritt. T ist die Umlaufdauer der Masse auf dem Kreis mit T= 2 


Man kann über die Bewegung l eine gradlinige Bewegung v/, = const lagern (Bewegung II), 
zu deren Aufrechterhaltung keine Kraft gehört. Die Bewegungen I und II können beliebig ge- 
wählt werden. Um eine zykloidische Bahn zu er- 


halten, wählen wir 


In Bild 1 ist die zykloidische Bahn dargestellt. Die 
Bewegung wiederholt sich in der Zeit 7, die zum 
Durchlaufen der Strecke 0 bis 4 benötigt wird. Zur 
Bild1. Zykloidische Bahn für sv= v/, bei der Zusammen- Zeit T/2 hat die Geschwindigkeit Bet 
setzung der Bewegung I der Masse m auf dem Kreis2?zs Wert © + sv—_= d oz in X-Richtung un vw 
sr Se ek We NE an Y-Richtung. Zur Zeit T= 0 und T _ T ist 

5 die Geschwindigkeit in X-Richtung null. In Y- 
Richtung tritt die größte Geschwindigkeit -+ d, ZU den Zeiten 7/4 und 3 7/4 auf. 


2. Die beiden Uhren > . 
Wir können eine beliebige Erregerkraft Z=Z, sin» T voraussetzen, ZU der der Zeitmaß- 


stab T = er in Umläufen der Masse m gehört (Uhr a). Die Uhr a läuft mit konstanter Fre- 


el — const um, wobei 9 = 9, 5in (#T) = 9, 5in Qist. Die beiden Faktoren 7 und », aus 


aT 
denen der Wert g unter dem Sinus gebildet ist, können unabhängig voneinander beliebig 
gewählt werden. Wir können z. B. T für die Betrachtung konstant setzen (Uhr a) mit dem 
Ergebnis, daß der Winkel @ verhältnisgleich mit der Frequenz v anwächst. 

Zur Aufstellung der dynamischen Grundgleichung 


a EN 20 


benötigen wir noch eine 2. Uhr b, die die Uhrzeit tin sec angibt. Da die Lage des Massenpunkts m 
durch den Winkel @ bestimmt ist, lesen wir aus der Uhr 5 ebenfalls den Winkel p ab, mit 


Zu = Ze, 608 p = 25608 (WE). Na aan ee le se (7) 


Die Beschleunigung b, von m infolge der Kraft Z, ist d2x/dt? (gemessen in der Uhr b). Die Zeit- 
maßstäbe T und tsind voneinander unabhängig. Wir können aber in Gl. (7) die beiden Faktoren 
v und t (ähnlich wie für 9 = »vT) beliebig wählen. Wenn wir insbesondere für t einen beliebigen 
konstanten Wert annehmen, wächst p verhältnisgleich mit v. Da wir sowohl 7 als auch t als kon- 


stante Werte annehmen können, können wir insbesondere in Sonderfällen 73 —= const setzen. 


Das soll im nachfolgenden für 2 besonders einfache Beispiele durchgeführt werden. 


3a. Die Bewegung von min X-Richtung unter der periodischen Kraft 
Z,=Z,cos»vT 


"Wir betrachten die Bewegung I für sich, setzen also v; =0. Wir beschränken die Be- 
trachtung auf die Bewegung der Masse m in X-Richtung. Zur Zeit T=0 hat die Kraft Z, 
ihren größten Wert Z,,. Sie nimmt bis 7/4 auf den Wert 0 ab, wird negativ, erreicht ihren größten 
negativen Wert — Z, zur Zeit 7/2 oder =. Zur Zeit 3 T/4 hat v, den größten negativen 
Wert —v,,. Zur Zeit T ist der Ausgangszustand wieder erreicht: Mitv, =0istZ =Z,. Esist: 


| sn T =, SÄiDP.: ren RS (8) 
Die Beschleunigung ist %e . &% = 04,9, wobei sich de, mit » ändert. De Weg x aus der 
Nullage ist: dp aT dp 
= %,608vT = 8,005 9; en EN KR ER (9). 
Nach dem Impulssatz ist der Größtwert v,,: 
T/& 
> 2. [vsin SAT = ar... wege #98), 
() Yy 


u u" C 
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vo ist der Größtwert, den die X, Komponente der Geschwindigkeit vsin 9 = sv,sin pzur Zeit 7/4 
erreicht und der gleich dem gegebenen Wert v» =», der Kreisbewegung ist. Aus Gl. (9) folgt 
dr WET (vie n,dnsino..- (10). 

Um die Vorzeichen, die bei der Differentiation auftreten, kümmern wir uns nicht. Der 
Richtungssinn kann in jedem Augenblick durch eine einfache Überlegung festgestellt werden. 
e Wir erhalten eine geradlinige Bewegung um den Koordinatenausgangspunkt O, bei der die 
Größtwerte Z,, und v,, um »/2 phasenversetzt sind. Als Koordinatenanfangspunkt O wählen 
wir für die X-Richtung die Lage von m zur Zeit T/4 oder für den Winkel g =n/2. Der Aus- 


schlag x der Masse m vom Nullpunkt aus folgt aus der Energiebetrachtung unter Berücksichtigung 
von Gl. (9) und (10) 


2/4 Ta 
mv? ın2 n/2 
E= 5 — [2 cos pda = Zu, | cos psinydp=Z,% _ | (Lu) 
ö 0 
Mu Mn, Seren 
ns Z.. en A u N ET (12). 


Die Größen Z,, m und v» = v, sind gegeben. 
So gut wir aber den Zeitbeginn so gelegt haben, daß Z, nach Gl. (8) seinen größten Wert Z,, 


‘ und » den größten Wert x, zur Zeit T = 0 hat, können wir auch den Zeitbeginn um 7/4 ver- 


schieben. 
3b. Die Bewegung von m unter Z,=Z,„ sinvT 


Die erregende Kraft Z, wächst im Zeitintervall O0 bis 7/4 von 0 auf den Größtwert Z,, 
an. Im Zeitintervall 7/4 bis 7/2 nimmt Z, wieder auf den Wert 0 ab. Die größte Geschwindig- 
keit v,, hat m zur Zeit 7/2. Im folgenden Zeitintervall 7/2 bis 7 wird v spiegelbildlich auf den 
Wert 0 abgebremst. Auf einen vollen Umlauf der Uhr a, d.h. auf den Zeitintervall 7 erfährt m 
die Verschiebung 2rs. Wenn man die Messung im Zeitmaßstab # durchführt, ist der zurück- 
gelegte Weg 4s. 

Wenn wir den Bewegungsvorgang im großen und ganzen bei oftmaliger Wiederholung 
der periodischen Bewegung betrachten, ist es gleichgültig, inwelchem Maßstab wir die Betrachtung 
durchführen. Die beiden Zeitmaßstäbe 7 und it sind nur durch einen skalaren Faktor vonein- 
ander verschieden. Für die Aufstellung der Gleichungen ist es gleichgültig, ob die Umlaufzeit 7 
eine sec oder 1/10 sec oder 1/n sec beträgt. 

Wenn wir aber die Bewegung des Massenpunkts x = x, sin vt der Masse m in einem Teil- 
stück 'betrachten, spielt der Unterschied in den Zeitangaben der Uhren a und 5 eine wesentliche 
Rolle. Es ist d7/dt zu berücksichtigen. 

Wir können den Weg, der bei der zusammengesetzten Bewegung lin einer Periode zurück- 
gelegt wird, einerseits in der Uhr a mit der Zeit T messen. Die Uhr a legt den Winkelweg 2 
zurück. Die Masse m legt deshalb auf eine Periode 7’ bzw. 2/» die Entfernung 2 auf der 
Kreisbahn zurück. Wir können aber auch den Weg, den m zurücklegt, in der Uhr b messen, 
weil m nach der dynamischen Grundgleichung in & Richtung beschleunigt wird und hierfür die 
Zeitangaben der Uhr 5 gültig sind. Aus der Uhr b entnehmen wir, daß sich die Periode £ wieder- 
holt, wenn der Winkel tv = (p)$ = 2 durchlaufen ist. Aus der Beschleunigung des Massen- 
punkts, die aus der dynamischen Grundgleichung berechnet werden kann, folgt, daß in einer 
Periode £=2/v, d.h. nach Durcheilen des Winkels @ = 2r der Weg 4 s zurückgelegt wird. 

Welche Uhr hat recht, die Uhr «a oder die Uhr 5? Beide Uhren machen richtige Angaben. 
Es folgt daraus, daß für 2 gegeneinander mit der Geschwindigkeit v, — const bewegte Systeme 
die Längenmaßstäbe im Verhältnis ,/2 voneinander verschieden sind, je nachdem, ob wir die 
Messungen im Zeitmaß 7 oder im Zeitmaß t durchführen. 

Es gibt nur die beiden Zeitmaßstäbe, die von den Uhren « und b angegeben werden. Der 
eine Zeitmaßstab T folgt aus der Überlegung, daß auf eine Masse eine periodische Kraft 
Z, = Z,,605 v T = Z,, cos peinwirkt, und daß wir = »T aus zwei Faktoren von verschiedenen 
Dimensionen zusammensetzen können, von denen der eine (z.B. T) eine beliebige Funktion von & 
sein kann. Der andere Zeitmaßstab t folgt aus der Überlegung, daß sich eine Masse m mit kon- 
stanter Geschwindigkeit v; = const ohne Einwirkung von Kräften gradlinig bewegt, wobei der 
Weg s; = v;t (gemessen in der Uhr b) zurückgelegt wird. Die Uhr a und b beeinflussen sich 
gegenseitignicht. Die Bewegung I hat aberzur Folge, daß die Masse eine Geschwindigkeit vo, = /( T) 
erhält und sich mit dieser Geschwindigkeit relativ zum festen Erdboden bewegt. Die Geschwindig- 
keit v, relativ zur Erdmasse m hängt aber nicht nur von der Komponente der Drehgeschwindig- 


2”8 \p. Neben T tritt der Zeitmafstab t auf, da die 


keit sv sondern auch von v, const = 
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dynamische Grundgleichung die Beziehung zwischen Weggeschwindigkeit und Beschleunigung 
der beiden Massen m und my relativ zueinander angibt. Die Einführung der beiden Zeitmaß- 
stäbe T und t hat deshalb zur Folge, daß die Wege der beiden Massen relativ zueinander im 
Verhältnis z/2 stehen. 


4. Zusammengesetzte Bewegungin X- und Y-Richtung 


Unter 3. haben wir gesehen, daß der Zeitbeginn T wesentlich für den Verlauf der Bewegung 
ist. Imeinen Fall (3a) erhalten wir eine Bewegung um den Nullpunkt O, wobei v, = fb T)=f(p) 
durch eine Cosinusfunktion dargestellt ist. Im andern Fall (3b) wird m zykloidisch um das 


2 
Wegstück 275 (Uhr a) oder um 4 s (Uhr b) im gleichen Zeitintervall T =t = ze fortbewegt. 


Wenn wir aber die dynamische Grundgleichung (4) für die gekoppelte Bewegung in X- und 
Y-Richtung aufstellen, sind wir an einen bestimmten Zeitbeginn T=0, den wir etwa für die 
Einwirkung der KraftZ, auf die Masse m (X-Richtung) angenommen haben, auch bei der Auf- 
stellung der dynamischen Grundgleichung für die Y-Richtung gebunden. Für die eine Richtung 
ist die Beziehung zwischen Z, und der Verschiebung x eine gradlinige Schwankung um den 
Nullpunkt. Für die andere Richtung (Y) liegt eine Verschiebung des Massenpunkts um 2 
vor. Das hat man bisher nicht beachtet. Die Bewegung des Massenpunkts m auf einer Kreis- 
bahn ist eine gekoppelte Bewegung, die man nicht entkoppeln kann. 


Die Bewegung der Masse m unter der Einwirkung der Kraft Z, ist durch die.Zuordnung 
von v, zu jeder Stelle x angegeben. Wir erhalten eine für die Darstellung v, = f(x) zykloidische 
Abhängigkeit, wenn wir für diese Darstellung die Uhr a benützen und wenn wir den Zeitbeginn 
T=0 so legen, daß Z, = Z, sin» T ist. 

Wir können aber die Kurve v, = f(x) auch unter Berücksichtigung der Uhr b darstellen. 
Wenn & = x,sin vt ist, wird die Abhängigkeit zwischen x und 9 nicht durch eine Zykloide 
sondern durch eine Sinuskurve wiedergegeben. Dabei bezieht sich die Kurve auf eine Gerade X, 
die durch den Punkt T = 0 hindurchgeht und deren Extremwerte + v,, und —v,, sind. Die 
Relativität des Zeitmaßstabes hat also zur Folge, daß wir die Bewegung x = f(t) nach der Uhr « 


- als Zykloide und nach der Uhr b als Sinuslinie für die Darstellung v; =f(T) bzw. vu, = ft) 


erhalten. z 


Im nachfolgenden beschreibt der Verfasser ein Anwendungsbeispiel, das zeigt, wie man 
mit Hilfe der vorausgehenden Betrachtungen unter Verwendung der 2 Uhren T und t eine ganz 
einfache Aufgabe aus der Technischen Dynamik lösen kann, die man mit den bisherigen An- 
nahmen der klassischen Dynamik nicht auflösen konnte. 


5. Aufgabe 


Ein Eisenbahnwagen mit der Masse m steht.auf einem ebenen geraden Gleis. An seinem 
Schwerpunkt greift die periodische Kraft Z,=Z,,cos»T in Gleisrichtung an. Es soll der 
Größtausschlag x, von m unter der Voraussetzung berechnet werden, daß die Bewegung von m 
unter der Kraft Z, nicht gedämpft ist. 

Die Betrachtung wird mit der Uhr a durchgeführt. Unter den gegebenen Größen tritt 
nur der Zeitmaßstab T unmittelbar auf. Die erregende Kraft Z, ist von T abhängig. Wir 
können die Berechnung zunächst im Sinne der Bewegung 3a durchführen unter der Voraus- 
setzung, daß zu Beginn der Zeit (T = 0) die periodische Kraft Z, = Z,, cos» T ihren größten 
Wert hat. Wir erhalten eine periodische Schwankung von x, die nach der Zeit 7/2 um «, gegen 
den Ausgangspunkt in X-Richtung fortgeschritten ist, wenn wir den Nullpunkt &=0 zur 
Zeit T/4 ansetzen. Zur Zeit T/2 hat die an m angreifende Kraft Z, den negativen Größtwert — Zz 
mit Größtausschlag — x,. Zur Zeit Tist der Ausgangswert wieder erreicht. Wenn wir zur Zeit T/d 
oder 3 7/4 die erregende Kraft plötzlich wegnehmen, ist die Masse m in Ruhe. Der Größtaus- 
schlag x, der Bewegung folgt aus Gl. (12). 


Mit dem gleichen Recht können wir aber auch den Zeitbeginn nach 3b um r/2 verschieben. 
m legt im der Zeit T a den Weg 4 s gemessen mit dem Zeitmaßstab t zurück. Wenn wir die 
Zuordnung # = f(T) innerhalb des Zeitintervalls T angehen wollen, müssen wir die Abhängig- 
keit der beiden Uhren a und 5 berücksichtigen. 

In beiden Fällen können wir eine bekannte Geschwindigkeit v5 = const über den Be- 


wegungsvorgang lagern und insbesondere den Vorgang 3a betrachten, wobei wir voraussetzen 
daß dem Absolutbetrag nach ist: 


14 
‚ 
3 
j 
£ 
H 
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Bei der Bewegungsumkehr an den Stellen 7/4 bzw. 3 7/4 ist einmal die Geschwindigkeit », 
des Massenpunkts v;, +sv =2v, und das anderemal 1, — sv —=0. 


Nach den Regeln der klassischen. Dynamik sollen die beiden Bewegungen mitsamt ihren 
Kräften superponiert werden können, wobei die dynamische Grundgleichung (6) erfüllt ist. 
Im folgenden Abschnitt wird ein Meßgerät beschrieben, mit dessen Hilfe man die Geschwindig- 
keit v, = const des Fahrzeugs vom Eisenbahnwagen aus bestimmen kann, ohne die Relativ- 
geschwindigkeit abzugreifen. Die Wagenmasse m gegenüber der Erdmasse my soll unendlich 


klein sein. 


. geführt und trägt die Masse m, die infolge der Vorspannung A s der 


‚lagekraft A=cAs von a auf m über- 


6. Das Meßgerät 


In Bild 3 ist eine dünne zylindrische Scheibe « dargestellt, die mit der in der Vertikalen Z 
liegenden Achse Z starr verbunden ist. Durch die vertikale Anordnung der Achse wird die Gra- 
vitation ausgeschaltet. Die Achse läuft in den Lagern c,, «,, die mit dem Fahrzeug (Eisenbahn- 
wagen) starr verbunden sind. An der Achse ist das eine Ende der 
Feder e befestigt. Das andere Ende ist durch einen Schlitz von a 


Feder mit der Kraft e As auf a gepreßt wird (ce — Federkonstante). 
Bei nicht umlaufender Vorrichtung (v = 0) wir die Auf- 


tragen. Man kann die Vorrichtung 
bei ©; =0 und beliebiger Drehfrequenz 
veichen, indem man zu jeder Vor- 
spannung As die Umlauffrequenz v» 
bestimmt, bei der sich m eben vona 
abzuheben beginnt (4 =). - 

Den gleichen Versuch kann man 
bei dem mit der Geschwindigkeit 
v, = const fahrenden Fahrzeug eben- 
falls bei den verschieden Vorspan- 
nungen As durchführen. Eine aus- 
gezeichnete, Frequenz m liegt bei Ti NE a daumen ame 
der Vorspannung As, vor, wenn ten 1und3 ist.ale, Auflgehratt  uhnaenges, rejativ zur fasten Rrde, 
A = ist, d.h. wenn sich m gerade in Mr Z-Achse = Vertikale 
den beiden Augenblicken des Durch- a 
gangs der Masse m durch die X-Achse abhebt (Bild 2). Zum ausgezeichneten Wert », gehört 
der ausgezeichnete Wert der Vorspannung As,. Aus der Eichung des Meßinstruments bei v,— 0 
kann man 


2 


Pe 
w = const 


Bene (14) 


entnehmen, d.h. das zur ausgezeichneten Vorspannung As, gehörige v, bestimmen. Es folgt 


_ die Geschwindigkeit v, = const des Fahrzeugs: 


ae N) ur a a el): 


MYg 


‘Die Größen c, m und die beiden zusammengehörigen Werte As, und »,, die man auf dem 


' Eisenbahnwagen gemessen hat, sind bekannt. Man kann also v, berechnen. 


sind nicht eingegangen. 


Referate der Vorträge: | 
. D. Graffi: Sulle oseillazioni forzate nei sistemi nonlineari. 


A. Signorini: Stereodynamische Anwendungen einer Erweiterung der Cullmann- 
schen Ellipse. 


- 


= 
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Über die mehrdimensionalen Erweiterungen der Airyschen 
Spannungsfunktion 


Von W. Günther in Braunschweig 


Ein n-dimensionaler nichteuklidischer Raum läßt sich im allgemeinen weder auf einen 
n-dimensionalen euklidischen Raum abbilden (d. h.: er ist i. a. nicht konformeuklidisch), noch in 
einen (n + 1)-dimensionalen euklidischen Raum einbetten (d. h.: er ist i. a. nicht von der Klasse 
eins). Beide Eigenschaften kommen allgemein nur den 2-dimensionalen Räumen zu, den gewöhn- 
lichen Flächen. Eine solche Fläche kann stets als ‚„‚Spannungsfläche“ in dem Sinne aufgefaßt 
werden, daß ihre Krümmungsgrößen in direkter Beziehung stehen zum Spannungstensor eines 
ebenen Spannungszustandes. Die Frage nach denjenigen mehrdimensionalen Spannungszustän-- 
den, die nur von einer Spannungsfunktion abhängen, führt auf die Untersuchung solcher ‚,‚Span- 
nungsräume‘“, die entweder konformeuklidisch oder. von der Klasse eins sind. Es sind dies zwei 
Sonderfälle innerhalb einer allgemeinen Geometrie der Spannungsfunktionen, die sich im An- 
schluß an die von der Relativitätstheorie entwickelten Vorstellungen aufbauen läßt. 


Eine Bemerkung zur Theorie der Stabilität erzwungener 
Schwingungen elastischer Körper von Herrn Mettler 


Von Wolfhart Haacke in Braunschweig 


Die Untersuchung der Stabilität erzwungener Schwingungen elastischer Körper führt Herr 
Mettler!) auf die Diskussion der Stabilität der Lösungen folgenden Differentialgleichungs- 
systems zurück: 


@Yyx +02Y + 8005 & I Q,m(®) Ym — 0 GEL Er) a EN 
m=1 - 


Die Stabilität wird in Abhängigkeit von dem Parameterpaar w, e untersucht. Es gilt a,, = a; .. 
Die Konstanten @,, ihrer Bedeutung nach Eigenschwingungsfrequenzen, sind positiv, verschieden 
und seien der Größe nach geordnet (w, > ,+1). Herr Mettler löst dieses System durch einen 
Ansatz mit verallgemeinerten Fourierreihen nach Floquet, deren Koeffizienten dann formal 
nach-dem kleinen Parameter & entwickelt werden. 


Wie an anderer Stelle gezeigt wird ?), ist ein Potenzreihenansatz 


Yan 2 010% @) &? | ...N) 

für “ E 
oo oe—1 
e< E=s =__ e-e(4nn M+1) 
Dar 
absolut konvergent. Dabeiist M — 1, > BL N 
2, max a, ;| 
0<A,SoSA,< + » variiert. Die y,, sollen den Anfangsbedingungen 
vamadn lt (A LEE Dee 

genügen. Sind alle Wurzeln z,, der Stabilitätsgleichung 

det ya) EB (3) 
reell, verschieden und dem Betrage nach kleiner eins, so ist (1) stabil. Die Koeffizienten der al-. 


gebraischen Gleichung (3) sind absolut konvergente Potenzreihen in e für e <E. Mit «= 5 
ergeben sich als Koeffizienten der Potenzreihen der Yar = 


‚ wenn win dem Intervall 


m=1 


0410 = 644 C05 Qu %, ara Zen [in 9 (E— 8) cos Eomavı(E)dE Sri 


1) Mettler, E.: Allgemeine Theorie der Stabilität i i 
AT 1a’ 8, Fe r ® erzwungener Schwingungen elastischer Körper. 
2) Erscheint demnächst in der ‚‚Mathematischen Zeitschrift“, 


RER 
% 


behandelt. 


BE rn _ F Pisstomechenik und Plastizitt = 


‚„‚Instabilitätsbereiche erster Art‘ der w, e-Ebene können in Punkten & _ die &-Achse 


berühren, ‚‚Bereiche zweiter Art“ in Punkten © — [00x 4 = os) Del) 


Im weiteren wollen wir uns auf den Fall N = 4 beschränken, bei dem alles Wesentliche 
zum Ausdruck kommt. Wenn wir (4) in (3) einsetzen und höchstens Glieder bis zu &? berück- 
ou — 0, — w? 


sichtigen, so erhalten wır mit duv == (0! — 0 — @B? — 4 08} 


folgende Stabilitätsgleichungen?) : 


1 1 
)aıprtı h+z. ktz 


9 . 
2, = cos Ang, + en (a,d,,ta,d)+0(e) u=12 u=#r. 
u 
1 
b) 1 te +, „= 5: 
„= —1-— = - + 0 (e°), die andere Wurzel wie bei a). 


rl, +5: 


41577 sin 270g, 


20 Yoo, 
4% sin 2 sin Ang, 


2— eos ing —E7, Er YE1+0(e). 


1=cos2ng + e YF1+0(e), 


2a 20 4 a? 
E a ale 2 0 et 2\ 
u=—1-e 51 |ı 3 2 — y: 3 +00) 


3 
f) 1 tR=2 GT 5> 9% 


r a ee ann (ehrt aut _ DR nn 12 + yD +0(®) mit- 
127 1 A 


4 © 20 — 
» (add + Aladaı @2 dp + et ai, | a EL ) $ ) 
I ©, a 20 — o, 0; mn —w 308 —20 


5 Wi 
Ds rn: 


na, 


a=—14+0@), a=—-1-8 5400). 


1 
hn+ga=!1>2, er5 
4 = cosirqı He (ad + Heda) +0 (E*) , 


2, = cos 2 F e? nanzat (5 + Mc) +0 le). 


1 
NR a 
22 
a=—- 140), a=-—-1-e2 4 +06). 


= ird nicht 
2) Es sei qı > q. angenommen. Der Ausnahmefall g, = (1 _ 3) qı bzw. RT (1 + 1)% wird nic 


4) Die beiden Vorzeichen in D beziehen sich auf die verschiedenen Zeichen von q, +9, = 2. 


u; + 


26 57 se 
rn 


u nr a 


ri 
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Daraus folgern wir: 

Es gilt keine Instabilitätsbereiche erster Art der Ordnung &, der Ordnung & nur bei = 2, 
und ©=2o,. Bereiche erster Art bei ® = — für 1 >1 sind von höherer Ordnung klein. 
Bereiche zweiter Art für Z=1,2 treten nur bei @ +@,—=lw-auf. Für 1=1 ergibt sich bei 
A), + 0 (keine Entkoppelung des Systems) immer ein Bereich zweiter Art der Ordnung e. Bei 
l—2 erhalten wir einen Bereich zweiter Art, wenn D<0 ist. Für hinreichend kleine und 
große a?, ist D>0. Dazwischen gibt es ein Intervall auf der aj,-Achse, das einen Bereich. 
zweiter Art bewirkt, 

1. wenn “4 hinreichend groß oder klein ist, 

Op 1. .: 
2. für alle Paare a,,, Q, genau dann, wenn1l<q,<1-+ zy10 gilt. 


Hier erkennt man, daß auch bei den ‚‚Summenoberfrequenzen“ nicht notwendig Bereiche zweiter 


Art auftreten. : : 
Ist 1>2, so sind die fraglichen Bereiche zweiter Art von höherer Ordnung klein. Gilt 


2 0% 
m 
zusammen, so treten immer Bereiche erster Art mit den ihnen typischen Eigenschaften auf. 


An für ein Paar ganzer Zahlen m, I, fallen also Bereiche erster und zweiter Art 


Über die Drehfrequenzbereiche mit instabilen periodischen 
Torsionsschwingungen bei Kurbelwellen 


Von E. Hölder in Leipzig 


Zu dem Gegenstand, derin C.B.BiezenoundR. Grammel, Technische Mechanik, 
Kap. XII, $9, dargestellt ist, werden begriffliche und mathematische Verschärfungen mit- 
geteilt, die auf zurückliegenden, von mir bei Prof. Blenk in Braunschweig angestellten Über- 
legungen und Rechnungen beruhen. Vor allem wird hervorgehohen, daß Anlaß zu Frequenz- 
bereichen mit instabilen nichtlinearen periodischen Torsionsschwingungen auch die halb- 
ganzen Werte des Periodenquotienten 


. . Frequenz & der Torsionsschwingung 
- Frequenz » der Umdrehung 


geben (genau so wie im Verzweigungsfall die ganzzahligen 7). Die Berechnung solcher Um- 
drehungsfrequenzen » j' ungerade, könnte in der Weise erfolgen, in dr Grammel 


von der Zweitaktmaschine zur Viertaktmaschine übergeht. [Sie sind, wie Herr Grammel 
in der Diskussion mitteilte, zwar theoretisch bekannt, aber nie beobachtet und erwähnt worden. 
Ihr ‚‚Gefährlichkeitsgrad‘, vgl.1.c.$ 7, 32, ist offenbar verschwindend klein, wobei dieser von 
Grammel mit Vorbehalt eingeführte Begriff, ‚die Energie, welche eine durch Störung von 
der periodischen Schwingung abweichende Bewegung in sich hineinfrißt“, wohl noch genauer 
mathematisch untersucht werden sollte]. 

Auszugehen ist von einer exakten periodischen Lösung e® des vollständigen mechanischen 
Problems der wirklichen Maschine — zu konstruieren nach der Theorie der ‚‚Verzweigung perio- 
discher Lösung nichtlinearer Schwingungsgleichungen“ von Iglisch, Math. Ann. Bd. 114 


: (1937) — insbesondere im Verzweigungsfall, den wir aber, da er tiefgreifende Modifikationen 


gegenüber der Betrachtung von Biezeno- Grammel erforderte, hier beiseite lassen wollen. 

Die Theorie bedarf zu dieser Anwendung aber einer weiteren Ausgestaltung, indem 
zum Parameter g(=e), mit dem die Veränderlichkeit der Massen der Kurbelgetriebe sowie 
die äußeren Tangentialmomente multipliziert sind, noch ein Zuwachs p(=4W) der Periode 
einzuführen ist, genauer 

1 
— —— 0 

3 BED; A=%-+p. 
Zu der gewonnenen periodischen Lösung sind dann die homogenen linearen Variations- 
gleichungen aufzustellen, von denen multiplikative (komplexwertige) Null-Lösungen 


Xy) mit X(2n) =uX (0) 


interessieren. 
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‚ Diese Variationsgleichungen unterscheiden sich von der Linearisierung der Bewegungs- 
gleichungen, die Gramme]l vornimmt, s.u. 

Sodann werden zunächst die Stabilitätsgrenzen ermittelt durch das bei w=—l 
entstehende selbstadjungierte Eigenwertproblem, den (nichtlinear eingehenden) Eigenwert p 
als Funktion des Störungsparameters q zu ermitteln. Nach der üblichen, auch von Biezeno- 
Grammelangewandten, von Rellich, Math. Ann. 113 (1936) mathematisch begründeten 
Theorie der Eigenwertstörung erhält man für die Nachbarschaft von q = 0 eine Aufspaltung 
des zweifachen Eigenwertes g=0 in zwei Eigenwerte p=9;(9), i=1,2, aus der Deter- 
minantengleichung 


det! . 2JHRK, Z)dy +qJ HP(R,, X) dy 2: ( 
1+3p[HO (X, X, X)dy +39 [HOX, X, X)dp+--- J 
Bei, 


wobei für die Verdrehung 9°*(y) und den zugehörigen Impuls n%(y) der periodischen Aus- 
gangsbewegung e° nur kurz 


(9° (p)) = PXoly) +9X,p) + 


‚gesetzt ist. 4® sind die Glieder zweiter Dimension, H(? diese Glieder nach g differentiert, 


soweit es explizit vorkommt, A) die Glieder dritter Dimension in der Hamilton-Funktion 
der Kurbelwelle. 

Daß schließlich der die g-Achse in der Mitte enthaltende Winkelraum zwischen diesen 
Grenzen p=p;(g) wirklich ein Drehfrequenzbereich mit — formal und wirklich —insta- 
bilen periodischen Lösungen ist, folgt aus einer Zusatzbemerkung von mir ‚Über die Viel- 
fachheiten gestörter Eigenwerte‘“, Math. Ann. 113 (1936), die sich eben auch funktionentheoretisch, 
hier in einem Vorzeichenwechsel von u? längs p = p,; (q) auswirken. » ; 

[Im Vortrag hatte ich voreilig von doppelt so großen Instabilitätsintervallen gesprochen 
(falls das Differentialgleichungssystem in der kanonischen Form angenommen wird), die überdies 
beim wirklichen Differentialgleichungssystem, wie (*) zeigt, auch von dem — unbekannten — 
nichtlinearen Zusammenhang zwischen Verdrehung und Torsionsspannung abhängen. Herr 
Grammel wies jedoch auf eine von ihm angeregte Untersuchung von Koiter, Proc. Acad. 
Amsterdam (mir noch nicht zugänglich) hin, der das sehr bemerkenswerte und merkwürdige Er- 
gebnis gefunden hat, daß zu Grammels ursprünglicher Rechnung z wei weitere Glieder 
hinzukommen, die sich gegenseitig gerade wieder wegheben. Auf Grund dieses (von mir veri- 
fizierten) Ergebnisses erweist sich die von Biezeno- Gramme] berechnete Aufspaltung 
dennoch als exakt — wenigstens bei Zugrundelegung des linearen Hooke’ schen Zusammen- 
hanges. So lautet die dortige Formel (41, 11) (bei Hinzufügung der Drehmassen ©, des glatten 
Ersatzsystems, die sich ja bei der innomogenen Maschine keineswegs wegheben, sowie bei Aus- 
rechnen der U,(u3) = 7j?w, durch die Amplituden u, der ungestörten Torsionsschwingung) 
unter Verwendung der Amplituden q A,,2;, q B.,2; der Massenveränderlichkeit schließlich 


1 02 0% . 1 3 
„= 3 30 ml2 9w A,2;) + (I 0. PRCT) u) Ulla aRRE: Mr en 


Der elastisch plastische Körper 
Von Hans Jung in Stuttgart 
Wird ein isotroper Körper durch ein Kräftesystem beansprucht, so geht der Punkt pP des 
Körpers mit den Koordinaten %, Y, in den Punkt P mit den Koordinaten @, Y, 2 über. Bezeichnet 
man mit 
u(&, y, 2) die Verschiebung in der x-Richtung, 
v(z, y,z) die Verschiebung in der y-Richtung, 
w(x, y,z) die Verschiebung in der z-Richtung, 


° soist der Zusammenhang der Koordinaten gegeben durch 


== —u(z, y,2), Yey-vıy2, Zeus y2) . - - : (1) 


ae 


260 F. Elastomechanik und Plastizität a nn O0 
BEER EEE EEE GER EEE EEE EIER EESEEEREEREEBBRER GEBE ERS on nn 5 oe sen won 


In bekannter Weise erhält man aus (1) die Verzerrungen 


ne, | 


Eh (2). 
0 „00 RO 0U FE 
mut ey my 0 ay °’ 
Die Verzerrungen (2) sind die Komponenten des Tensors 
1 AR 
s-(} yo), 1, = Ri ee er (3). 


Die Invarianten dieses Tensors werden mit Jj, J,, J, bezeichnet. Der Tensor (3) wird in einen 
Kugeltensor ©, und den Deviator © zerlegt. 
Die am Körperelement im verzerrten Zustand angreifenden Spannungen bilden den Tensor 


0% Toy Tz 
Bi 1 7.00 Tue 1 TER u Se ru 
Tpz Tyz 7 


Der Tensor (4) hat die Invarianten P,, P,, Pzund wird in einen Kugeltensor ®, und einen Deviator 
PB aufgespalten. Mit den Tensorinvarianten aus (3) und (4) lautet das Gesetz von RoS- 


Eichinger 
ST -ı/1 
BysRa—r=rlyays a) N (5). 


Bei Übergang zu kleinen Formänderungen muß das nichtlineare Elastizitätsgesetz in die Be- 
dingungen von Mises 


1 
9% 5 (9 +0, +05) RE RI A 2, Pe 
mit der Nebenbedingung J, = 0 übergehen. Dies wird durch den Ansatz in Tensorform 


? „u OF Re (6) 


y3 Vs P—Pp, ya) rer 


erreicht. Mit (6) und (5) erhält man die Tensorgleichung 


3-84) et). 7), 


wobei die Schubfunktionen 


Ex ER 
el F (r - Vs J—J: ) 
s Y 2 V3 J ng 2 
eingeführt wurde. Mit den Invarianten aus (3) erhält man die Volumdehnung 


U VER I (8). 


Da in (2) nur die quadratischenGlieder berücksichtigt wurden, läßt sich der Zusammenhang zwi- 
schen dem hydrostatischen Spannungszustand und ©, in der Form angeben r 


| : - 
=&f & 7) RR (9). i 

- Mit (7) und (9) erhält man das nichtlineare Elastizitätsgesetz in der Form ; £ 
w 4 

2 1 fl rk 

En (& 3 oa 3r-r.) ta arg): BR 


RK 4/1 | Ber. 
malt: Zu) Fb Dr 


5 ER 
N F 
7% 


schaft. 
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Durch Spezialisierung läßt sich aus (10) das Hookesche Gesetz, das Gesetz von H. Kau- 
derer und die klassische Plastizitätstheorie gewinnen. 


Die durch (10) gegebenen Spannungen müssen noch die Gleichgewich Ing 
N ss ichtsb - 
füllen. Geht man mit dem Ansatz ee ur 


uv=eaWwWw+tdu | 


RE 
| 


vd +@v, rl LESER 

v_5wtrawtagw,::: 
in (2) ein, wobei c,%g, 6, %9, 63%, die 
aus der klassischen Elastizitätstheorie 
erhaltenen Verschiebungen bedeuten, 
so lassen sich mit (10) die Span- 
nungen bestimmen. ı Die : Gleichge- 
wichtsbedingungen führen dann auf 
lineare Gleichungssysteme zur Be- 
stimmung der %,, v;, und w,. Mit den- 
selben Vereinfachungen, die auch in 
der technischen Biegelehre gemacht 
wurden, läßt sich für den durch das 
Bild gegebenen Belastungsfall die Dif- 
ferentialgleichung für die Durchbie- 
gung w aufstellen 


Rn u 
aha a 


wobei M(x) das auf die Platte ausgeübte Biegemoment bedeutet und A, sowie B Materialkon- 
stanten sind. 


Die Beugung elastischer Wellen an der Halbebene 
Von A.-W. Maue in Freiburg (Br.) 


Das Problem der Beugung elastischer Wellen an der spannungsfreien Halbebene, das für 
die Diskussion der Spannungsverhältnisse in der Umgebung einer Bruchfläche in elastischem 
Material von Bedeutung ist, wird nach einer von Clemmow und von Hönl angegebenen 
Methode mit Hilfe zweier Integralgleichungen formuliert. Diese lassen sich auf die“ Aufgabe 
reduzieren, eine vorgegebene verzweigte Funktion so in zwei Faktoren aufzuspalten, daß sich 
die Verzweigungspunkte in: vorgeschriebener Weise auf beide Faktoren verteilen. Die Auf- 
spaltung gelingt mit Hilfe des Cauchyschen Satzes. Eine ausführliche Darstellung des Gegen- 

 standes erscheint demnächst in dieser Zeitschrift. 


Strukturmodelle und Systematik der Fließvorgänge‘) 
Von Ernst Mewes in Braunschweig 


Die Fließeigenschaften vieler Stoffe werden an Hand von mechanischen Modellen veran- 
schaulicht, die aus Federn und Dämpfern aufgebaut sind. Sofern Fließgrenzen auftreten, sind 
Glieder mit fester Reibung hinzuzunehmen. Das allgemeinste Modell, das sich aus parallel ge- 
richteten Federn und Dämpfern mit linearen Gesetzmäßigkeiten aufbauen läßt, ergibt sich durch 
Hintereinanderschaltung von einer Feder und einem Dämpfer, d.i. zusammen ein sog. Max- 
well- Glied, und beliebig vielen Paaren parallel geschalteter Federn und Dämpfer, d.s. sog. 


_ Kelvin-Glieder. Dynamisch gleichwertig hiermit sind Modelle, die aus einer entsprechenden 


Zahl vonMaxwell-Gliedern in Parallelschaltung aufgebaut sind. 
Für beide Modelle ergibt sich die allgemein gültige Differentjalgleichung zwischen Verfor- 
mung y und Belastung z: 


Y5 a’ ++ ty Tr the tb tr Hm 


mit den zeitunabhängigen Beiwerten a und b. nist die AnzahlderMaxwell- und Kelvin- 
Glieder. Diese Beziehung sieht etwas anders aus als die bisher im Schrifttum dafür angegebenen 


+ Aus dem Institut für Landtechnische Grundlagenforschung in der Forschungsanstalt für Landwirt- 
| 18 
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Beziehungen (vgl. M. Reiner, Twelve Leetures on Theoretical Rheology. ae 
hlahing Co., Amsterdam 1949; ders., Deformation and Flow. Lewis, London 2 Eins 
ZDAATTr 8% ‚ Mechanical Behaviour of High Polymers. Interscience Publishers, New r Br \ 
auch: The Methods of Speeifying the Properties of Viscoelastic Materials. J. Appl. Phys. 16 ( ); 
S. 700). ? 
j ichteten Federn und Dämpfern in 
Es wird nun festgestellt, daß auch jedes aus parallel geric = 
beliebig. komplizierter Weise aufgebaute Modell nur ein Sonderfall des geschilderten allgemein 


Wa 


(2 0 4ı 
Ct, 4 


III = 
II = 


@2 (7 du 


ten Modelles ist. Dazu wurde jeweils die Differentialgleichung unter Zuhilfenahme der Opera- 
toren-Methode nach He aviside aufgestellt. Für die Modelle im Bild ergibt sich 


47 +,y=-r+bi+brt 


mit (c = Federhärten, u = Viskositäten) 


a = Mat U = u, = u + ln 
1 .) 
%,=(% +%)b; — 1/0 = hür iR Rn r 
YHı , & , KM _Mı, Kı 
bi = [Wa(65 +6) + Mp(Ce + Ca) + HelCa+ &)\/do = ” Bun 4% = ee ne 
_ Kıkla u Brenn 
b; IA (Ha U nr Up We Ir WeHa)/ bo — 6% Cr rk 


bes 95-46: 


Die Aufstellung von komplizierten Modellen in der Art des Systems im Bild links (s. J.M. 
Burgers, Mechanical considerations concerning viscosity and plastieity. Aus: First Report 
on Viscosity and Plasticity. Verhandl. d. Kon. Akad. v. Wetensch. Amsterdam, Afd. Naturk. 
(1. Sect.) DI XV Nr. 3 (1935)) gibt nichts Neues hinsichtlich des Fließvorganges. 

Durch geeignete Zusammenstellung von Spezialfällen des systematisch aufgebauten allge- 
meinsten Falles (s. 0.) und Betrachtung der Eigenschaften dieser Modelle (Auftreten von elasti- 
schen und plastischen Verformungen, Kriechen, Nacheilen, Nachwirken, Erschlaffen, endliche 
oder verschwindende Fließfestigkeit u. a.) kann eine bessere Systematik der Fließvorgänge auf- 
gestellt werden, als sie in Nature 149 (1942), S. 702, bzw. Schweizer Archiv 9(1943), S. 98 und 
Z. Kunststoffe 34 (1944), S. 26 angegeben ist. 


Selbstangefachte Schwingungen von Maschinenwellen 
Von Eduard Pestel in Hannover 


Von den verschiedenen Arten selbstangefachter Schwingungen von Maschinenwellen wird 
der Fall behandelt, bei dem die Selbsterregung infolge der Preßsitzreibung zwischen Welle und 
aufgeschrumpfter Scheibe eintritt. Der Entstehungsmechanismus möge am einfachen Beispiel 
der vollkommen ausgewuchteten senkrechten Welle erörtert werden, die eine zentrisch aufge- 
schrumpfte Kreisringscheibe trägt, deren Symmetrieebene auch bei Biegung der Welle in ihrer 
ursprünglichen Ebene bleibt. Ihre Drehgeschwindigkeit 2 sei höher als die der 1. Wellenkriti- 


schen. Durch irgendeine Störung werde nun die Scheibe in eine Zirkular-polarisierte Schwingung 


versetzt, deren Richtung mit der Drehrichtung der Welle zusammenfalle. Vom Standpunkt eines 
Beobachters in Scheibenmitte dreht sich nun die Scheibe nur mit einer Drehzahl entsprechend 


der Differenz zwischen Wellendrehzahl 2x2 und Schwingungsfrequenz 2rw. Er stellt nunfest, 
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daß ein Punkt auf der Welle, der infolge der Drehbewegung der Welle z. B. zur Zeit i gerade die 


Stelle der größten Zugbeanspruchung passiert, zur Zeit t-+ a die Stelle der größten Druck- 


beanspruchung durchläuft. Infolge der mit der Zugbeanspruchung verbundenen Längung und 
der den Druckbeanspruchungen entsprechenden Kürzung der Werkstoffasern wird in dem 
Schrumpfsitz Reibungsarbeit geleistet und somit mechanische Energie nichtumkehrbar in Wärme 
umgewandelt. Versuche, von Kimball') haben ergeben, daß bei sonst gleichen Verhältnissen 
die Verzehrung mechanischer Engergie hier etwa proportional dem Quadrat der Durchbiegung e 
der Welle an der Preßsitzstelle ist. Sein Schluß, daß dieser Energieverlust E auch unabhängig 
von der Drehgeschwindigkeit 2 der Welle ist, erweist sich experimentell wie mathematisch als 
unrichtig. Vielmehr ist E unabhängig von der Größe der Relativdrehgeschwindigkeit Q—o. 
Damit verzehrt die Reibung im Preßsitz während einer vollen Relativdrehung der Welle folgen- 
den Betrag mechanischer Energie 


B=2nße. 


Die Beziehung zwischen der Konstanten ß und dem zusätzlichen äußeren Antriebsmoment M, 

das aufgebracht werden muß, um trotz des Energieverlustes E die Drehgeschwindigkeit 2 der 

Welle konstant zu halten, ergibt sich wie folgt. Die sekundliche Leistung von M beträgt M-Q. 
27 


. . Die Zeitdauer der Relativdrehung ist 7 = — —. 


O—ai 
Während der Zeitdauer T leistet also das zusätzliche Antriebsmoment M die Arbeit 


A=M-9-T=-M ng. Wegen E = A folgt 


M=-pete. 


Indem man das Wirken einer geschwindigkeitsproportionalen Dämpfungskraft (angreifend in 
Scheibenmitte) mit der Dämpfungskonstanten b voraussetzt, läßt sich folgendes einfaches Sta- 
bilitätskriterium aufstellen 

2—o 


a EIER UN 


Man erkennt, daß für verschwindende Dämpfung stabiles Verhalten nur bei unterkritischem 
Lauf, also für Q< & vorhanden ist. Dabei ist die Kreisfrequenz w der Störschwingung praktisch 
gleich der der ersten Wellenkritischen. 


Unter Benutzung der Versuchsergebnisse von Kimball gewinnt man durch eine ein- 
fache Modellbetrachtung eine rohe Abschätzung für die Größenordnung von ß 


d? B? re 
ß = 6000 7a ; 
Darin bedeutet B die Preßsitzbreite, d der Wellendurchmesser und / die Wellenlänge. 
Als wirkungsvollste Gegenmaßnahme zur Bekämpfung derartiger Schwingungen hat sich 

der Einbau gefederter Lagerschalen erwiesen. Mit den Bezeichnungen 

c, Federzahl der Welle, 

c, Federzahl der Lagerfederung, 

m Auf die Schrumpfsitzstelle reduzierte Masse des Wellen-Scheibenaggregates 


ergibt sich dann folgendes Stabilitätskriterium, das für c,—co wieder in (1) übergeht. 


e 


nn IE EEE) 


= \ % und og le 1 
G+6 m 


| 1) Kimball, A.L., Measurement of Internal Friction in a revolving Deflected Shaft. General Electric 
Review 28 (1925), S. 554—558. 
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t i Be Z. angew. Math. Mech. 
264 - F. Elastomechanik nnd Plastizität Bd. 32 Nr. 8/9, Aug./Sept. 1952 


Zur Festigkeit im vielkristallinen Werkstoff 
Von A. Slibar in Wien 
wird angenommen, daß die metallischen Werkstoffe sowohl homogen 
noch die zweite Annahme ist auch nur ser en 
ss li 2 i Grenzproblemen, die bisher mit der Spannungstheorie des Kon inu- 
a Een AR race des Versagens im realen Aufbau der Werkstoffe sucht. Die 
rechnerische Berücksichtigung des tatsächlichen Aufbaues und der Glesteigenschaiten des Poly- 
kristalls erlauben nur die Methoden der Statistik. Gemeinsam mit F. Vito vec!) wurden vom 
Verfasser die Einschnürung einer zentralsymmetrischen Zugprobe behandelt. Bei Berücksichti- 
gung der tatsächlichen Verformungsmöglichkeiten des Werkstoffes folgt für die Kontur des ein- 
geschnürten Stabes die Verteilungsfunktion der Einzelgleitungen 
N ERRTTRE 
p(a) == "€ 
7 
deines Querschnittes von Stabmitte, q das vom Werkstoff, dessen Vorbehand- 
ne Ed aan ae abhängige Einschnürmaß bedeuten. Die entwickelte Funktion 
erlaubt die analytische Erfassung des gesamten Konturverlaufes. > 
In einem vielkristallinen Werkstoff behindern sich die Kristallite gegenseitig an der Ver- 
formung. Die Behinderungsbeiträge für ein bestimmtes Bezugskorn hängen ab A 
1) von der Anzahli der zwischen den beiden aufeinander einwirkenden Kristalliten liegenden 
Schichten, 
2) von der Zahl der möglichen Gleitkombinationen des Kristallsystems (K) und der mitt- 
leren Orientierungsabweichung der Kristallachsen dieser beiden Kristallite zueinander (Aw) P 
3) vom Einfluß der Korngröße und ;Korn- 
oberfläche (x). 

Der Ansatz für die Wahrscheinlichkeit des Be- 
hinderns eines Bezugskornes durch eine Kugel- 
schicht mit dem Durchmesser 3-3 (3 mittlerer 
Korndurchmesser) ergibt sich daher zu 


Aw 
Ve ne 3 . (2) 


K 


n, ist in (2) die Anzahl der Kristallite in der als 
Kollektiv mit vorläufig unendlicher Merkmal- 
menge zusammengefaßten behindernden Schicht. 
Die gesamte Stützwirkung, auf ein Korn folgt 
als Summe der Einzelbeträge 


In der Festigkeitslehre 
als auch isotrop sind. Weder die erste 


v()=x- 


x —- 2 K-Aw+.it.x 
Wi) = [. 405218) 
rat 7 = 140) 
Aomzahl ı bıs Oberflache DE E 5 R a 
Verlauf der Verformungsbehnderung ın Abhängigkeit rom wenn K, Aw, i und x die bereits eingeführten 
Kornabstand von der Oberfläche Größen, Av das umgebende Volumelement und 


V («) den gesamten Behinderungsbereich als Funktion der Anzahl der umgebenden Schichten 
darstellen. Die rasche Konvergenz von (2) bzw. (3) erlaubt ein Abbrechen der Summenbildung 
bei einer Grenze «* = 10 bis 15 Kristallitschichten. Bedeutet 3 den mittleren Durchmesser eines 
Kristalliten, so ergibt die Größe ö* - 3 die von L. Föppl1?)nach Versuchen vorgeschlagene neu 
einzuführende dritte Materialkonstante. Der Verlauf der Elastizitätsgrenze über den Querschnitt 
folgt bei Berücksichtigung der Störung des Behinderungsbereiches bei Annäherung an die Ober- 
fläche (Bild). | 
Weitere durch Oberflächeneinflüsse verursachte Erscheinungen wurden von den Verfassern 
behandelt und deren Ergebnisse werden in nächster Zeit zur Veröffentlichung gelangen. 


.. )F.VitovecundA.Slibar, Bestimmung der Form der beim Zugversuch sich au bildend in- 
BEIDE: Öst. Ing. Zub IV (1950). ys: 5 N 
F. Vitovec und A. Slibar, Die Einschnürung einer Zugprobe als statistisches Probl $ 
Vierteljahresschrift d. Univ. Wien, Bd. 4 (1951). < ji ee 


°?) A. Slibar und F. Vitovec, Über den Wirkungsbereich der $ sti 
Arch. 16 (1950). er gsbereich der Bpannungsverfestigung. Schw. 


F.VitovecundA.Slibar, Über den Einfluß der freien Oberfläche auf die Verformungsbehinderung 


im a N Schw. Sri 16 (1950). Er 
.SlibarundF. Vitovec, Über den Einfluß der freien Oberfläche auf di izitätsgrenze viel- 
kristalliner Werkstoffe. Federhofer-Girkmann Festschrift, Wien 1950. Rn ; ee Hr 
°»)L.Föppl, Drang und Zwang. III. Bd., Leibniz Verl., München 1947, 8. 178. 
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Hydrodynamische Theorie fester Stoffe 
Von ©. Torre in Wien 


Die geologischen Erscheinungen werden nach den Gesichtspunkten der Werkstoffmechanik 
gruppiert und mit Hilfe dieser die Spannungen in einem Gebirge berechnet. Wir haben bisher 
insbesondere das klüftige Gebirge (Fels) eingehend mathematisch untersucht. Dabei wurden 
jene Kräfte berechnet, die bei der Entstehung der Klüfte geherrscht haben. Da in einem schon 
beruhigten Gebirge so große Kräfte nicht mehr herrschen, haben wir näherungsweise die 
berechneten Grenz beanspruchungen als fern wirkende Kräfte aufgefaßt und auf diese Weise 
die am Bauwerk derzeit wirkenden Kräfte ermittelt. Wir haben noch Gleichungen für die 
Berechnung des Zähigkeitskoeffizientes fester Stoffe angegeben “und die Ergebnisse mit den 
Versuchen verglichen. Eine zweite geologische Erscheinung, die näherungsweise und etwas 
kürzer als das klüftige Gebirge behandelt wurde, ist ein zwischen zwei Berge zusammengepreßtes 
geologisches Gebiet, das seine Verwirklichung am Semmeringpaß bei Wien findet. 

Der ausführliche Bericht des Vortrages erscheint demnächst in ‚‚Ingenieur-Archiv“. 


Spiegelung beim Torsionsproblem 
Von (©. Weber, Schlewecke/Nette 


Der Querschnitt eines tordierten Stabes sei durch zwei Kreislinien begrenzt, so daß ent- 
weder ein Querschnitt mit exzentrischem Loch oder, falls die Kreise sich berühren oder schneiden, 
ein sichelförmiger Querschnitt entsteht. Die Spannungsfunktion wird durch abwechselnde 
Spiegelung einer Singularität an den Kreislinien gefunden. Besonders anschaulich wird das 
Verfahren durch die konforme Abbildung des Querschnittes auf einen unendlichen Streifen. 
Die Untersuchung zeigt, wann Potenzreihenentwicklungen möglich sind. Auch noch nicht be- 
kannte geschlossene Lösungen werden mit Hilfe des Verfahrens angegeben. 


Bieseschwingunsen des Stabes mit nichtlinearem Elastizitätsgesetz 
Von F. Weidenhammer in Clausthal 


Die Bewegung eines Stabes wird durch. gekoppelte, nichtlineare partielle Differential- 
gleichungen für die Verschiebungsfunktionen v(z,t) und w(x,t) in Längs- und Querrichtung be- 
% 

2 
lehre hinausgehend nichtlinear berechnet. In folgerichtiger Näherung hat man daher auch die 
Nichtlinearitäten in der Spannungs-Dehnungs-Gleichung des Werkstoffes zu berücksichtigen, 
was durch Ergänzung des Hookeschen Gesetzes durch ein kubisches Glied im Anschluß an 
R. Grammel [1] näherungsweise geschehen soll: 


o()=Ee(l + A). 


schrieben, wenn man die Stabdehnung e= w,+ zW,„ über die klassische Schwingungs- 


_ Hierin bezeichnet E den Elastizitätsmodul und A eine zusätzliche, für technische Baustoffe meist 


negative Werkstoffkonstante. Zur Aufstellung der Bewegungsgleichungen hat man zunächst 
die kinetische Energie und die Formänderungsarbeit zu berechnen.: 


I 


[irda + VER). 


0 


I i ı l 
EF wz\2 EJ EJ,A 
T-# [(+wi) de; 4-3 [et 3)e+ [endeı I 
\ 0 


Hierin bezeichnet J, die Querschnittgröße [ zdF-—=x12J (x dimensionslose Zahl, © Flächen- 
trägheitsradius). In A; sind Summanden mit höheren als sechsten Potenzen von w fortgelassen 
worden, wobei zur Abschätzung der Größenordnungen die nachträglich kontrollierbare Annahme 
einer ortsunabhängigen, mit dem Hookeschen Gesetz berechneten Stablängskraft zu Grunde 


| gelegt wurde (vgl. E.Mettler[2]). Mit dieser Annahme ergibt sich nämlich die Stablängskraft 


: ey 2 
S-Er(w +) =" | Fan BRzsan () 
- 0 


des an beiden Enden unverschieblich (u = 0) gelagerten Stabes als quadratisch in w. Das 
- Hamiltonsche Prinzip liefert dann mit (1) zunächst für die Querbewegung eine nichtlineare 


E- Br, 066 
= partielle Differentialgleichung, in der die Längsverschiebung nur in der Verbindung u wi/2 


vorkommt, so daß man mit (2) eine nichtlinerare Integrodifferentialgleichung für w allein mit 
den Randbedingungen WW = 0 für = 0 und «= !herleiten Band: 


EI Wanna HUT = ar (dia Wa + BJ; An en u ER 


Für den- Grenzfall eines streng linearen Werkstoffgesetzes ist A=0 und (3) Bent Hand! in die 
bereits von S. Woinowsky-Krieger[3]und von D.Burgreen [4]auch experimentell 
untersuchte Integralgleichung der Stabschwingungen über. Sie läßt dann spezielle strenge 


- 
Lösungen w—= W,(t) sin T x zu, wobei W,„(t) jeweils einer nichtlinearen Schwingungsgleichung 


genügt und jede Eigenschwingungsform für sich allein bestehen kann. Im allgemeinen Faleder 
vollständigen Gleichung (3) ist dies nicht mehr der Fall, wie bereits R Grammel[1Jmit 

Hilfe eines Verfahrens der Störungsrechnung für nichtlineares Werkstoffgesetz mit inarem 
ae Bunpeanspsz gezeigt hat. Es liegt daher nahe, zur Lösung von (3) den Ansatz 


(,— & 1W,() sin T x zu verwenden und die dimensionslosen Zeitfunktionen W,auseinem nicht- EN - 


linearen System gewöhnlicher Differentialgleichungen zu berechnen. Um den Einfluß der Ab- 2 RER 
weichungen vom Hookeschen Gesetz auf das Schwingungsverhalten des Stabes in erster 
 Näherung zu erkennen, genügt ER bereits die Betrachtung einer einzigen, etwa der ersten Fer 
Zeitgleichung er - er 2 
wi + @ı W. m,(i A en 1* „| wi) =0 NT AR TI DAT (4. Fer 
‚Hierin bedeuten w, die erste Biegeeigenfrequenz, A den Schlankheitsgrad, a Querschnitts- hie 
® ‚größe zur Kennzeichnung von J,und A die zusätzliche Werkstoffkonstante. FüreinHoo ke- a er 
sches Werkstoffgesetz ist A = Ö, ‚so daß dann der Vorfaktor des kubischen Rückstellgiees 
Bin (4) stets positiv ist. Die freien Biegeschwingungen werden in diesem Falle also näherungs- Ser 
_ weise durch die Gleichung des Einmassenschwingers mit überlinearer Federkennlinie beschrieben. _ 
- Für technische Werkstoffe ist A jedoch von Null verschieden und negativ, so daß dieser Vor- 
faktor verringert wird, jedoch für normale Werkstoffe und Stababmessungen noch positiv bleiben 
Mi. dürfte. Der Charakter der Schwingungen als überlinear ändert sich also nicht, wenn auch de 
 Frequenz-Amplituden-Gleichung etwas andere Zahlenkoeffizienten erhält. — Nach Lösung von SR 
a) läßt sich die Annahme (2) auf gleiche Weise wie bei E. Mettler [2] rechtfertigen. Ins-| » 
besondere erkennt man, daß der ortsabhängige Anteil der Stabkraft (2) den Faktor u2jAR erhält ve 
nd mithin bereits für Stäbe mittlerer Schlankheit YernagpPasighar klein; ist. 
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Sagen äherte een aai nr Schalentheorie B Bra kic 
Fr mit Anwendung auf die Zylinderschalen er 
1 Von W. Zerna in Hannover u 


Die beherrschenden Gleichungen der Elastostatik dünner. Schalen RER] bekar 
hgewichtsbedingungen am Schalenelement und das Elastizitätsgesetz, das di 
wischen Schnittgrößen und Verschiebungen vermittelt. Während sich nu 
dingungen in mathematisch exakter Form ableiten lassen, müssen zur 
izit ätsgesetzes anschaulich physikalische Vorstellungen zu Hilfe 
h gewisse Annahmen eingeführt werden. Diese Annahmen bedin 
als erste 'Näherung anzusprechen ist. Es erhebt sich damit die Frage 
R treng abgeleiteten Gleichgewichtsbedingungen vielleicht Glieder e 
erste Näherung hinausgehen und damit logischerweise zu vernachläss 
benfalls angenäherte Gleichgewichtsbedingungen ergeben. 
Diese Frage soll nachfolgend behandelt und gezeigt werden, 
herung an Stelle der exakten Gleichgewichtsbedingungen angenä 
Zunächst wird die Betrachtung in ganz allgemeiner Form 
deutung is Ergebnisses am Beispiel a RER aufgezeigt In 
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Bezeichnen n°® den Längskrafttensor, q* den Querkraftvektor, m? den Momententensor, 
: pP, p die Komponenten der Belastung und b; die gemischte Form der Koeffizienten der zwei- 
ten Fundamentalform der Mittelfläche !), so lauten die exakten Gleichgewichtsbedingungen 2) 


me, = g°+ pP ER EN SE RE ET (1), 
a le 7 
m, — 1? —( Et a Far A en hs (3). 


Die Dehnungen bzw. Gleitungen der Mittelfläche werden beschrieben durch 
1 
av = 5 (Out Yon — Bar %) I N ee te (4), 
worin v,, w den Komponenten des Verschiebungsvektors entsprechen. 


Für die Verkrümmungen bzw. Drillungen der Mittelfläche ergeben sich in erster Näherung 
die Ausdrücke 


On Va a ee Sl: 
Nach dem Prinzip der virtuellen Arbeit gilt 
BUN AO A EL er Or 


Darin bedeutet öU die bei einer virtuellen Verrückung entstehende elastische Energie, ö*A, die 
virtuelle Arbeit der längs der Randkurve c angreifenden Randkräfte bzw. Momente und ö*A, die 
virtuelle Arbeit der Belastung. 


Die elastische Energie läßt sich durch | 
Ne ST er 


ausdrücken, wenn d,S das Flächenelement der Mittelfläche bezeichnet. Wird darin (3), (4) und (5) 
eingesetzt, so ergibt sich nach einigen Zwischenrechnungen 


Str. +P) tw b,+ ga tp)öu} dE=0 ....... 8) 

Da öv,, öw willkürliche und unabhängige virtuelle Verrückungen sind, so folgt aus (8) 
3A=E a en 
a 0a en ee end ach a EO: 


Diese Gleichungen treten an Stelle der Gleichgewichtsbedingungen (1) und (2). Ein Ver- 
gleich mit ihnen zeigt, daß (2) mit (10) übereinstimmt, während in (9) das Glied mit g* fehlt. So- 
mit stellt also (9) die angenäherte Form der Gleichgewichtsbedingungen (1) dar. 

Die Bedeutung dieses Ergebnisses soll jetzt am Beispiel der Zylinderschale gezeigt werden. 

Mit den üblichen Bezeichnungen lauten die (1) entsprechenden exakten Gleichgewichtsbe- 
dingungen der Zylinderschale ?) 


NEN. +tZ=U-N + N: -utrY=09-::% -: Srr ll) 
Punkte bedeuten Ableitungen nach 9, Striche nach =, wenn x und p Zylinderkoordinaten und a 


den Krümmungsradius angeben. 
Die (9) entsprechenden angenäherten Gleichungen sind 


ee ae ar er (12) 


die sich in einfacher Weise wie in der Scheibentheorie durch eine Spannungsfunktion identisch 
befriedigen lassen, was für den Aufbau der Theorie der Zylinderschale von Vorteil ist. 


3) i i hi he Buchstaben durchlaufen die Zahlen 1 und 2, über doppelt auftretende gleiche 
| Da akeätrend der Summationsregel zu summieren; kovariante Differentiation in Bezug auf die 


i äche ist durch einen senkrechten Strich angezeigt. 5 
RA E.G re und W.Zerna, Quart. Journ. Mech. and Applied Math. 3 (1950), 8. 9. 


®) Vgl.z.B,W, Flügge, Statik und Dynamik der Schalen, Berlin 1934, S. 110. 
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Ausführliche Referate folgender Vorträge: : 
E. Mettler: Bericht über die Theorie der nichtlinearen Schwingungen elastischer 
Körper. 
H. Richter:'Elasto-plastische Reflexion eines Stabes. 
H. Schaeffer: Die Spannungsfunktionen des dreidimensionalen Kontinuums 
und des elastischen Körpers. 
F. Schultz-Grunow: Greensche Funktionen für elastische Platten. 
G. Sonntag: Die Rohrschelle, ihre Beanspruchung und Wirkung auf das Rohr 
unter Berücksichtigung eines anfänglichen Spieles 


werden in den nächsten Heften dieser Zeitschrift erscheinen. 


Referate der Vorträge: 
H. Föp pl: Betrachtungen über den Einfluß kleiner elastischer und plastischer Ver- 
formungen auf die Versprödung von Stahl. 
L.N. Persen: Über Einflußfelder der Kreisplatte und der unendlich ausgedehnten 
Kragplatte 


sind nicht eingegangen. 


G. Strömungslehre 


Schmiermitteldruck und Randverformung des elastischen Rollenlagers 
Von Johannes Dörr in Darmstadt 


Um zu Aussagen über die Lebensdauer von Rollenlagern kommen zu können, muß man die 
im Betrieb auftretenden Schub- und Normalspannungen in den Rollen, der Lagerschale und der 
Welle kennen. Diese Spannungen lassen sich berechnen, wenn man die an den Lagerelementen 
angreifenden Randkräfte kennt. Diese Randkräfte sind aber identisch mit den im Schmiermittel- 
spalt zwischen den Rollen einerseits und der Schale bzw. Welle andererseits auftretenden Druck- 
und Schubkräften des Schmiermittels. 

Zur Berechnung der Schmiermittelkräfte genügt die Kenntnis der Spalthöhe als Funktion 
der Längenkoordinate. Bei einem als unverformbar angenommenen Rollenlager läßt sich die 
Spalthöhe aus dem Radjus der Rolle und dem Radius der Welle bzw. Schale unmittelbar an- 
geben, sofern man für die Mindestspalthöhe einen Wert vorgibt. Die Ermittlung der Schmier- 
mitteldrücke stellt dann lediglich eine Integrationsaufgabe dar. 

Beim verformbaren Rollenlager ist die Form des Schmiermittelspaltes nicht von vornherein 
bekannt. Die Spaltform wird nämlich wesentlich beeinflußt von den Normalverschiebungen, die 
die Oberflächen des Rollenrandes und des Wellen-" bzw. 
Schalenrandes unter dem Einfluß der an ihnen ‚wirkenden 
Kräfte erleiden. . : 

Zur Vereinfachung der folgenden Überlegungen führen 
wir einige unwesentliche Vereinfachungen ein. Wir nehmen 
nur die Rollen als elastisch an und ersetzen das wirkliche 
Rollenlager durch das in Bild 1 gezeigte Ersatzsystem. Die 
hierfür gefundenen Ergebnisse lassen sich mit einfachen 
Beziehungen auf das wirkliche System übertragen. [Ferner 
ZZ werde das System in Richtung der Rollenachse als beid- 

seitig unendlich lang angesehen, um mit ebenen |Verfor- 
mungszuständen und Druckfeldern rechnen zu können. 
Die Berechnung der Normalverformung eines mit perio- 
Re dischen Druckkräften der Periode x belasteten Rollenrandes 
läßt sich mit Hilfe der speziellen Green schen Funktion 
End 
G(g) = le cos pIn ale: 


ıE +2 5 7 king] 


streng durchführen. Sie stellt die Radialauslenkung durch zwei auf demselben Durchmesser 


Wirkliches System: 


Ersatzsystem: 


wZ 


liegende, zum Mittelpunkt weisende Einheitskräfte dar. Die Winkelkoordinate p wird von diesem 


Durchmesser aus gemessen. 


Da sich die Schmiermitteldrücke auf sehr kleine Winkelbereiche konzentrieren, dürfen wir 
tg (9/2) linearisieren. Außerdem dürfen wir in einem kleinen Winkelbereich den Rand der un- 
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verformten Rolle durch eine Parabel ersetzen. Mit den üblichen, von Reynolds genauer 
begründeten Vereinfachungen erhält man dann für den Schmiermitteldruck die Beziehung 


% 


b . 
OEL 1 BA) 5 (ia) 


wobei die Spaltweite durch 


© 
2 Sl —n (a 
B@=s + un Son ze +ı)a.... . m) 


beschrieben wird, und s ein frei vorgebbarer Parameter ist, der beim unverformbaren Lager die 
Mindestspaltweite darstellt. Durch n wird das Zähigkeitsmaß gekennzeichnet, Y, und V, sind 
die' substantiellen Geschwindig- 
keiten der Spaltränder, bezogen 


NN 
SINN) 
Schmiermitteldruck q(X N n& no) 


auf die Stelle engster Spaltweite. Se 
Alle andern Bezeichnungen ge-' \ N RR 
hen aus Bild 2 hervor. War SI x = 

.. der Randpunkt der Druckkurve Bon N \ A Sn ER 
%, vorgegeben wird, so bestimmt +5 FE ARER Ba BER RT URARTRENREEERENER Bir. Rose 
sich der andere Randpunkt ©, Y/ GH DSL DEE: WE, 0, Ti = 
u ERDE Rollenlager in der Umgebung der Stelle kleinster Spaltweite 


(2) = Il) =0 _ (le). Bild 2 


Bei starker Schmierung ist &,> x&,, wodurch der Grenzfall x, — oo mathematisch gerechtfertigt 
wird. Wir behandeln hier nur diesen technisch sehr wichtigen Grenzfall. 

Zur Lösung der durch (1) formulierten Aufgabe kommen nur .Iterations- oder Probierver- 
fahren in Frage. Das nächstliegende Iterationsverfahren besteht darin, für g(z) eine plausible 
Funktion anzunehmen, hiermit H(x) nach (1b) und die verbesserte Druckkurve nach (la) und 
(le) zu berechnen, um mit dieser verbesserten Druckkurve das Verfahren zu wiederholen. Die 
Schwierigkeit besteht dabei darin, daß sich (x, — &,) stark mit der relativen Größe der Spaltver- 
formung ändert. Das wird überwunden durch eine Koordinatentransformation 


Bit, güctm)=r, Hü+z)=NQ), 


durch die der Punkt x, in den Punkt £ = 0 und der Punkt x,, bei geeigneter Wahl von A in den 
Punkt & = 1 überführt wird. Die Beziehungen (1) lassen sich dann in folgender Form schreiben: 


h'(£) eh 
3 »0) la Fargo) ER SENSE EHE 2 (2a), 
£ 
h2(0)  [ [A°(0) A°(0) Pag 
Pe) wi | Kr ” 2) EEE (2b), 
Ko= FF NE (2c) 


Damit p(£) an der Stelle — 1 sein Maximum hat, muß gelten: 
1 
% 
h(0) =D) ER a= —1 on) K(&)d DR SE (3). 


Die Konstante b bestimmt sich aus der Forderung p() =0. Der Parameter » charakterisiert 
die Stärke der relativen Verformung und wird bei dem Iterationsverfahren vorgegeben. Durch 
den Quotienten K(£)/K (0,4) erreichen wir eine Normierung der Verformungsfunktion an der 
Stelle & = 0,4. 

‘Das Iterationsverfahren läuft wie folgt. Mit einer angenommenen Druckkurve nullter 
Näherung wird K(£) nach (2e) und a nach (3) berechnet. Durch Probieren und Interpolation 


> wird dann der Wert b festgelegt, indem zuerst h(£)/h(0) nach (2a) und dann die Druckkurve 


» 


ee re Fe 


EEE Pr a 
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1. Näherung nach (2b) mit der Eigenschaft p(oe) = O ermittelt wird. Mit den verbesserten Druck- 
kurven wird das Verfahren wiederholt, bis ausreichende Konvergenz eintritt. 

Bei diesem Verfahren ist die Berechnung der ‚‚Verformungsintegrale“ X (2) sehr lästig wegen 
des singulären Integranden. Diese Schwierigkeit wird überwunden durch Einführung von ele- 
— mentaren Druckfunktionen 
p,(£), die die Eigenschaft 

j 5 . haben, daß sich ihre zugeord- 
Druckkurven bei gleichem Rollenradius R neten Verformungsintegrale 
und gleicher Rollenlast pP K,(€) streng berechnen las- 
sen. Als 9,(£) wurde die 
Druckfunktion des unver- 
formten Rollenlagers gewählt. 
« als Verformungsparameter Für die übrigen 2,(£) mit 
v=1,2,3,... erwiesen sich 
Funktionen der Eorm 


ml — 2) 
Pre) O THE 


mitm =2 oder 3und n=3, 
4, 5 als geeignet. Die durch 
Iteration gefundenen Druck- 
kurven werden durch Reihen 
der p,(£) angenähert, wodurch 
sich die zugeordneten Verfor- 
mungsintegrale sofort als Rei- 
hen der K,(£) ergeben. 


Das Verfahren erwies sich bis zu relativ großen Verformungen als brauchbar. Bei sehr 
großen Verformungen ergaben sich Schwierigkeiten, die aber im Wesen der Aufgabe begründet 
sind, da sich bei starken Verformungen die Spaltweite als kleine Differenz großer Zahlen ergibt. 
Bild 3 zeigt die bisher ermittelten Druckkurven. Über weitere Einzelheiten, insbesondere über 
den Zusammenhang zwischen dem Verformungsparameter und den Betriebsdaten des Rollen- 
lagers wird bei anderer Gelegenheit berichtet werden. 


Bild 3 


Eine neue Reihenentwicklung für laminare Grenzschichten 
Von A. Görtler in Freiburg i. Br. 


Die bekannte Blasiussche Reihe für stationäre, ebene, laminare Grenzschichten hat den 
Vorteil der einfachen Anwendbarkeit, weil die erforderlichen Koeffizientenfunktionen sich rasch 
als Linearkombinationen von vorweg ein- für allemal zu vertafelnden universellen Funktionen 
berechnen lassen. Trotz der heute weit fortgeschrittenen Vertafelungsarbeiten kann man jedoch 
wegen der geringen Konvergenzgüte kaum hoffen, unter Verwendung der mit praktisch erreich- 


. barer Ordnung abgebrochenen Reihe eine Grenzschicht vom vorderen Staupunkt bis zu einer 


eventuell vorhandenen Ablösungsstelle mit ausreichender Approximation zu berechnen. Ent- 
sprechendes gilt für Abwandlungen oder mögliche Verallgemeinerungen der Reihe sowie für die 
zugeordnete Reihe für die Querkomponente der Geschwindigkeit bei Schräganblasung des 
Zylinders }). 


Diesem Versagen der Blasiusschen Reihe liegen vor allem zwei Tatsachen zugrunde: 
1. Die äußere Randbedingung, wonach die wandparallele Geschwindigkeitskomponente asympto- 
tisch in die vorgegebene äußere Geschwindigkeitsverteilung U(z) (x Wandbogenlänge) überzu- 
gehen hat, wird durch gliedweises Abgleichen der Blasiusschen Reihe mit der Reihenent- 


wicklung von U(x) erfüllt — genügend’genau also erst bei Berücksichtigung von Gliedern ent- 
sprechend hoher Ordnung —;-2. Die benutzten Koordinaten — neben x der senkrechte Wand- 


abstand y — haben, da Krümmungseinflüsse in der Theorie vernachlässigt sind, nur die Bedeu- 
tung cartesischer Koordinaten in der Ebene der Abwicklung der Kontur auf die x-Achse und sind 
nicht den speziellen Daten des jeweiligen Problems in besonderer Weise angepaßt. 


1) Einige bierzu gehörige universelle Funktionen für einen Sonderfall findet man in Diagrammen (kei 
Zahlentafeln) dargestellt bei W.Sears, Journ. Aeron, Sei. 15 (1948), 8.49 —52. Sorgfältig he Yale 
tafeln für den allgemeinen Fall werde ich in Kürze im Arch. Math. veröffentlichen, EN Eh 


s . 
‚r s j en 27 
. Z. angew. Math. Mech, ee. 
Bd. 32 Nr. 8/9 Aug./Sept. 1952 G. Strömungslehre 971 


Es wird-eine neue Reihenentwicklung vorgeschlagen unter Verwendung der unabhängigen 
. Koordinaten ?) 


8 [Uga, n=U) 7 I [va | 


(v kinematische Zähigkeit). Die Reihe stellt eine Entwicklung nach Potenzen von & dar mit von 
abhängigen Koeffizientenfunktionen, die wieder auf Linearkombinationen universeller Funk- 
‚tionen zurückgeführt werden können. Das erste Glied der neuen Reihe erfüllt die äußere Rand- 
bedingung exakt, die weiteren Glieder bringen nur Korrekturen im Inneren der Grenzschicht. 
Bereits das erste Glied der neuen Reihe gibt eine brauchbare Approximation der Grenz- 
schicht vom vorderen Staupunkt (bzw. Vorderkante bei Anströmung von Konturecken oder 
-spitzen) bis weit stromabwärts, wo das entsprechende erste Glied der Blasiusschen Reihe 
überhaupt keine Approximation mehr liefert. Dieses erste Glied ist im wesentlichen ein Har- 
tree-Profil, jedoch mit der obigen Variablen n, welche das örtliche U(x) und die Vorgeschichte 
berücksichtigt. Die Vertafelungsarbeiten für die Glieder höherer Ordnung sind in Angriff ge- 
nommen. 
Aus Raumgründen muß hier darauf verzichtet werden, auf verschiedene im Zusammenhang 
mit der neuen Reihe stehenden theoretischen und praktischen Fragen auch nur andeutend ein-' 
zugehen. Ein ausführlicher Bericht wird vorbereitet. Von zentraler Bedeutung ist eine gewisse 


Ü 


Funktion ß(£) (deren Entwickelbarkeit in I ß,&” vorausgesetzt wird, worin ß,der Hartree- 
N | 


3 v= d 
Parameter ist und durch die Verhältnisse an der Vorderkante bestimmt wird). Als einziger 
Koeffizient der Differentialgleichung in den neuen Variablen enthält diese Funktion die Daten 
des Problems, während die übrigen Koeffizienten und die Randbedingungen universell sind. 


Zur Energieverteilung in strömenden Medien 
2 Von G. Heinrich in Wien 

Die Gleichung für die stationäre Bewegung der homogenen und reibungsfreien Flüssigkeit 
läßt sich bekanntlich auf die Form bringen: | 

- ee ; KR j v2 
ie 0x (ex) =r[ 
4 z ae Lane 
mit P als Druckfunktion und U als Potential der Massenkräfte. 

" _ Wir setzen vx(7/xb)=bxmw=t(m ist Wirbelvektor). Da t die Translation des 
_  Wirbelvektors mißt, soller, nach einem Vorschlag von L.Flamm, als Translator bezeichnet 
werden. Der Translator spielt hier die Rolle eines Energiedichtegradienten und ist daher für 
die homogene und reibungsfreie Flüssigkeit wirbelfrei, während er im allgemeinen Wirbel 
es Nach dem Stokesschen Satz gilt also für homogene und reibungsfreie Flüssigkeiten 
ei stationärer Bewegung: alle a 
ne  /frxd d=dt-da=0. EEE 3 
$ Linien-Integral dt: dx des Translators läßt sich aber leicht mit der Zirkulation [= $ ddr 
in Zusammenhang bringen. Bildet man nämlich den substantiellen Differentialquotienten. 


4240) 


ne k £ at N 
ee DET Di dr Di m 
ausder Um formung 000.0 au R 
ER a ein &ch ae | 
EBEN TE 


y. Soc g5 (1948), 8. 545—576 und spätere Arbeiten in dergleichen Zeitschrift] 
rsuchungen Koordinaten, die er aus Potential und Stromfunktion 

ter Weise bildet. Unsere Koordinaten sind im wesentlichen 
| 1 hier exakt eingeführt, mit dem Vorteil der un- 
g der Meks ynschen Koordinaten bei jedem 
_W. Mangler (vgl. W. Tollmien, 
ew o Science 1939—1946, Band ‚‚Hydro- und 
erer Variablen & vorgeschlagen, benutzt aber zur Bil- 
tion der Grenzschicht, die sich erst mit der Grenz- 
mständlicher ist. Den Hinweis aufdenMang- 
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für stationäre Strömung: 


Dr _ Ko (®®) 2004 balt) = br: ar baum. 
Ir) arpai +p a) 
Das zweite Integral verschwindet längs eines geschlossenen Weges, daher verbleibt: 


DT 
$t-d= a 
Das Linienintegral .des Translators längs einer geschlossenen Kurve ist also gleich dem 
negativen substantiellen Differentialquotienten der Zirkulation längs derselben Kurve. Daraus 
kann man auch sofort auf den Thomsonschen Satz von der Konstanz der Zirkulation längs 
einer geschlossenen flüssigen Linie bei reibungsfreien homogenen Flüssigkeiten schließen. Bei 
inhomogenen Flüssigkeiten folgt aus dem Satz vonBjerknes: 
L 


Pxt=r(O) xVp. 


Der Translator t = x iw hat die Form des bekannten Poyntingschen Vektors der 

Elektrodynamik: 
6 
\ S=-,,(Ex 9: 

Dies legt die Frage nach einem Energietransportvektor der strömenden Flüssigkeit nahe. Der 
Vektor dx mw hat jedoch nicht die Bedeutung eines Transportvektors der Energie. Ein solcher 
ist, soweit mir bekannt, für den Fall einer Flüssigkeit, die den Stokes-Navierschen Glei- 
chungen genügt, noch nicht eingeführt worden. Der Energietransportvektor © soll so definiert 
werden, daß dE =(©& do) dtdie in der Zeit dt durch das Flächenelement do transportierte 
mechanische Energie ist. Diese setzt sich zusammen aus der durch do hindurchtretenden Be- 
wegungsenergie, der potentiellen Energie und der mechanischen Arbeit, die von der durchtreten- 
den Flüssigkeit während di geleistet wird. 


Nach dem Stokes-Navierschen Ansatz ist der Spannungsaffinor p gegeben durch: 

p = —I7+nW,0p +v,V/)—AS(V »d). 
® ist hierin der Einheitsaffinor, V, vd und v, V bedeuten dyadische Produkte. Nach Stokes ist 
= an doch wollen wir diese etwas problematische Annahme hier nicht einführen. Das 


Skalarprodukt p - do ist somit der Spannungsvektor der im Element do auftretenden Schnitt- 
spannung. Die Arbeit, die die durch do hindurchtretende Flüssigkeit während dt leistet, ist gegeben 
durch d44A= —(p -do)-d.dt, oder, wegen der Symmetrie von p: dA=—(p -v):dodt. Die 
gesamte durch do transportierte Energie ist daher: | 


dE (+ U )ov- ddp -d)-dodi 


und daher der Energietransportvektor:. 


S-( +U)en—r:v. 


Durch Einsetzen des Wertes von p ergibt sich 


&= (2 ee O)en nd: (W040, + Anl): 


/ 


Für die reibungsfreie Flüssigkeit (n = 0 und A = 0) hat somit © die Richtung von d, für 


die zähe Flüssigkeit schließt © i. a. einen Winkel mit d ein. 


Bildet man div © so erhält man, bei Benützungder Stokes-Navierschen Gleichung. 


in der Form 


90 pi2 1 | | 
Dr 53 x WR) = PP PU—PER + Pr (7-b) 


den Ausdruck: 


\ () we a I 
ivo-rp.s- (ehr eo) + pr 0) nW,0) (9,0 +0, P)+AW-d)% 
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div. © stellt aber die Quellung der Energie pro Volums- und Zeiteinheit dar, und diese ist gleich 
der Abnahme der Bewegungs- und der potentiellen Energie, ferner einem von der Volumsdilata- 
tion abhängigen Glied, das von einer thermodynamischen Umsetzung herrührt und der Energie- 
dissipation-durch innere Reibung. 
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Indirekte Bestimmung von Lüfterkennlinien 
Von W. Richter in Frankenthal 


Wenn in einem einfachen geschlossenen Leitungssystem mit einer Queiverbindung zwei 
Lüfter arbeiten, so können drei wesentlich voneinander verschiedene Aufgaben vorliegen: voll- 
ständige Hintereinanderschaltung, Parallelschaltung und teilweise Hintereinanderschaltung. 
Wenn im passiven Teil des Leitungssystems verschiedene Drosselzustände eingestellt werden 
können und einer der beiden Lüfter zumindest in einem begrenzten Bereich regelbar ist, so kann 
man aus einer Messung dieses Lüfters auf die Kennlinie des zweiten Lüfters schließen. Die Auf- 
gabe führt jeweils auf eine gewöhnliche lineare Differentialgleichung erster Ordnung. Für deren 
graphische Lösung kann eine Leitkurve (E.Czuber, Z. Math. Phys. 44 (1899), S. 41—49) ver- 
wendet werden, deren Parameterdarstellung angegeben wurde. Aus der Möglichkeit, im normalen 
Lüfterarbeitsbereich die Kennlinie eines Lüfters gut durch eine Parabel zu ersetzen, wurden rech- 
nerische Näherungslösungen hergeleitet. Diese graphischen und rechnerischen Verfahren wurden 
an praktischen Beispielen erläutert und deren Ergebnisse miteinander verglichen. Eine prak- 
tische Verwendung empfiehlt sich wohl nur dann, wenn man die Kennlinie des zweiten Lüfters 
direkt nicht bestimmen kann. Insbesondere erhält man auf diese Weise zwei verschiedene Mög- 
lichkeiten, den natürlichen Auftrieb einer Grube durch Messungen am Hauptgrubenlüfter zu be- 
stimmen. Ihre Anwendung setzt voraus, daß die Messungen am Hauptgrubenlüfter mit hin- 
reichender Genauigkeit möglich sind. 


Ein neues Verfahren zur Berechnung turbulenter Grenzschichten 
Von J. Rotia in Göttingen 


Das neue Berechnungsverfahren, dessen theoretische Grundlagen in vorausgegangenen 
Veröffentlichungen !) niedergelegt wurden, gilt für turbulente Grenzschichten inkompressibler 
Strömungen mit beliebigem Druckverlauf. Für die Bestimmung der Impulsverlustdicke wird die 
bekannte Impulsgleichung benutzt. Für die dabeiin dem Verhältnis H,, der Verdrängungsdicke 6, 
zur Impulsverlustdicke ö, ausgedrückte Profilform und den örtlichen Reibungsbeiwert c; sind 
noch weitere Beziehungen erforderlich. Eine dieser Beziehungen liefert die Energiegleichung, die 
das Gleichgewicht zwischen der der Grundströmung entzogenen Energie, der dissipierten Energie 
und der mit der Grundströmungsgeschwindigkeit fortgeführten Turbulenzenergie beschreibt. Die 


‘ übrigen Zusammenhänge werden im wesentlichen durch drei grundsätzliche Annahmen vermittelt, 
_ die für das Verfahren kennzeichnend sind: 


1. Die Grenzschicht kann aus einem wandnahen und einem äußeren Teil zusammengesetzt 
werden. 
2. Die äußeren Teile der Geschwindigkeitsprofile können als einparametrige Familie auf- 
gefaßt werden. RE 
3. Der Zusammenhang zwischen den Profilen der Grundströmungsgeschwindigkeit und den 
Verteilungen der statistischen Schwankungen kann als eindeutig angesehen werden. 
Auf Grund der ersten und zweiten Annahme läßt sich der Reibungsbeiwert c; als Funktion 
des gebräuchlichen Profilparameters H,,; der Reymoldszahl und der Wandrauhigkeit dar- 
stellen. Für den Reibungsbeiwert bei glatten und sandrauhen Oberflächen wurden Kurven- 


tafeln angefertigt. 


- 


1) J. Rot ea , Mitteilungen aus dem Max-Planck-Institut für, Strömungsforschung, Nr. 1, “Göttingen 


1950. J. Rotta , Schubspannungsverteilung und Energiedissipation bei turbulenten Grenzschichten, Ing.- 


Arch. 20 (1952), 8.195. _ | 
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Nach Einführung einer Gesamtenergie-Verlustdicke 


U Ei: nn 
= I — |1—1|;r) — dy, 
| 0, | 0, nase 
bei der die Umwandlung von Grundströmungsenergie in Turbulenzenergie nicht als Energiever- 
lust betrachtet wird, läßt sich die Energiegleichung unter Benutzung der Impulsgleichung in eine 
Differentialgleichung für die Änderung des Formverhältnisses 43, = 63/6, umwandeln. Dabei ist 
U die Grundströmungsgeschwindigkeit, U, deren Wert außerhalb der Grenzschicht und u, dv, W 
die Komponenten der Geschwindigkeitsschwankungen. Man kann diese Differentialgleichung in 
einer Fluchtlinientafel darstellen (Bild), wobei die Energiedissipation durch die dimensionslose 


Dargestellte Beziehung 


p.: statischer Druck 
? 


9,7% 


0 -4 
Fluchtlinientafel | 
zur Bestimmung von 6> dH%ldx 


Größe c) berücksichtigt wird. Reynolds zahl und Wandrauhigkeit werden dabei durch den 
örtlichen Reibungsbeiwert ausgedrückt, so daß diese Tafel sowohl für glatte als auch für rauhe 
Oberflächen gilt. Auch für Rotationskörper und Keilströmungen hat diese Tafel unverändert 
Gültigkeit. Das Formverhältnis 35 ist eine Funktion von A,, und c5. Dieser Zusammenhang 
kann näherungsweise ebenfalls als Fluchtlinientafel dargestellt werden, die auch wieder in gleicher 
Weise für glatte und rauhe Oberflächen gilt. Die Zusammenhänge ergeben mit der Impuls- 
gleichung zusammen ein System von 2 Differentialgleichungen, das sich mit gebräuchlichen Inte- 
grationsmethoden lösen läßt. : ; 

Die Festsetzung der Anfangswerte geht davon aus, daß Impulsstrom und Energiestrom 
beim Übergang vom laminaren zum turbulenten Zustand keinc Änderungen erfahren, so daß Ö,, 

3 und somit auch A435, unverändert bleiben. Liegen diese Daten der laminaren Anlaufstrecke bis 

zum Umschlag vor, so sind hinreichende Angaben über die Anfangswerte der turbulenten Rech- 
nung bekannt. | 

Das Verfahren ist durch folgende Besonderheiten ausgezeichnet: 


1. Es ist für glatte und rauhe Oberflächen anwendbar. 

2. Der Einfluß des Grenzschichtdickenverhältnisses H,, auf den Reibungsbeiwert wird 
gemäß experimentellen Beobachtungen berücksichtigt. 

3. Die allen bekannten Verfahren zugrundeliegende Annahme der Einparametrigkeit der 
Geschwindigkeitsprofile ist durch den Aufbau der Grenzschicht aus 2 Teilen verfeinert, so daß das 
Grenzschichtprofil nicht durch das Verhältnis 7, allein, sondern durch H,,und c; beschrieben wird. 

4. Die Frage der Anfangswerte findet durch den stetigen Verlauf von Impuls- und Energie- 
strom eine befriedigende Beantwortung. ; 

5. Die Grundlagen des Verfahrens sind theoretisch begründet, so daß man die Grenzen der ° 
Leistungsfähigkeit überblicken kann. 


l 
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Bemerkungen zu dem Begriff „Kritischer Turbulenzgrad“ a 
Von H. St. Stefaniak in München 4 


"Bekanntlich führten die Versuche von G.B. Schubauerund H.K.Skramstad!), . 
welche zur Nachprüfung der Tollmien-Schlichtingschen Theorie der Stabilität dr 
Grenzschicht an der längsangeströmten ebenen Platte durchgeführt wurden, zu-einer vorzüglichen B. 
Übereinstimmung zwischen Theorie und Experiment. Bei dieser Gelegenheit wurde eine Meß- :B 


reihe aufgenommen, bei welcher der Einfluß des mittels Hitzdraht festgestellten Turbulenz- 


= grades Ty = Un/U. auf die Lage Re, = U, - »/v derjenigen Stelle untersucht wurde, von 

% welcher an der Umordnungsprozeß von dem laminaren Strömungsvorgang in den turbulenten 

je einsetzt. (U. = Anströmungsgeschwindigkeit, %y7 = Wurzel aus dem mittleren Quadrat der 

E mittels Hitzdraht festgestellten turbulenten Zusatzgeschwindigkeiten, © — Lauflänge der Grenz- 

E schicht, » = kinematische Zähigkeit.) Das Bild gibt den Befund wieder. Man erkennt, daß sich 

g die kritische Stelle, falls man von großen und kleinen Ty-Werten fortschreitet, immer r weiter 

Be stromabwärts verschiebt, und zwar ist die 

4 Abhängigkeit Re,(T7;) sehr gut linear bis Al SSR 

e _ herunter zu dem ‚kritischen Turbulenzgrad‘?) T @ E: | ] 

rn Ta = 0,1%. Normalerweise würde man ver- Wrlalenen Yra | = 

=: muten, daß sich die Gerade, wie gestrichelt nn 

4 Fr eingetragen, bis zu dem Punkt Re,; bei 

N Ta =Ofortsetzt, bei welchem die Grenzschicht ü 

* beginnen würde, spontan umzuschlagen. So- 

0. weit konnte sich aber bei den Versuchen die 

; a Grenzschicht nicht entwickeln; denn unterhalb al 
E: des kritischen Turbulenzpegels knickte die wer 

Gerade plötzlich um und verlief ohne Über-- „ Br 
Pr influß des Turbulenzgrades auf die Lage der Übergangs- 

u, + . gang bis Ta — 01 in Höhe von Re, = 2 2 8: = zone in der Grenzschicht an der längsangeströmten ebenen 

Br parallel zur Äbszisse, Platte ohne Druckgradient [1] 
Be. . Wenn man sich die Frage Worlket; wie man diese Erscheinung deuten kann, so erkennt man — 


Pr 


n EN als Voraussetzung für einen bis Ty-= 0 ununterbrochenen geradlinigen Verlauf der Grenzkurve. 
£ u die Tatsache, daß außer der einen, von dem Hitzdraht erfaßten Störung keine weiteren mehr vor- 
_ handen sind, welche einen vorzeitigen Umschlag herbeiführen können. 


ES IE n Beidem Erklärungsversuch, welcher zur Diskussion gestellt werden soll, wird angenommen: 


e 1. Es existiert tatsächlich noch eine zweite Quelle von Störungen, welche von der Hitz- 
_ drahtanordnung nicht registriert wird. 

2. Diese Störungen lösen auch einen Umschlag aus; jedoch ie der innere Mechanismus, 
welcher zur Instabilität Anlaß gibt, anderer Art als der bei der Windkanalturbulenz vorhandene. 
2 3 Die Ursache der andersgearteten Störung ist in der un Unruhe zu suchen. 

. Zur Begründung wurde angeführt: | ; 
wi ‚Die Stoßzahlen der Moleküle in der Luft vom Norsalsoständ liegen in der Größenord- 
10°—1010 1/s. Es wurde daher die unregelmäßige thermische Bewegung von der Hitz- 
htanordnung, deren Wiedergabebereich sich auf a zwischen 1 bis einige 10°Hz er- 
| mit erfaßt. BiN- 
nen Schwankungen haben endlic B e Aaipltude, Es beginnen daher die 
san ansteigenden Schwingungen bereits mitendlicher Ansse 

\ renzschi ht ‚nach Durchlaufen einerendlichen Wegstrecke a 


‚hat immer Molekülhaufen zum Träger i im Gegensatz zur thermi- 
von Back an ist, a man eh a auch die 


1, BER Anber der durch die Sana urhnler 
ite, YOmdjeReE, a Dbaneis —im Ei a nunkten h 


iges Ende setzt. Pi 
tätskurven (Bild) läßt sich eine Ver- 
== Turbulenz Re, =üg :L/vy und. 
fo miteinander ae E = 
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die aus Korrelationsmessungen sich ergebende Länge (= scale of turbulence), I = freie Weglänge, 
% = mittlere Geschwindigkeit.) Aus ®) lassen sich für Re, für die durch 0) gekennzeichneten Meß- 
punkte mit einiger Sicherheit die Grenzen 0,7< Re, < 6,3 und als wahrscheinlichsten Wert 2,1 
erschließen. Mit » = 16/3 - 78/(8 x) errechnet sich der von den Zustandsgrößen des Gases unab- 
hängige Wert Re; = 27/8 = 3,26, welcher sich als etwas größer als die Reynol ds sche 
Zahl der bei den Versuchen vorhandenen unterkritischen Turbulenz erweist. Es liegt daher nahe 
anzunehmen, daß der ‚‚kritische Turbulenzpegel“ durch die Gleichheit von Re, mit Re, definiert 
ist. Aus den angestellten Überlegungen folgt aber, daß zur vollständigen Klärung dieses Phäno- 
mens Versuche bei gleich geringem 7’ aber verschiedenen Z und dann vor allem mit Flüssigkeiten 
durchgeführt werden müßten; denn die Ableitung, welche zu dem Zahlenwert Re, = 3,26 führte, 
stützte sich nur auf die kinetische Theorie der Gase. 


Berechnung der Schaufelform von Kaplan-Laufrädern 
bei vorgegebener Geschwindigkeitsverteilung 
auf der Schaufeloberfläche 


Von M. Strscheletzky in Friedrichshafen 


[e” Die Umströmung der Laufschaufeln von Kaplan-Turbinen wird durch gewisse gegen- 
seitige Beeinflussung des Laufrades, der Leitvorrichtung und des Saugkrümmers gesteuert. Ist 
diese Gegenwirkung rechnerisch erfaßt, so entsteht die Aufgabe, die Schaufelform des Laufrades 
bei vorgegebener Geschwindigkeitsverteilung auf der Schaufeloberfäche zu berechnen. 

Die dreidimensionale, wirbelbehaftete Strömung im Arbeitsraum der Turbine läßt sich 
ö rechnerisch erfassen, wenn man die Wände des Turbinenhohlraums und die Schaufeloberfläche 
: als Stromflächen reibungsloser Strömung betrachtet, welche unter dem Einfluß von Wirbeln, 
Br Quellen und Senken erzeugt wird. | 
Br Es seien: w,(r, 4) die lokale Geschwindigkeit der rotationssymmetrischen, wirbelbehafteten 
Be ' Strömung, v(r, Y, ©) die von den gebundenen Wirbeln der Laufschaufeln induzierte Geschwindig- 
keit, und Av(r, y, ©) die von den Quellen und Senken der Leitschaufeln erzeugte Geschwindigkeit. 

Die lokale Geschwindigkeit &, + ist zu der Skelettfläche der Schaufel tangential gerichtet. 
Für genügend dünne und schwach gewölbte Schaufel kann man den Gegeneinfluß der Wölbung 
und der Schaufeldicke vernachlässigen. Die Umrechnung der an der Skelettfläche auftretenden 
Geschwindigkeit auf die Oberfläche der Schaufel wird dann nach der angenäherten Methode von 
Riegels durchgeführt. Die lokale Geschwindigkeit auf der Schaufeloberfläche ist 


u (2) = [wo(a) + a) —— +5 
+1/2 
- Ef MED a 


FERLA, 
Vi) 
ya dyy.\2 G— % dz 
| | ex | 
m y: + ( e 

—/2 
In den Stellen + y, der Druck- bzw. Saugseite der Schaufel, welche dem Aufpunkt %, der 


Skelettlinie im betreffenden Schnitt der Schaufel entsprechen, sind die lokalen Geschwindigkeiten 
w, und w, vorgegeben. Damit ergibt sich aus (1) die lokale Intensität der gebundenen Wirbel 


I u 
BEERAUTN Gi 


\ 


y(%) = (w, — wu)‘ Yırl DA (2). 


An Hand von y ermittelt man die Geschwindigkeit %. Da w ur 1 AL Der 
ist, ergibt sich aus (1) und (2) die Integralgleichung (0) 2 (ws +0) + 9 (w, — w,) 
| +1/2 


1, wtw 1 day wat) 1 m dy d 
ee re er ea 


—j2 


®)H.L.Drydenand G.B.Schub Th i | 
Tunnel Turbulence, Journ. Aeron., ae 14 1947), 22 net a 


vw. 
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wobei der Index A der Austrittskante entspricht und die Integration längs der betreffenden Strom- 
linie @,+- ® auf der unendlich dünnen Skelettfläche der Schaufel durchzuführen ist. 


Die Gl. (2) kann man nach verschiedenen Methoden lösen. Im Vortrag wurde ein Iterations- 
verfahren angeführt, welches die Lösung mittels der von Truckenbrodt errechneten Bei- 
werte ermöglicht. 


Die im Vortrag angeführten Beispiele zeigen, daß der dreidimensionale Charakter der Strö- 
mung im Arbeitsraum der Turbine für die Schaufelform von entscheidender Bedeutung ist. 


Die Berechnung der Auftriebverteilung an Pfeilflügeln 
Von E. Truckenbrodt in Braunschweig 


Eine anschauliche Auffassung des Tragflügels erhält man, wenn man ihn aus sogenannten 
Elementarflügeln, d.h. Flügeln vom Spannweitenelement dy’, zusammensetzt. Gezeigt ist eine 
schematische Darstellung für einen schiebenden Pfeilflügel. Jedem Elementarflügel ist aufgrund 
der Geometrie des Flügelgrundrisses ein ganz bestimmter Platz im Gesamtverband des Flügels 
zugeordnet. Es lassen sich mit dem beschriebenen Verfahren Flügel beliebigen Grundrisses, 
beliebiger Profilwölbung und beliebiger Verwindung längs Spannweite behandeln. | 


Die Verteilung der Wirbeldichte % längs Flügeltiefe wird aus zwei Anteilen zusammen- 
gesetzt. Die erste Verteilung, die Grundverteilung A, entspricht der ersten Birnbaumschen 
Normalverteilung; sie ist proportional dem Auftriebsbeiwert c, des Flügelschnittes y’. Die zweite 
Verteilung, die Nullverteilung A,, ergibt sich durch Überlagerung der zweiten Birnbaumschen 
Normalverteilung mit einer flächengleichen ersten Normalverteilung; sie ist proportional dem 
Momentenbeiwert c,, des Flügelschnittes y’. Es gilt also für die Verteilung der Wirbeldichte 
eines Elementarflügels 


2 
ka, y) = aly') hl, y) +m(y') hola’, Y)- 


Unter Anwendung des Biot-Savartschen Gesetzes lassen sich für diese Wirbel- 
verteilung die induzierten Normalgeschwindigkeiten am Ort der Flügelfläche berechnen. Die 
Bestimmung wird entscheidend vereinfacht durch Verwendung von sogenannten Einflußtafeln 


EEE 


4 (® 9, y’)und x, (@, Y, y’)(entspre- 
chend den vorgegegebenen Vertei- 
lungen A und A,). Diese Einfluß- | 
tafeln sind universell und gelten 
für alle vorkommenden Flügel- 


ie 


grundrisse, Elementarflügel| ———]|_ Gebundenes 
Die Abwindbedingung, die iv) | le 

fordert, daß die aus der Anströmung =Kax 

herrührende Normalgeschwindig- 

keit entgegengesetzt gleich der vom 

Wirbelsystem des Flügels induzier- 

ten Normalgeschwindigkeitist, wird [Ei 

wegen der zwei Unbekannten c, und ı) ing are 


C„ an zwei Stellen eines jeden Flü- 
gelschnittes yerfüllt. Als besonders 
charakteristisch haben sich hierfür 
die Hinterkante und der Einviertel- 
punkt ergeben. | > 

Damit erhält man zur Bestimmung der Auftriebsverteilung und der Momentenverteilung 
zwei Integralgleichungen, deren Aufbau der Prandtlschen Abwindformel ähnlich ist. Hat 


b=2s 


“ man die Werte der Einflußfunktionen x(z, y, y') und x%.(&, , y') mittels der Einflußtafeln 


bestimmt, so löst man die beiden Integralgleichungen unter Anwendung der bekannten Ver- 
fahren zur Berechnung des Prandtlschen Abwindintegrals schrittweise. Am geeignesten hat 


sich dabei die Multhoppsche Quadraturformel erwiesen. 


19 
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Zur Theorie der Strömung in Kaplan- und Propellerturbinen. 
Von E. Weinel in Jena 


Bei der üblichen Berechnung des Strömungsfeldes in Kaplan- und Propellerturbinen sieht 
man die Strömung in axialen Zylinderschnitten als ebene Gitterströmung an und umgeht auf 
diese Weise die Schwierigkeiten des dreidimensionalen Problems. Da von Seiten des Strömungs- 
maschinenbaus immer wieder betont wird, daß eine theoretische Erfassung der räumlichen 
Strömung notwendig und wünschenswert ist, möchte ich über einige einfache potentialtheo- 
retische Rechnungen berichten, die sich für ein schematisiertes Turbinenmodell (Bild 1 und 2) 
durchführen lassen. 


1. Wir betrachten ein konzentrisches Rohr, das den Arbeitsraum der Turbine idealisieren 
möge, und einen Wirbelstern, der in erster Näherung die Schauflung ersetzen soll (Bild 2). Der 


Bild 1 Bild 2 


Geschwindigkeitsvektor B%=d*-+d setzt sich additiv zusammen aus der durch den Wirbel- 
stern induzierten Geschwindigkeit d* und einer Zusatzgeschwindigkeit b, die bewirkt, daß die 
Zylindermäntel r, und r, Stromflächen werden. 


2. Das Strömungsfeld eines s-strahligen Wirbelsterns ergibt sich unmittelbar aus dem 
Biot-Savartschen Gesetz. Für die Komponenten des Geschwindigkeitsvektors b* erhält 
man in Zylinderkoordinaten: 


3, sin s 
u LE erh ae Dark 
4 Pe Cof sA —cos sp ’ 
Z I Sin 


nr Anr BojsA—cos sp ’ 
a ERERE sin sp 


ee BE al IN at 
* Anr CofA LofsA—cossp ' 
Hierin ist A = X ©in (z/r), 7, bedeutet die Schaufelzirkulation und I'd i i i 
aus Eintritts- und Austrittsdrall. Br ee. 


3. Das Strömungspotential ® der Zusatahrarn irdi i : 
der Forma ung wird in Zylinderkoordinaten dargestellt 


9 £, oo 
&(r,9,2) = Ei sin | {An (a) Ju(iar) + B. (0) Y, War)}da 
0 


N= 


% 
. 
$ 


i 
t 
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mit den Se 5 > 


BE” 2 4% =0 - ürernzon ud r=n. 


Zur Kroithlüng der en A„(a) und B,(&) entwickelt man v* in die Fo uriersche Reihe 


2 sI, 


nel 32 De a sin np (n =0 mod 5) 


und erhält dann vermöge der Randbedingungen die beiden Enkel clngen 


ER 6) 


Aniar, 00 } 


An (&) In(war,) + B,(a) Yn ee — e®P cos (ro Sin B)dP, 


&.>} 


u 
Sa ER? cos (r, a Sin P)dß. 


Anıar, 0% J 


Anl) da (tan) + Bao) Yaliarı) = 


Die rechts stehenden Integrale vom Typus 


E Bes y Je ”P cos (x Sin B)dB (n = 0 mod s) 


< erden; je nachdem, ob n gerade oder ungerade ist, entweder durch Schläflische Polynome in Ver- 
_ bindung mit Zylinderfunktionen oder durch Anger- -Webersche Funktionen mit imaginärem 
Argument dargestellt. eh 

R: . Hiermit ist das Geschwindigkeitsfeld vollständig bestimmt. Grobe Überschlagsrechnungen 
- lassen vermuten, daß der Einfluß der Zusatzströmung i i. a. nicht sehr erheblich sein wird, so daß 
eine auf dem Prinzip der tragenden Linie bzw. tragenden Fläche aufgebaute Theorie auch bei 
SE Vernachlässigung der Wandeinflüsse bereits einen ersten Einblick in das dreidimensionale 
=; Problem vermittelt. z ’ re: 


_ Ausführliche Berichte sender Vorträge: Nur 

a “ : KR,Ja ekel: Zur Theorie der Profile geringer Dicke und Wölbung. 
5 e ie en ‚< R. Re i chhar dt: Über die Energieversorgung der Wandturbulenz. 

i SEES. auer: Fortschritte der Theorie der. Überschallströmung seit 1945. 
1 E* IE H. w undt: "Wärmeleitung i in pulsierenden Rohrströmungen vr 


u 


f 5 fi werden i in den nächsten Heften dieser Zeitschrift erscheinen. 


_Roferate der Vorträge: 2 ER ED } & 
> Fe Be AH yra: Zu neueren Problemen der atmosphärischen Strömung Aber‘ ein Gebirge: 
DE Frl Schu h: Zur Berechnung der furbulenten Grenzschicht bei eb 
‘ Druckverlauf | en: 


angen. 


sind solche Butter der RENTE ER an denen 
der von zwei beliebigen Normalen »,; ».EN ein- 
- ‚geschlossene Winkel nicht stumpf ist. Das skalare ; 
Produkt (v,».) ist dann nicht negativ. { DR 
 d=ule, 9,2) = u(P) sei eine skalare reelle zwei- 
‚mal stetig differenzierbare Funktion der Koordinaten 
; ‚2inB. - 
12 2 nehme u die Werte van: le) — “0, 


in beschränkte Funktion auf der Berandung. 
) Ableitungen werden abgekürzt gerchien 


hal 


2 ee 
6 AB; Aufgabenstellung 
zZ „der + Diterenigleichung in®: 


Au: — Una + Uyy' 2 U = 0 
.19* 
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mit der Randbedingung auf T': 

„t+ou=ß . Ion. (1) 
mit «(Q); ß(Q) als reell gegebenen Funktionen auf /, 
wobei &< 0 vorausgesetzt sei, sei eine Näherungs- 
lösung U(P) bekannt, die in ® zweimal stetig dif- 
ferenzierbar ist und der Differentialgleichung AU = 0 
genügt; dagegen erfülle U(P) die Randbedingung nur 
genähert: 


U(Q) sei beschränkt; durch U(P) und &(Q) ist dann 
die Funktion ß*(Q) auf der Berandung I’ definiert. 
Gesucht wird eine Abschätzung des Fehlers v in 
B-+HT: 

v(P)=U(P)—u(P).... (3). 


Hiernach ist v(P) zweimal stetig differenzierbar; es 
genügt der Differentialgleichung Av = 0. Nach (2) ist 
außerdem v(Q) beschränkt. 


C. Abschätzung des Fehlers 

Wir verwenden bei der Herleitung den Satz (4): 
Jede Funktion f(x, y, 2), die in einem beschränkten 
Bereich einer gewissen homogenen elliptischen Dif- 
ferentialgleichung 2. Ordnung, z.B. der Potential- 
gleichung genügt, nimmt ihren größten und ihren 
kleinsten Wert auf der Berandung des Bereiches an. 
(Siehe z.B. Collatz: Fehlerabschätzung zur 
l. Randwertaufgabe bei elliptischen Differential- 
gleichungen, Z. angew. Math. Mech. 32 (1952).) 

Auf IT’ genügt nach (1), (2), (3) die Fehlerfunktion 
v(P) der Randbedingung: 


+09 = (Ua, ta -W=P*—B. 


Aus den bekannten Funktionen &, ß, ß* läßt sich eine 
Funktion y(Q) auf T’ berechnen: 


„-E-f!n+s. 


& 
Damit lautet unsere Behauptung: 


Minys Min vs Max v<sMaxy. 
auftT ind+T in®ß+T auf! 


D. Beweis 


Auf der Berandung I’ nehme die Funktion %(Q) 
ihren größten Wert im Punkt Q, an; nach (4) ist dann: 


v(P) SPA); VA) SvA). 
Notwendige Bedingung dafür ist, daß für die par- 
tiellen Ableitungen nach dem Randparameter s, wenn 
sie im Punkt Q,(s,) existieren, gilt: 
sn us, + 2 — v8) <0; lim v8) — an _ > 0 


(5). 


‘ Auf der Berandung I' nehme die Funktion y(Q) ihren 


größten Wert im Punkt Q, an: 


a) SYr(Q). 


Wir führen unseren Beweis indirekt, indem wir die 
gegenteilige Annahme widerlegen. 
Annahme: Es sei das Maximum von 5(Q) größer als 


das Maximum von y(Q): 


iQ) >r@Q@)zrlaı) . . . 


Aus den Definitionen folgt dann: 


ra), 


1 
ZA) = YA) —5Q)<0 nach (6) 


und da nach (1) R </ 0 vorausgesetzt ist, muß sein: 


BHO SEOL 2: ga Ey 


Wir treffen nun die Fallunterscheidung: 

Fall a): Wenn im Punkt Q, die Normale » nicht ein- 
deutig bestimmt ist und das skalare Produkt von 
zwei möglichen Normalen (v,v.) = Oist, so sind diese 
Normalen »,, v, Tangenten an die. Berandung im 
Punkt Q,. Wählen wir als Randparameter z. B. 
sin Richtung »v,, so ist mit (7) eine Bedingung (5) 
verletzt und damit unsere Annahme (6) widerlegt. 

Fall b): Wenn im Punkt Q, die Normale » eindeutig 
bestimmt ist, oder wenn sie nicht eindeutig bestimmt 
ist und dann das skalare Produkt von je zwei be- 
liebigen möglichen Normalen (v,»,;) > 0 ist, so weisen 
diese Normalen in das Innere des Bereiches ®. 
Nehmen wir jetzt einen Punkt P, an, der vom 
Q, in Richtung der Normalen um die Strecke & ent- 
fernt im Innern von ® liegt, so gilt nach (3): 


vB) = 59) +5) + Zemr(Bo), 


wobei P,ein innerer Punkt der Strecke Q, P, ist. Durch 
hinreichend kleine Wahl von g, 


%(9ı) 
vr(P;) 
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erreichen wir also, daß B 

v(P) >%(Qı) 
ist im Widerspruch zu Satz (4). Damit ist unsere An- 
nahme (6) widerlegt. 

Durch Änderung der Vorzeichen von v und y 
läßt sich zeigen, daß Entsprechendes von den Minima 
gilt. Wir haben damit die behaupteten Fehler- 
schranken bewiesen: 


Miny < Minv < Maxv S Maxy. 


E. Hinweis auf die Möglichkeit einer 
Erweiterung des Anwendungs- 
| bereiches 

a) Bei n-dimensionalen Problemen treten an Stelle 
der Veränderlichen &, y, z die Veränderlichen z;; 
j=1...,n. Zugelassen sind Bereiche mit stück- 
weise glatter Berandung, die als Unstetigkeiten nur 
Knickstellen mit der Einschränkung 


ZU EN 
7 


hat. Hierin sind „ die n Komponenten der Normalen 
y, in der Knickstelle, 
b) Die Fehlerabschätzung ist anwendbar all- 


gemeiner bei elliptischen Differentialgleichungen 
2. Ordnung der Gestalt: 


L(u) = &%r Urjz; = z bjür, +cu=0. 


Inhomogene Differentialgleichungen können in üb- 
licher Weise durch Subtraktion einer speziellen 
Lösung hierauf zurückgeführt werden, ; 

Die Koeffizienten der Differentialgleichung unter- 
liegen den Einschränkungen: 


(ajx) pos. def., d.h. Zur viyk >0,.. .(8), 


worin vi die Komponenten eines willkürlichen Vektors 
sind, und 
c<o0,. 


Durch L(u) ist auf Z’ eine Konormale o definiert; ihre 
Komponenten berechnen sich aus den Komponenten 
der Normalen: 


ok = Zajkvi, 
; 


Nach (8) weisen die Konormalen bei glatter Be- 
rahdung in das Innere des Bereiches: 

Eokyk= Zajevivk>0. 

k j,k | 
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Allgemein-sind Bereiche mit stückweise glatter Be- 
randung zugelassen; die Berandung habe höchstens 
Unstetigkeiten ihrer Normalen mit der Einschränkung 


ad, 
zes in >0, 


worin v! die Komponenten der möglichen Normalen 
in einem Knickpunkt sind. 

In der Randbedingung und in der Fehlerabschätzung 
ist an Stelle der Normalen die Konormale einzuführen. 
Beim Beispiel der Potentialgleichung ist die Ko- 
normale gleich der Normalen. 

c) Wenn auf einem Teil des Randes T,, genannt 7}, 
als Randbedingung vorgeschrieben ist: 


u(Q) = 6(0) Qer, 


und auf dem übrigen Rand, genannt T;, die bisherige 
Randbedingung vorgeschrieben ist: 


&,(0) + au(Q) = B(Q) QE,, 


so ist für die Fehlerabschätzung die Funktion y(Q) zu 
definieren: 


für 


für 


v(Q) = 8*0) — 80) für GEN, 
Pr) -PA 5„. 
y(Q) = BETEN im OEL, 


worin 6*(Q) = U(Q) definiert ist. 

d) Zuletzt noch ein Hinweis auf eine Fehler- 
abschätzung, wenn auf einem Randteil 7‘, die par- 
tielle Ableitung 


% 


ger, 


vorgegeben ist. 

In einigen Sonderfällen des Differentialoperators, 
z.B. L(u) = Au, läßt sich nachweisen, daß die Fehler- 
funktion v entweder konstant ist oder ihren größten 
Wert nicht in einem inneren Punkt QET', auf dem 
Randteil I’, annimmt. 

Die Abschätzung läßt sich dann auf die obigen 
Fälle zurückführen. 


F. Ein numerisches Beispiel 
Das Problem 


Au(z,y)=0 für l.l<1, yI<1 
w—cgl’u=0 für - jej=1 
W—-U=0 für vI=1 

hat die strenge Lösung 
u=— cos 2 Cojy 


und etwa die harmonische Näherung 
U = — 0,356 + 0,184 (x? — y?) 

Die Fehlerabschätzung ergibt 

— 0,183 <U — us 0,063 
Die strengen Schranken sind 

— 0,0494 S U — u S 0,0281. 
Der Fehler wird hier also etwa 3fach überschätzt. 
Hannover. H.J. Grünsch. 


_ Über ein Beispiel aus der Statistik. 


Die hier betrachtete Fragestellung kann man kurz 
so kennzeichnen: X,, X,, ..... seien unabhängige nach 
einer Normalverteilung mit Mittelwert 0 und Streu- 


ung o? verteilte zufällige Variable (kurz: nach 


N (0, 0?) verteilt). Welche Verteilung hat die Summe 
der Abweichungen vom Mittelwert & X, oder auch 


iSE, 
. die Summe der Fehlerquadrate 2 X?, wenn &un- 
i< 
abhängig von den X, nach einer geometrischen Ver- 


teilung mit dem Parameter g verteilt ist, d.h. 


VEe=n)=-pg(n=12..) p+qg=1 
Außerdem sollen Fragen der Parameterschätzung 
für q und auch für o? erörtert werden. 

Diese Problemstellung scheint in dieser speziellen 
Fassung noch nicht behandelt worden zu sein, jedoch 
sind in der (mir nur unvollkommen zugänglichen) 
mathematisch-statistischen Literatur Verteilungen 
des hier betrachteten Typus hinsichtlich gewisser all- 
gemeiner Eigenschaften untersucht !)?2), Die oben 
aufgeworfenen Fragen gehen auf Unterhaltungen mit 
Statistikern der Praxis zurück, worin folgendes 
äquivalente Problem aufgeworfen wurde: Es laufen 
zwei Produktionsgänge I und II, die zwar in wahr- 
scheinlichkeitstheoretischem Sinne als unabhängig 
angenommen werden können, jedoch in bestimmter 
Weise miteinander gekoppelt sind ®). Die Ergebnisse 
von I: &,&,,... seien identisch alternativ verteilt mit 
den Wahrscheinlichkeiten 9, ,, p+q=1 und etwa 
p >gq. Der Prozeß II lasse sich charakterisieren durch 
zufällige Variable X,, X,, ..., welche unabhängig 
nach N (0, 0?) verteilt sind und es seien weiter €, und 
X, aufeinander bezogen. Gefragt wird nach der Ver- 


N N 
teilung von z{D= % X, oder auch von 2@) = N X}, 
1 i=1 
wenn n das erstmalige Auftreten des unwahrschein- 
licheren Ereignisses in der Reihe der &; anzeigt. Außer- 
dem sei es als untunlich angenommen, den Prozeß I 
direkt zu beobachten. Es mögen nun eine Reihe von 
Ergebnissen hinsichtlich der Größen z{1) oder x(2) vor- 
liegen und daraus Schlüsse über die Parameter 
q, 0<gq<1l oder o?, 0< 0?<{ oo gezogen werden, 
genauer Schätzfunktionen oder Konfidenzintervalle 
angegeben werden. 


15 aD, a, 19 0) seien unabhängige, zufällige 
Größen, die einer mit der Dichte f(z() verteilten 


N 
Gesamtheit entnommen sind. Hierheiist «dD= N X,, 
i=1 
wobei die X; unabhängig nach N(0, 0?) und die n un- 
abbängig davon geometrisch verteilt sind, so daß 


Wna=N=pAgp+gqg=1. 


Wir zeigen: Konsistente Schätzfunk- 
tion für gist: 


Hierzu sei bemerkt: Die Verteilung von x() ist durch 
©o 1 „2 
fizid) = = yPTcme-n<all)<o (2) 
2 


TOn=1 /)R 


gegeben. Damit fängt man nicht viel an, doch erhält 
man leicht für die charakteristische Funktion 


en 


gu = 


—o 


1) Feller, Ona General Class of, ,‚Contagious‘‘ Distributions, 
Ann. Math. Statistics 14 (1943), p. 389—400, zit. in Feller, An 
Introduction to Probability Theory and ist Applications I. New 
York — London 1950. 

2) Robbins, The asymptotic distribution of the sum of & 
random number of random variables. Bull. Amer. Math. Soc. 54 
(1948), S: 1151—1161. 

2) Es wurde z.B. die Frage aufgewırfen, ob es möglich sei, die 
'Wahrscheinlichkeit für Stockungen in einer Strangmaschive bei 
‚der Zigarettenerzeugung an dem unschwer festzustellenden Aus- 
schuß zwischen je zwei Stockungen zu schätzen. Tatsächlich 
besitzen jedoch die modernen Maschinen dieser Art stets Zähl- 
einrichtungen so daß eine direkte Schätzung ebenso leicht durch- 


führbar ist, ; 
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so daß man für f(x(1)) auch 


+% 
IT, 
1 ge 2 - ital) 
1-26 
erhält. 
Natürlich ist Z(z(d) = 0 und 
8 
E(zW2) — = DES (4), 
3 0% 
N) herr ta 


wie man sofort durch Differentiation von (3) oder 
kürzer durch Aufsuchen der Koeffizienten von t? bzw. 
{4 in der Potenzreihenentwicklung von (3) findet. 


mM 

I) 
Sei U) = en . Dann ist E(#(1)) = 0 und nach (4) 

2 

E(zU2) = — . Es ist zwar prinzipiell möglich ver- 
möge (3) die Verteilung von z{l) aufzuschreiben, aber 
für die praktische Verwertbarkeit bei der Parameter- 
schätzung kommt dieser Ausdruck kaum in Frage. 
Aus eben diesem Grunde dürfte sich auch die Ver- 
wendung des Maximum Likelihood Prinzipes ver- 
bieten. 

Nimmt man zunächst an, daß o? bekannt ist, bleibt 
als letztes Auskunftsmittel immer noch die Tsche- 
byscheffsche Ungleichung: Sei 0< a <{ 1 vorge- 
geben: Es ist 


c 1 
#(lzW| < ——  —— & 
( = VYmg yı —)E 
und daraus erhält man zu vorgegebenem Konfidenz- 
koeffizient >« ein Vertrauensintervall für g: 


0? 
ISIS zum ah) 


Für hinreichend großes m und bekanntes 0? erledigen 
sich die Fragen der Schätzung von q natürlich durch 
den Grenzwertsatz. 

Um für q und o? Schätzfunktionen zu gewinnen, 
kann man ausgehend von der Pearsonschen 
Momentenmethode etwa so verfahren: Es ist nach (4) 
und (4) 


B> „(2 
Nach bekannten Sätzen konvergiert daher ER 
m 


Zal4 
1 


i co 301 
stochastisch gegen — und gegen —-(2—9). 


Aber auch jede rationale Funktion dieser Größen 
konvergiert nach Slutsky stochastisch gegen 
die entsprechende Funktion der Grenzwerte 4), 


Führt man also durch 
PEULEER 
De, 


Zalllı 30 A 
= ir —4q) 


und 


*) Vgl. Oramö6r, Mathematical Methods of 8 - 
Ani Heabn. of Statistics, Prin 
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Schätzfunktionen 0°? und Z ein, dann ist 
m Z a{1)4 
(Ze) 


und damit schon erwiesen, daß (1) konsistente Schätz- 
funktion für q ist. Wir verzichten darauf den Aus- 
druck für 9? hinzuschreiben. 


2. Es ist klar, daß die Fragestellung für nach (2) 
verteilte Gesamtheiten größere praktische Bedeutung 
hat, als wenn man die Verteilung von x(2) zu Grunde 
legt. Es zeigt sich jedoch, daß man die Untersuchung 
dieses Falles (insbesondere nach Anbringung einer 
leichten Modifikation) unter Heranziehung schon 
längst bekannter Prüfverteilungen durchführen kann. 
Weiterhin sei o® als bekannt angenommen, etwa 
o= 1. 


2-93= 


Es seien also die Größen X, unabhängig nach 


N(0,1) verteilt und n davon unabhängig geometrisch 
und es wird die Größe x(2) = £ X} betrachtet. Für 
ihre Verteilung ergibt sich: er 

_z(2)j2 


qp[® —1/2 
ir (z(2)b) Er 


„(2) 
b (d—1/2) z(2) et 
flr9) = +y2 e | yi dt+ \(6). 
0 
ee n ; z>0(D & 2) 
0; z22)<0O 


Die charakteristische Funktion ist 


VYI-2it+p 
(2 we, RER AFFE 
p(2) (t) In pa 247 
der Mittelwert Z(x(2)) = re die Streuung 
2(29 ROTE 
q q 


Hierzu beachte man: Da die X, unabhängig nach 
n 

N(0, 1) verteilt sind, ist z(2) = 5X? bei festem n 
i—1 


= 
bekanntlich nach Helmert-Pearson mit n 
Freiheitsgraden verteilt. Dies liefert für die Dichte 
von x(), wenn nun n unabhängig von den X; als 
geometrisch verteilt angenommen wird: 


| Se! „ealr | 
fir) = 2 PILFUTTE prlg 2M>O (8). 


0 2) <0 


Die Übereinstimmung von (6) und (8) folgt lei 
der Bemerkung, daß (6) (8) folgt leicht aus 


A 


el? (xd) 1% 
y=e a  , 
a1 T'(n/2) 


71/2) 
n 
wegen I' (2) -(7 _ 1) r(z _ ) der einfachen 


Differentialgleicbung genügt: 


#2 —1/2 
ER 2]2 (zb) +4)» 


T(17) P) cd. 


EM ITE 


u nt ETLEE EEURETETN Te 


A 
a 


(10) sofort: 


für # >0 durch 
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Löst man (9) unter der Bedingung, daß das von 
0 bis oo erstreckte Integral über (8) den Wert 1 erhält, 
ergibt sich (9). 


. % 
Es seierwähnt, daß man F(x) = f f(z(2)) dz(2) leicht 
0 


numerisch berechnen kann, wenn eine Tafel für die 
unvollständige Gammafunktion für die Parameter Us 
und 1 vorliegt sowie eine Tafel für das Ga u ßsche 
Fehlerintegral und die Dichte der Normalverteilung. 


EI) 
Wegen E (2%) ee 


yı-2it 


= 1 
Bern 
pt) ara am? q 


erhält man (7) sofort 


Liegt nun eine Stichprobe vom Umfange m hin- 
sichtlich der Größe x(2) vor, d.h. m unabhängige 
nach (6) verteilte Größen x”), ..., x), dann kann 
man bei Verzicht auf die Kenntnis von (1) jedenfalls 
wieder mittels der Tsehebyscheffschen Un- 
gleichung Vertrauensgrenzen für q gewinnen. 

Man kann jedoch mehr im Anschluß an die heute 
üblichen Prüfverfahren Parameterschätzung vor- 
nehmen, wenn man die Zuordnung zwischen den 
Prozessen I und II etwas abändert. Hierzu setze man 
den ersten Produktionsgang wie oben voraus, für den 
Prozeß II seien die Variablen X,, X,,... unabhängig 
nach N(0, 1) verteilt und jedem £&, aus I seien jetzt 


zwei Variable X2x—-ı, X2x zugeordnet. Die Frage- 


stellung heißt also jetzt zunächst: Wie lautet die Ver- 
teilung von 2n Fehlerquadraten, deren Variable un- 
abhängig nach N(0, 1) verteilt sind, wenn n davon 
unabhängig geometrisch verteilt ist ? 

Es ist ganz einfach zu zeigen, daß unter diesen 


; 2n 
Voraussetzungen 2 = Y x; nach ; 
i-1 


e>o| 


ae 
== ee 


verteilt ist, BIEN, 
Weiterhin hat nun x nur noch die hier festgelegte 


Bedeutung. ? 
Um (9) zu bestätigen, bemerke man nur, daß g(z) 


o 2 nl gel? _®D 
a —1 
lm) P 


gegeben ist. 
g(x) ist also eine Pearson -TypIII-Verteilung, 
_ somit ge nach Helmert-Pearson mit 2 Frei- 
 heitsgraden verteilt. Dies darf vielleicht deswegen 
gewisses Interesse des Statistikers beanspruchen, 


- weil an bei der Beobachtung von x im Mittel 


eine große 


ahl von Freibeitsgraden ‚‚verlorengehen‘“. 
Man hat ja im Mittel-2/9 Beobachtungen zu erwarten, 
um den aufgetretenen Wert von x zu erhalten. 


Es soll nun eine Stichprobe vom Umfange m von 
Größen %,,... ., m vorliegen, die alle nach (10) verteilt 


sind. Für die Verteilung von 


tent tm 


erhält man wegen der Reproduktionseigenschaft von 


(a\m em! Ag ? wre 
(2 a al). 


Wie man leicht kontrolliert, erhält man ET 
teilung auch dadurch, daß man die Summen zZ 2 


ct - = . uB Be BEE 

betrachtet und die Anzahl der Summanden n nach 

iner (etwas modifizierten) P as calschen Verteilung 
Sn ö EN E r .2 E 


= u 


Ba Fa £ f Ben er Re 


unabhängig von den X; verteilt ansieht, für die den 
Merkmalen 0, 1,..., m — 1 die Wahrscheinlichkeit 0 
zukommt. Dies muß so sein, da sich die so definierte 
Pas calsche Verteilung als Summenverteilung von m 
unabhängigen der hier betrachteten geometrischen 
Verteilungen ergibt. 

Um nun Konfidenzgrenzen für q zu gewinnen, steht 
eine Reihe von Möglichkeiten offen, die auf bekannten 
Prüfverteilungen beruhen und nur kurz gestreift 
werden sollen: 

a) Der Vollständigkeit halber erwähnen wir auch 
hier die Tsehebyscheffsche Ungleichung. 
u 


Sis- 
Fenll2), 


Ein Vertrauensintervall für q zum Koeffizienten >« 
(0<a<1) wird somit 


2 Rn 
= —_. ee ; 
: Im(1—a) 2 Ym(1—a) 

b) Engere Konfidenzgrenzen sind zu erwarten, 


wenn man (11) heranzieht und %,, x, so wählt, daß 
Xı 


N ne ea 


Be: 
Im) 2” 
Kı=" 
Dann erhält man aus Wu, <qg&<x,)=a oder 


W ER <gq< X2\ x ein Konfidenzintervall zum 
mi M% 


Koeffizienten a. 

q& — 2m 
t —— 

2/m 
N(0,1) verteilt und daraus erhält man ein asymp- 
totisches Vertrauensintervall 
md 
m 5a 


<a)=s, 


Für großes m is asymptotisch nach 


q& — 2m 
V2 m 


w(-4< 
wenn ; 


era. 


1 d 
il 


 d geht wesentlich langsamer gegen © für a —1 als 


ner Eür& = n gilt z.B. d = 1,96 und für 
—& 


den Wert aus a) Bash 4m) ; 
3 


c) Die in b) erhaltenen asymptotischen Konfidenz- 
intervalle sind auch die asymptotisch besten, wie man 
erkennt, wenn man die Likelihood-Funktion der 
Stichprobe #, ..., zm heranzieht, d.h. 


q m _ tr tm) 
2) = (3) e 2 RE 


IE =. 
\m 


.(13) 


GE 

Weiter ist 

Blog Di er N loop D\ 1. 
An -0q \= me? 2 ( 9 ) mg 

Bekanntlich ist aber ae 2 asymptotisch 

ER ; öq m € 


nach N(0, 1) verteilt. Diesliefert natürlich für m > 


dieselben Konfidenzgrenzen wie in b). 


i 
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Von (13) ausgehend kann man im Anschluß an 
Cramör?°) für jedes m gültige Konfidenzgrenzen 
für q gewinnen. ENT 

Das Maximum Likelihood-Prinzip liefert 
nämlich die Schätzfunktion 9 = 2/#. Vermöge (11) 
rechnet man leicht nach: 


mg 
= 


2(m —]) 
mE 


d.h. = 


Eo)=- 


ist unverfälscht. 
Für die Streuung erhält man: 


2 
Eo— = 1,5 (m>3). 


Aus (11) gewinnt man durch Einführung der neuen 
Variablen z = 1/2 die Verteilung von 1/2. Betrachtet 
man nun die standardisierte Variable 


2(m — 1) — qmz& —; 
wi 5 me 2 
„ qam“E fi i 
dann ergibt sich nach leichter Rechnung die folgende 
parameterfreie Verteilung: 


ho) = 
| 1 KMm-Y(m-23, 
—) Tim) + ma"! 
— Ym—- 2< v4 =) 
0 v<—\m— 2 (14). 
Nun wähle man v, %, >—\m— 2 und so, daß 
120) dv=x. Dann erhält man aus W(v, <v<v,)—=& 


v+/m-2 


(m—1) Vm—2 | 


m>3 


v, 
das folgende Konfidenzintervall für qg zum Koeffi- 
zienten &: 


Am-h/m-2_, 2(m—1)/m —2 
mä(lu +Ym—2) mil, +VYm— 2) 


3. Bekanntlich kann man die Poisson-Ver- 
teilung als Grenzverteilung Pascalscher Ver- 
teilungen gewinnen. Wir. wollen dies zum Anlaß 
nehmen, um noch ähnliche Fragen wie eben zu 
streifen: X; seien unabhängige nach N(0,1) ver- 
teilte Größen, n sei davon unabhängig nach einer 
Poisson- Verteilung verteilt mit W(n = 0) =0 
und dem Parameter a, Gefragt ist die Verteilung der 
Summe einer geraden Anzahl von Fehlerquadraten: 


2n 2 
Y = yX; . 
i=1 


Man erhält: 


oo vr el? ak-ı ee 
are. | 
fy) = gyl2 ay .. (15), 
et lfer)e>i 
0 .9<o0 


wobei I, wie üblich die Bessel- Funktion rein 
maginären Arguments bezeichnet. 


Die charakteristische Funktion für die hier be- 
trachtete Poisson -Verteilung ist durch 
6% git guet 
gegeben, also erhält man für die Fouriertrans- 
ormierte von (15): 
1—2it 


Daraus: 


Biy)=2(1+0), B(y—2(1+a))=4(1+ 2a), 
®) L.c. *) p. 505. 


Yu. +», Ym sei eine Stichprobe aus der nach (15) ver- 
teilten Gesamtheit. Für die Schätzung von a bei 
kleinem m stehen z. B. folgende zwei Wege zur Ver- 
fügung: : 2 
Füra > »ist nachBRobbins?)y asymptotisch 
nach N (2(1-+ a), 41 + 2a)) verteilt, also kann für 
Yıt9y+t:':t9m 


mit 
ei 


hinreichend großes ay = 


N(2(1-+a), EN + a) verteiltangenommen werden. 


Für kleines «>0 steht jedenfalls wieder die 
Tschebyscheffsche Ungleichung zur Verfügung 

Die Kenntnis von (15), (16) und damit von der 
exakten Verteilung für Y dürfte für Zwecke der 
Schätzung von a nicht von Belang sein. 

Wien. L. Schmetterer. 


Eine graphische Richtungsfeldkonstruktion 
für den Phasenplan nichtlinearer freier 
Schwingungen. 

Da für nichtlineare Schwingungen — es soll sich 
dabei um Systeme mit einem Freiheitsgrad handeln — 
im allgemeinen keine Aussicht auf eine geschlossene 
Integration der beschreibenden Differentialgleichung 
& — f(x, z) = 0 besteht, so begnügt man sich viel- 
fach mit einer qualitativen Diskussion des Vorganges 
im „Phasenplan‘“, Hier wird in einem cartesischen 
Koordinatensyttem die Geschwindigkeit 
& = yin Abhängigkeit vom Ausschlag z (anstatt 
von der Zeit t) betrachtet, so daß die Differential- 


gleichung der „Phasenkurven“ SU — say) lautet. 
Diese Behandlungsweise wird als „‚‚topologische 
Methode‘ bezeichnet. Wie man nun unmittelbar er- 
kennt, entsprechen densingulären Punkten des 
Phasenplans y =f(z,0) = 0 (verschwindende Ge- 
schwindigkeit und Beschleunigung), die also auf der 
x-Achse liegen, die Gleichgewichtslagen 
des Systems (deren Stabilität oder Instabilität mit 
dem geometrischen Charakter der Phasenkurven in 
der Nachbarschaft der Singularitäten zusammen- 
hängt), während einer geschlossenen Phasenkurve 
ein periodischer Schwingungsvorgang, oder einer 
spiraligen Phasenkurve eine gedämpfte bzw. an- 
gefachte Schwingung entspricht. Für einen orien- 
tierenden Überblick über die Gesamtheit der Phasen- 
kurven (Phasenplan) und damit über den Schwingungs- 
verlauf bei allen möglichen Anfangsbedingungen im 
Falle einer allgemeinen Dämptıngehunktioh p(x) und 
einerlinearen Rückstell (Feder-)kraft hat sich die 
graphische Konstruktion für das Richtungsfeld der 
Phasenkurven von A. Li&nard eingebürgert !), 
Im folgenden wollen wir eine andere und, wieYuns 
scheint, recht einfache Konstruktion für den 
viel allgemeineren Fall einer nichtlinearen 
Rückstellkraftundeinsskombinierten 
DEUPIDEBEESREANTE angeben. Ergänzend 
seien hierzu noch die Verfahren vonH.Heinrich?) 
und W. Richter?®) erwähnt, die als Hilfsmittel 
Nomogramme benötigen; es ist indessen fraglich, ob 
für den vorliegenden Zweck dieser Aufwand lohnt. 
Ferner sei für den Fall erzwungener Schwingungen, 
wo auf der rechten Seite der Differentialgleichung eine 
explizit zeitabhängige Funktion steht, etwa auf das 
Verfahren von E. Braun) hingewiesen, 

*) A. Li6&nard: Btude des oscillations entretenues. Rev. 
G6n. d’ Elect., 1928. — Vgl.hierzu die modernen Werke über nicht- 


lineare Schwingungen; etwa J. J. Stoker: Nonlinear Vibrations, 
Interscience Publishers, Inc., New York 1950. 

») H. Heinrich: Ein abbildungsgeometrisches Verfahren zur 
Darstellung von Richtungsfeldern und zur Erzeugung von Tan- 
gentenverwandtschaften. Dtsch. Math. Bd. 3 (1938), 8. 353—389. 
ae a = Ir hy a N nn BDIeBBIDEwERung in einem 

o els nomographischer fsmittel. . 
Bd. 11 (1940), 8. 437—450. nn er: 

*) E. Braun: Über die graphische Lösung der Differential- 
a SB ee bei beliebigem Gesetz 

pfung, ckste und Antriebskra „Archi 
Bd. 8 (1937), 8. 198—202. I 5 E 
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Wir befassen uns zunächst mit Schwingungs- 
gleichungen vom Typus 


2+ pi) +Fla)=0 [Fa)=atz +ba® +] (la) 


Die der Rückstellkraft proportionale Funktion F(x) 
wird bei den Anwendungen meist durch ein Poly- 
nom approximiert oder ist eine Kreisfunktion. Als 
Dämpfungsfunktionen können vorkommen die ab- 
schnittsweise analytischen Funktionen (3) = 
(sgn &) » r (Coulombsche (trockene) Reibung), p(#) = 
(sgn &) r&? (quadratische Dämpfung durch Luft- 
widerstand) oder bei selbsterregten Schwingungen 
z.B. o(&) = — 0? + ß?%° (Anfachung bei kleinen, 
Dämpfung bei großen Geschwindigkeiten). Wird in 
(la) @ = y gesetzt, so hat man für das Richtungsfeld 
der Phasenkurven die Beziehung 


dy _— poly) — Fa) 
= een 


Die singulären Punkte $ (y = 0)sind durch — F(x) = 
9(0) bestimmt. 

Zur Ausführung der Konstruktion (vgl. Bild 1), für 
diesich dieVerwendung von Millimeterpapier empfiehlt, 
werden vorerst die „Grundkurven“ y = — F(x) und 
y = p(x) eingetragen (am einfachsten mit der gleichen 
Zeicheneinheit für beide Koordinatenachsen). Dem 
zweckmäßig im Schnittpunkt zweier Netzgeraden des 
Millimeterpapiers gewählten jeweiligen Feldpunkt 
Pı(%,Yı) wird die „Trägergerade‘“ (Netzgerade) 
x = y, zugeordnet, welche die p-Kurve in Q, schneide. 
Die Netzgerade x = x, durch P;, schneide die F-Kurve 
in R,. Durch Projektion von Q, auf die y-Achse erhält 
man als Fußpunkt den ‚Pol‘ Y,, durch Projektion 
von R, auf die Trägergerade den Fußpunkt X,,. Die 
Richtung der Verbindungsgeraden 
YıXı ist dann gleich der gesuchten 
örtlichen Feldrichtung. Das ist un- 
mittelbar aus Bild 1 abzulesen, da Q,Xı= 
— F(z,) — plyı) und 7,9, yı- 

Will man sich nun ein genügend dichtes Richtungs- 
feld. verschaffen, so faßt man alle Punkte P,, mit 
gleicher Ordinate y, zusammen, denen der gleiche 
Pol Y, ‚und die gleiche Trägergerade x = y, zu- 


- Jay 


# 
| 
Pa... En en Pu (19) y- yo) 
ig ' 167 
ENDET SE NE RE & 
s 
dy\ _ Pu) 
Kalb 57 


y--Fik) 
Bild 1. 


geordnet sind. Beim Übergang zu Punkten P,, mit 
gleicher Ordinate y„ ist im wesentlichen nur der 
neue Pol Yzu ermitteln, während man die Punkte X „x 
durch Herüberloten der bereits bekannten Punkte X,, 
auf die neue Trägergerade x = y„ sofort bekommt. 

Wir gehen nun zu Differentialgleichungen von der 
allgemeineren Form 


+) Ha) +Fa)=0..... (28) 


über, bei denen also Größe und Vorzeichen des 
Dämpfungsgliedes außer von der Geschwindigkeit 
auch vom Ausschlag abhängig sind. Von praktischer 


Bedeutung ist allerdings bislang wohl nur der Fall‘ 


2) = 4; ist ferner z.B. @(@)=— oa’ + 3 ß?a*, so 
t man es wiederum mit selbsterregten Schwingungen 
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zu tun. Hingegen ist der Fall 9(&) = (sgn 5) : #? hier 
von untergeordnetem Interesse, da er sich mit 
elementaren Mitteln erledigen läßt: die Substitution 
Fe 2 führt auf dielineare Differentialgleichung 
1 de 
Quadratur liefert t = {(x). 
Setzt man in (2a) x = y, so ergibt sich 
d — G(2)- —F 
z £ I N (2b). 
T Y 
Die singulären Punkte 8 sind durch — F(x) = 
G(z) : 9(0) gegeben; ist insbesondere 9(0) = 0, so sind 
sie die Nullstellen von F(x). Der Richtungsfeld- 
konstruktion (vgl. Bild 2) dienen als Grundbestand- 


und eine zusätzliche 


teile die Kurven y= — F(x), y=@(x), 2 =o(y) 
y-! | 
x= Yly) 
y=6(x) 
—————r 


| = Surf 


Pı 4 


Ya -F(x) 


und die Ordinatenparallele x = 1 (Zeicheneinheit). 
Als jeweiliger Feldpunkt P,,(zı, y„) empfiehlt sich 
natürlich wieder ein Knotenpunkt des Millimeter- 
netzes. Dann schneide die Netzgerade x = r, die 
G-Kurve in Qı und die F-Kurvein R;, die Netzgerade 
y = Yu die o-Kurve in T,„. Durch Projektion von Qı 
auf die Ordinatenparallele © =1 erhält man den 
Fußpunkt Z,, durch Projektion von R; auf die Träger- 
gerade x = y„ den Fußpunkt X,,. Verbindet man den 
Ursprung O mit Er, so schneidet die Verbindungs- 
gerade die Ordinatenparallele durch 7,, die im 
Spezialfall o(&) = & = y mit der Trägergeraden zu- ' 
sammenfällt, im Punkt 7’, (dessen Lage übrigens 
auch im allgemeinen Fall leicht zu schätzen ist). Be- 
stimmt man noch durch Projektion von T7/,, auf die 
y-Achse den Fußpunkt Y.., so stimmt die 
örtliche Feldriehtungin P,, mit der 


Richtung von YuaÄyu überein. Nach 
dem Strahlensatz.ist nämlich (vgl. Bild 2) 


E38,:0E} = TOT: OT) 


oder FR, 
OYur — Ri == G(z) "olyu) 


womit die Konstruktion auf die des Gleichungstyps 
(1) zurückgeführt ist. 
Ist speziell F(x) = wie: 

" to8%)+olr=0,......(), 
so sei angemerkt, daß sich der Sonderfallg($) = % der 
Gleichungsform (2a): 

ty) Itoey=0 ....(W 
sofort auf (3) zurückführen läßt; denn substituiert 
man y=&: 


+ 9G)z toi, 
‚so liefert, wenn g(z) = fo dt gesetzt wird, die 


Integration in der Tat 
ne (3). 


« 
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(Die Integrationskonstante c kann Null gesetzt werden 
da die Differentialgleichung die Transformation 
2 +c=7 zuläßt; dies bedeutet physikalisch, daß 
der Ruhelage der Nullausschlag zugeordnet wird.) 
Die formale Analogie darf jedoch nicht übersehen 
lassen, daß es sich bei «(t) und y(t) natürlich um ver- 
schiedenartige Schwingungsabläufe handelt. Die 
Differentialgleichung (3) ist aber bequemer zu dis- 
kutieren als die entsprechende (4), was nicht immer 
ausgenutzt wird. In dem von Lienard graphisch 
behandelten Spezialfall (3) tritt bei unserer Konstruk- 
tion noch die Vereinfachung ein, daß die Koordinaten 
der Punkte X, auf der Trägergeraden durch & = yu, 
y = — x, unmittelbar gegeben sind. Ist in (4) z. B. 
Ö(y) = — a®+3 ß?y?, so liegt die bekanntevan der 
Polsche Gleichung vor. 


Durch die äargelegten Verfahren kann man einer- 
seits einen raschen Überblick über die möglichen 
Schwingungstypen erzielen, andererseits aber auch 
genauer diejenigen Schwingungen konstruieren, für 
die man sich im besonderen interessiert. 


Göttingen. Manfred Schäfer. 


Der Energletransport in strömenden Medien 


Jede Strömung eines flüssigen oder gasförmigen 
Mediums ist mit einem Energietransport verbunden. 
In der Elektrodynamik hat sich seit langem zur Dar- 
stellung des Energietransportes im elektromagneti- 
schen Feld der Poyntingsche Vektor eingebürgert. 
In der Strömungslehre fehlt, soweit dem Verfasser 
bekannt ist, ein analoger Vektor. Nur für den Sonder- 
fall der reibungsfreien Flüssigkeit mit konservativen, 
zeitunabhängigen Massenkräften wurde ein solcher 
Vektor hergeleitet }). 


Es soll im folgenden, zunächst für reibungsfreie, 
inkompressible Strömungen, ein Energietransport- 
vektor © nach einer Methode hergeleitet werden, die 
eine Erweiterung auf beliebige Strömungen eines 
zähen und kompressibeln Mediums zuläßt. 

Ist do der Vektor eines beliebigen, raumfesten 
Flächenelementes, so soll: 

dE = 6©-dodi?) ll), 
die im Zeitelement di durch die Fläche do trans- 
portierte Energie der Strömung bedeuten. Damit ist 
auch die Bedeutung des Energietransportvektors fest- 
gelegt. 

Wir wollen uns im folgenden grundsätzlich auf die 
mechanischen Energieformen beschränken. Es sind 
dies die Bewegungsenergie der strömenden Flüssigkeit 
und die potentielle Energie der Massenkräfte. Dazu 
tritt noch die mechanische Arbeit, die beim Durch- 
tritt durch do in der Zeit dt von der durchströmenden 
Flüssigkeit nach außen geleistet wird?), 


Es ist demnach, wenn v der Vektor der Strömungs- 
geschwindigkeit, 0 die Massendichte, p der statische 
Druck und U das auf die Masseneinheit bezogene 
a ebupngige) Potential der Massenkräfte be- 

eutet: 


= (Ste: O)o-doat + pao- nat 


_ [po ; 
-($+ | +O)ev do dt. 


ı) Vgl. Cl.Schaefer, Einführung in die theoretische 
Physik, Berlin und Leipzig 1929, I. Bd., 8. 796 ff. 

%) Der Punkt bedeutet das skalare Produkt. 

®) Eine strenge Herleitung des Energietransportvektors 
für den allgemeinen Fall der zähen, kompressibeln Flüssig- 


keit aus den Stokes-Navierschen Gleichung soll als Er- 
gänzung dieses Aufsatzes folgen. 


Durch Vergleich mit Gl. (1) folgt daraus: 
e-[Pr PE (2). 
0 2 


Bei der reibungsfreien, inkompressibeln Strömung hat 
also der Energietransportvektor dieselbe Richtung wie 
der Geschwindigkeits- oder der Impulsvektor und 
man erhält ihn, wenn man den auf die Volumseinheit 
bezogenen Impulsvektor mit dem Betrag der Gesamt- 
energie pro Masseneinheit multipliziert. 

Der Energietransport durch eine geschlossene 
Fläche ist gegeben durch: f &©-dodt. do soll nach 

F 


außen weisen. Bezieht man ihn auf die Volums- und 
Zeiteinheit und schnürt die Fläche auf einen Punkt 
zusammen, so erhält man: 


JS ©: do dt 


F 
lim ep 
v0 vd E 


Man erhält also die Divergenz des Energietransport- 
vektors. * 

Aus der Eulerschen Gleichung der Hydro- 
dynamik in Verbindung mit der La m bischen Um- 
formung des substantiellen Differentialquotienten 
folgt aber: 

ou wi ? , Jol 
2 x(pxp) =— Ar; + 2) (3). 


Aus Gl. (2) folgt ferner: 


ARE Ars 


2 
+{@+ 3 +0)r-en. 


Da aber o = const. vorausgesetzt wird, ist 


vr .(eY)=0 
Es folgt somit aus Gl. (3): 


° 
7-6 =|- 5 + xwxu| en. 
Nun ist aber: 
[px(pxv)]-oev=0. 


7 SR RTEINN 
V s--z( a Race (4). 


Die aus dem raumfesten Volumselement pro Zeit- 
und Volumseinheit austretende Energie entspricht 
also der Verminderung der auf die Volumseinheit be- 
zogenen Bewegungsenergie in der Zeiteinheit. Die 
potentielle Energie er Änderung erfahren, da 
das Potential als zeitunabhängig vorausgesetzt wurde. 
Bei stationärer Strömung ändert sich die Energie 
natürlich nicht. 


Geht man zu zähen Flüssigkeiten über und nimmt 
zunächst noch 0 = const an, dann lautet bekanntlich 
der Spannungstensor p, derdem Stokes-Navier- 
schen Ansatz zu Grunde liegt: 


=—-Ip+tnlmo+oup) (8). 


% bedeutet den Einheitstensor, der Beistrich soll das 
dyadische Produkt andeuten, n ist die dynamische 
Zähigkeit. Die Arbeit, die beim Durchsetzen des 
Flächenelementes do in der Zeit dt nach außen ge- 
leistet wird, ist dann: 

dA=—(p-do).vd ....(6). 
Aus (5) ergibt sich aber: 


dA =pv-dodt — [n(p,v +p,p)]: do dt. 


Daher gilt: 
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Wegen der Symmetrie der Dyade p, v + v, ygiltnun: 
(WB +d,p)-do-vd=(p,o-+ov,p)-v-dodt 
und daher: - 
dA=[pv—n(p,v+p,p)-v]-dodt (6a). 
' Die mechanische Energie dE, die durch die Fläche do 
während dt transportiert. wird, ist demnach: 
dE — (5 be+ed) v-dodt +dA. . (7), 
wenn wieder ein zeitunabhängiges Potential der 
Massenkräfte vorausgesetzt wird. 


Aus (6a) und (7) folgt: 
ang + 2 +o)eu-nn,v+on:n]-aoa: 
und daher nach I. (1): 
o-(2+, +)en-nuinnton). . (8). 


Der Energietransportvektor hat nun im allgemeinen 
nicht mehr die Richtung der Stromlinientangente, 
sondern schließt mit ihr einen von Null verschiedenen 
Winkel ein. 2 

Bildet man wieder die Divergenz von ©, so erhält 
man für n = const ; 


a rl 
v:6=v(gtg +0):eo I% 


| . (9). 
nv le (vd +0, p)] 
Nun lautet die Stokes-Naviersche Glei- 
chung in Vektorform: 
Ob /\vl? 
atr(e) —bdx(pxb) 
—- A er *(10) 
eu Are ; 
mit - 
N 
OT 10a 
Hr (10a) 
Führt man Gl. (10) in (9) ein, so erhält man: 
er ö /o|vl? 
rv: = (5 mern enorm) 


: Man kann nun folgende Umformung ausführen; 
7 To (p, 045, M)]= (9) (040, W)+5-7°% 


"wobei die Beziehung p - » = 0 beachtet werden muß 
_ (-- bedeutet das doppeltskalare Produkt der beiden 
Dyaden). Daher folgt: R 


2.96 - = u a) (ot). 


z & r Zum Verlust an Bewegungsenergie tritt nun noch die 
- Energiedissipation 

np ne 
Man überzeugt sich leicht, daß diese Formel mit den 
in den Lehrbüchern angegebenen Formeln für die 


pro Raum- und Zeiteinheit vom 


X Energiedissipation übereinstimmt ®). 


1931, 8.6541f. 


w Per = . 
ee 4 or x er nn DR 
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Betrachten wir z.B. die Poiseuillesche Strö- 
mung in einem zylindrischen Rohr ohne Massenkraft, 
so gilt für die Geschwindigkeitsverteilung, wenn x in 
Richtung der Rohrachse gelegt wird, bekanntlich: 


Va a a nee 
Hier ist 3* das konstante Gefälle, g die Fallbeschleuni- 
gung, a = Rohrdurchmesser, r = Abstand von der 
Rohrachse. Sind e, und e, die Einheitsvektoren in 
Richtung von z und in Richtung des Aufpunktradius, 
so gilt nach Gl. (8): 


2 


® vn 
S=uonl2 +2) nen a 3 


Führt man die Werte aus Gl. (12) ein, so erhält man: 


ee Fer; 09? 3*2 (a? — r?)? 
S=& 4v (g AP+ 32»? 
par (13). 
sm 
+ er 892 r(a FEs 


Es besteht also nicht nur eine Energieströmung in 
Richtung der Flüssigkeitsströmung, sondern auch eine 
radial nach außen gerichtete Komponente des Energie- 
transportvektors. Diese Radialkomponente besitzt 


‘ein Maximum für r = = . Im Gebiet 0O<r< EG 


strömt demnach radial nach außen mehr Energie ab 


: = Ä N a. 
als von innen zuströmt, im Gebiet a > r > — ist es 


y3 
„wird also ständig 


= 


von der energiereicheren Zone0 <r < nz gespeist. 
y 


umgekehrt. Die Zone a>r > 


Bei reibungsfreien, kompressibeln Strömungen wird 
der Energietransportvektor nach wie vor durch Gl. (2) 
dargestellt, da die Kompressibilität weder die durch 
das Flächenelement hindurchtretende potentielle und 
kinetische Energie, noch die dabei nach außen ge- 
leistete Arbeit ändert. Die Divergenz von ©&, die die 
Energiesbnahme eines raumfesten Volumselementes 
bezogen auf die Volums- und Zeiteinheit mißt, nimmt 
hier allerdings eine andere Form an. Bei veränder- 
licher Dichte kann man schreiben: 


1 p p 
EN REN A LAN 
me I a ( 


Ferner gilt die allgemeine Kontinuitätsgleichung: 


ö 

tr -en=0 FAER 23 (15). 
An Stelle von GI. (3) tritt somit: 
au », 2, u! 
a wu + rer te). 


Man erhält sodann: 


Op r 
= =... + DxX (pxXb) — are) 


» , |oP AG 
Hat Ho) ren. 


Mit Hilfe der Gleichungen (14) und (15) läßt sich dieser 


- Ausdruck umformen zu: 


a N NR 
P:6=- leg +eo)+rr v0). - (16) 


N 


PP IET IT 


ERS, 


pr 
ei 
Re 
=, 
De 
E 
1 
. 
Fr 
Br 
1. 


up 


Er 


E 
5 1 
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Die Energieabnahme setzt sich also zusammen aus 
der Abnahme der kinetischen und potentiellen Energie 
und der Abnahme zufolge der Volumsdilation. Die 
letztere bleibt auch in stationärem Zustand bestehen, 
da hierbei eine thermodynamische Umsetzung statt- 
findet, bei der sich die mechanische Energie ändert. 


Liegt ein kompressibles und zähes Medium vor, 80 
tritt zum Spannungstensor Gl. (5) noch ein Glied, das 
von der Kompressibilität berrührt: 


p=—Ip+tnlprp+o,p)—AS(p:v) (1. 


Nach der Stokesschen Annahme ist A= - 


diese Annahme problematisch ist, wollen wir sie im 
folgenden nicht einführen. Für © ergibt sich nun 
folgender Ausdruck: 


2 
o-(5+ 9 +D)eo-nv-(no+u,n) | (18). 


+Ab(p-o) 
Auch hier ist © wieder gegen die Stromlinien geneigt. 


n. Da 


Nach der allgemeinen 'Gleichung von Stokes- 
Navier, die auch die Kompressibilität berück- 
sichtigt, gilt statt Gl. (10): 


2 
atr(a)-vrewxu-- rn yU 
@ 
(19). 


Herten) ren 


Mit Benutzung von Gl.(14) und (15) können wir 
daraus wieder die Divergenz von & ermitteln. Es 
ergibt sich: 


al +eo)+ rw. 


— n(p,0)- (1,0 +0,p) 
+A(p 0)? 


Außer den in Gl. (16) aufscheinenden Gliedern erhält 


(20). 


man noch den der Bee euipenian durch Reibung. 


entsprechenden Term nach .(11), und ein von 
der Kompression herrührendes dissipatives Glied. 


Wien. G. Heinrich. 


Zur Beurteilung empirischer Funktionen. 


Wenn man die Vereinbarkeit einer empirischen 
Kurve mit einer theoretischen Kurve beurteilen will, 
bedient man sich gewöhnlich der Pearsonschen 
ranNe: Da der Ausdruck y? eine Summenbildung 

er Quadrate der Abweichungen zwischen der empi- 
rischen und der theoretischen Kurve darstellt, hat die 
Methode Mängel, die recht störend wirken!. Einmal 
ist die an unempfindlich gegen das Vor- 
zeichen der Unterschiede zwischen den beiden Kurven. 
Zum anderen gestattet sie keine Prüfung auf reine 
Zufallsmäßigkeit der Unterschiede. In einem Fall, in 
dem zwischen der empirischen und der theoretischen 
Kurve rein zufallsmäßige Unterschiede auftreten, 
liefert sie denselben reinstimmungsgrad wie in 
einem Fall mit systematisch sich verändernden Unter- 
schieden, sofern der Ausdruck y? in beiden Fällen den- 
selben numerischen Wert liefert. Außerdem ist die 


*) Vgl Mittmann, O,, Erbbiologische Fragen in mathematischer 
B Sunndiang,. Walter de Gruyter, Berlin 1940, Kap, III, Abschn. 
‚8.148. 
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x?-Methode auch nur für empirische Kurven geschaf- 
fen, die als Häufigkeitsverteilungen angesehen werden 
können. Die folgenden Betrachtungen mögen dazu 
beitragen, die aufgezeigte Lücke schließen zu helfen. 

Haben wir irgendeinen Vorgang, etwa einen natur- 
wissenschaftlichen oder wirtschaftlichen Vorgang, 
nach einem fortlaufenden Argument wie z.B. nach der 
Zeit beobachtet und einen ‚Trend‘ ausfindig gemacht, 
dann wollen wir fragen, ob die Unterschiede zwischen 
dem Trend und dem empirisch beobachteten Ablauf 
reine Zufallsschwankungen sind oder eine nicht zu 
vernachlässigende systematische Komponente ent- 
halten. Wir können die Frage auch anders stellen und 
danach fragen, wie klein eine künstliche Klassenein- 
teilung des Beobachtungs-Argument-Abschnitts ge- 
wählt werden darf, damit die Unterschiede zwischen 
Trend und empirischem Ablauf noch als Zufalls- 
abweichungen angesehen werden dürfen. Der Trend 
selbst soll uns dabei nicht mehr interessieren, sondern 
nur noch die Unterschiede zwischen empirischer 
Funktion und Trend, die wir als Zufallsvariable y 
betrachten wollen. 

Zufallsvariable mit ganzzahligem Argu- 
ment. Wir gehen zunächst von einer Zufallsvariableny 
aus, die nach einem ganzzahligen Argument i beob- 
achtet sein möge. Die aufeinander folgenden Beob- 
achtungen seien Yı, Y» Y3; ---; Yx. Der Erwartungs- 
wert des einzelnen y;, den wir mit E(y;) bezeichnen 
wollen, soll sich systematisch mit i verändern dürfen. 
Die in y; enthaltene Zufallskomponente sei eine Zu- 
fallsvariable z;, die unabhängig von i sein möge. Der 
Erwartungswert E(y;) und die reine Zufallsvariable z; 
sollen sich additiv zusammensetzen, so daß 


i=Ely)+ta- 


Das einzelne z; möge sich aus unendlich vielen Teil- 
variablen additiv zusammensetzen, so daß die 2;-Ver- 
teilung eine Gaußverteilung (Normalverteilung) ist. 
Die Teilvariablen und mit ihnen die z; seien unab- 
hängig voneinander. Es soll ein Kriterium entwickelt 
werden, das zu prüfen gestattet, ob die Beobachtungs- 
serie der y; noch mit der Hypothese reiner Zufalls- 
mäßigkeit vereinbar ist. 


Falls reine Zufallsmäßigkeit herrscht, wäre Z(y;)=0 
für alle i, und alle Rohkorrelationen zwischen den y; 
wären gleich Null: R(y.,y;)=0 für i=+#+j. Falls hin- 
gegen keine reine Zufallsmäßigkeit herrscht, wären 
nicht alle Rohkorrelationen R(y:;,y;) gleich Null. 
Nach den gemachten Voraussetzungen sind alle 
R(yi,yi+m) bei gleichem m einander gleich. Die}Roh- 


‘korrelationen R(yi, ve sind praktisch natürlich 


durch ihre empirischen Repräsentanten 


1 
Amne'Tnn Zaun 


zu ersetzen. Es dreht sich dann praktisch darum, die 
Frage zu entscheiden, ob die 2, sämtlich mit dem 
Wert Null vereinbar sind oder nicht. 


Unter den oben gemachten Voraussetzungen sind 
wegen 
R(zi* 2, 2° 2) = R(zi, 2x) R(ej, zı) 
+ R(zi, 21) R(2j, 2x) 
alle R(Z, Zu) für m=#n gleich Null, und das Streu- 
ungsquadrat von Zy ist 
Z 
Mm) TN_m' 


Da die Z„, ebenso wie die 2;, einer Gaußverteilung 
folgen, wird die Wahrscheinlichkeit dafür, daß die 
normierten Variablen 


EEE en 
o(Z,)" 0(Z,) (2) o(ZX—ı) 


a 
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bestimmte Werte annehmen, dann einfach gleich 
De et 1. 3 
Do, en nD dam 
m=1 


Die erzeugende Funktion der auf den Punkt Null 
bezogenen Momente von Q lautet 
Bed) =u1-, O2, 

die Verteilung von Q ist die y?-Verteilung. Speziell 
ist der Erwartungswert Z(Q)= N —1 und das Streu- 
ungsquadrat 0?(Q)=2(N —1). Bei Benutzung der 
Pearsonschen y?-Tabelle ist zu beachten, daß die 
Zahl der ‚‚Freiheitsgrade‘ hier gleich N — 1 ist. 

Zufallsvariable mit stetigem Argument. Wir 
gehen nun zu einer Zufallsvariablen y(x) über, deren 
Werte sich nach einem stetigen Argument x ändern 
und im Intervall von e=0 bis #=1 vorliegen mögen. 
An Stelle von Z„, bilden wir jetzt den Ausdruck 

1— 


1 
Zi= | y(z) y(c+A)de 
0 
mit O0<A<1. Teilen wir das Intervall von 2=0 
bis @=1 in K gleichlange Teilintervalle (Klassen), 
soist das Streuungsquadrat von Z, entsprechend 
2 


Z 
= Run’ 
und @ erhält nun die Gestalt 


1 
9-5 aa [a-nza. 
A v 
0 


Bei unbestimmter Klassenanzahl K können wir K 
so bestimmen, daß Q@ mit der Hypothese reiner Zu- 
fallsmäßigkeit noch vereinbar ist. Wir benutzen hier- 
bei die Tatsache, daß die y?-Verteilung bei großer 
Zahl der Freiheitsgrade in eine Gau ßverteilung über- 
geht. Wir brauchen nunmehr bloß zu fragen, wann 

Q—-E(Q)|=ec:0(Q) mit c=-2bis3. 
Da E(Q)=K-—1l.und 0?) =2(K—]), ergibt 
sich der gesuchte K-Wert aus der Gleichung 


1 
ee | 1—1)Zu— (K—1)=cV2(K—ı). 
0 


Setzen wir K = J? + 1 und schreiben zur Abkürzung 
1 


S= (EAN ZUR, 
ö 


so kommen wir zur kubischen Gleichung 


7 - 23 
3 ® Pl re ER 
JE) ( 2) 12% 0, 


2 
woraus wir, wenn wir zur Abkürzunga =1 — 2 und 
ER 
b=—c (2% schreiben, sofort als Lösung erhalten 


b\2 a 
-_ eu 
undl=> EN BE mit a ee 
3] 3 Fol 
De 
3 
b 2 a 3 
falls (5 +(3) 6. 


Wählen wir eine Klassenzahl, die kleiner als die vor- 
stehend durch J berechenbare ist, dann ist die An- 
nahme reiner Zufallsmäßigkeit vertretbar. Das Streu- 
ungsquadrat der Klasseninhalte ist hier dasselbe wie 
das mittlere Fehlerquadrat der Klasseninhalte. Wählen 
wir dagegen eine Klassenanzahl, die größer als die 
vorstehend dureh J berechenbare ist, dann darf die 
systematische Komponente gegenüber der Zufalls- 
komponente nicht vernachlässigt werden. Das Streu- 
ungsquadrat der Klasseninhalte ist. dann nicht das- 
selbe wie das den Math. Statistiker interessierende 
mittlere Fehlerquadrat. 


Göttingen. Ot£frid M. J. Mittmann, 


BUCHBESPRECHUNGEN 


W.Nowacki (a0. Prof. a. d. Universität Bern) 
Fouriersynthese von Kristallenund 
ihre Anwendnngin der Chemie, (Lehr- 
bücher und Monographien aus dem Gebiete der exak- 
ten Naturwissenschaften. Chemische Reihe Bd. VI.) 
237 S. m. 120 Abb. Basel 1952. Verlag Birkhäuser. 
Preis brosch. 30,15 Fr., geb. 34,30 Fr. 

Die mehrdimensionale Fouriersynthese von Kri- 
stallen ist eine der Methoden zur Berechnung der Pa- 
rameter einer Kristallstruktur. Sie soll mittels vor- 


- liegender experimenteller Daten die Elektronendichte 


in einer Elementarzelle bestimmen. Außer dem kurz 


vorher erschienenen Buch von Booth ist das vor- 
‚liegende das einzige, das Theorie und Praxis dieser 


Verfahren behandelt. Esist aus einer Vorlesung des 


- Verfassers hervorgegangen. Die Leser dieser Zeit- 


schrift wird vor allem der zweite Teil des Buches in- 


 teressieren,in dem die Rechenhilfsmittel für die Durch- 


führung einer solehen Synthese besprochen werden, 
Hier werden zunächst die verschiedenen Streifen- 
methoden behandelt, ferner die Verwendung von Loch- 
karten, weiter die von mechanischen Rechenmaschi- 


nen und endlich die von elektromechanischen und 


elektronischen Rechenhilfsmitteln. Für alle diese Ver- 
fahren sind eine Reihe verschiedener Methoden vor- 
geschlagen, die hier sämtlich erläutert werden. Es 
folgen die optischen Methoden und, ganz kurz, die Ver- 
wendung von Mikrowellen zur Strukturbestimmung, 
Die von de Beauclair berechneten Phasenfaktoren- 
tafeln und die sonstigen vonihm vorgeschlagenen Ver- 
fahren sind noch nicht berücksichtigt. Ausführliche 
Literaturverzeichnisse am Anfang eines jeden Ab- 
schnittes erleichtern die eingehendere Beschäftigung 
mit den einzelnen Methoden. Ein Tabellenverzeichnis, 
ein Namenregister, ein Sachregister und ein Substanz- 
register machen das Buch auch zum Nachschlagewerk 
geeignet. 


Dresden. ‚Willers. 


Prof. Dr. W.H. Westphal, Physikalisches 
Wörterbuch. Zwei Teile in einem Band, VII + 
831 + 795 8, mit 1595 Abb. Berlin-Göttingen-Heidel- 
berg 1952. Springer-Verlag. Preis geb. 148,— DM. 

Das vorliegende Wörterbuch soll das vor ungefähr 
zwanzig Jahren von Berliner und Scheel herausgege- 
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bene ersetzen. Als Benutzer ist an Physiker gedacht, 
die sich über Fragen orientieren wollen, die nicht in 
ihr Spezialgebiet fallen, ferner an Mathematiker, In- 
genieure und Chemiker. Die einzelnen Artikel sind so 
abgefaßt, daß sie eine erste Übersicht geben und in 
vielen Fällen die Heranziehnug von Spezialliteratur 
überflüssig machen. Abgesehen davon, daß selbst- 
verständlich die Fortschritte der Physik seit dem Er- 
scheinen jenes Handwörterbuches berücksichtigt sind 
— ist doch z. B. die eigentliche Kernphysik erst seit 
jener Zeit entstanden — sind in dem neuen Buch die 
physikalische Chemie, die Astrophysik, in gewissem 
Umfang auch die Geophysik und die Biophysik be- 
rücksichtigt; ferner ist das aus der Mathematik auf- 
genommen, was der Physiker laufend gebraucht. Um 
den Umfang des Bandes nicht allzugroß werden zu 
lassen, hat man sich im wesentlichen auf die Grund- 
lagenforschung beschränkt und die angewandte Phy- 
sik nur soweit berücksichtigt, wie es ihrer Bedeutung 
für die Grundlagenforschung entspricht. 

80 Wissenschaftler haben sich an der Arbeit be- 
teiligt. Der Stoff ist in rund 10 500 Stichworte auf- 
gegliedert. Dazu kommen drei Anhänge, die eine 
kurze Geschichte der Physik, eine Zusammenstellung 
der Lebensdaten von Physikern, auch von noch leben- 
den, und Tabellen wichtiger Konstanten usw. ent- 
halten. Man muß dem Herausgeber dankbar sein, daß 
er sich der Mühe der endgültigen Üßerarbeitung und 
der Abstimmung der einzelnen Artikel auf einander 
unterzogen hat. Es scheint mir das gut gelungen zu 
sein. Das bei der Fülle des Stoffes einige Versehen 
vorkommen, ist nicht zu vermeiden. Stichproben 
haben z. B. folgende Kleinigkeiten ergeben: Die For- 
meln für die Koeffizienten der Newtonschen Inter- 
polationsformel sind nicht richtig. Bei der Figur zum 
Newtonschen Näherungsverfahren ist von einemPunkt 
ausgegangen, bei dem die Vorzeichen von f(x) und 
f’’(z) nicht gleich sind, so daß die Konvergenz nicht 
gesichert ist. Bei den Rollkurven ist die Abbildung 
der Kardioide falsch. In dem Artikel über Winkel- 
geschwindigkeit sind in Abb. 3 die Vektoren falsch ein- 
gezeichnet. Das Bild des Poinsot-Ellipsoides ist 
schlecht. Bei den Potenzreihen wird im Reellen aus 
der Konvergenz für = x, auf die für |2| <|x|, 
also auch für — x, geschlossen; im Komplexen wird 
dieser Fehler nicht gemacht. In Bd. II S. 699 Spalte 2 
Zeile 5 muß es heißen 2.,1 =... = 2%,1= 2. Aber 
das sind gegenüber dem vielen Ausgezeichneten, was 
das Buch bietet, unbedeutende Kleinigkeiten. Jeden- 
falls ist in dem physikalischen Wörterbuch für 
Physiker und Wissenschaftler der Nachbargebiete 
ein außerordentlich wertvolles Hilfsmittel geschaffen 
worden, das allgemein begrüßt und viel benutzt wer- 
den wird. 


Dresden. Willers. 


J. P. den Hartog (Professor am Mass. Inst. of 
Technology) u. &. Mesmer (Professor a. d. Techn. 
Hochschule Darmstadt.) Mechanische Schwin- 

ungen. Zweite Auflage. XVI--4278. mit 
299 Abb. Berlin-Göttingen-Heidelberg 1952. Sprin- 
ger-Verlag. Preis geb. 42,— DM. 

Die deutsche Erstauflage des Buches, die im Jahre 
1936 erschienen war, hatte großen Anklang gefunden, 
da Verf. es versteht, auch bei komplizierteren Vor- 
gängen den Leser auf anschauliche Art in die wesent- 
lichsten Zusammenhänge einzuführen. Es wird dabei 
zwar nicht ein höherer theoretischer Überblick ver- 
mittelt, aber in den meisten Fällen ein einfacher 
Rechnungsweg gefunden, der für die praktische 
Schwingungstechnik zu befriedigenden Ergebnissen 
führt. Insbesondere haben wohl auch die zahlreichen 
Anwendungsbeispiele und Aufgaben dazu beigetragen, 
dem Buch in weiteren Kreisen Freunde zu verschaffen. 
Nachdem inzwischen die dritte amerikanische Auflage 
erschienen war, wird nunmehr vom Bearbeiter und 

tsetzer eine zweite deutsche Auflage vorgelegt, 


wobei infolge Zeitmangel der sehr anschaulich ge- 
schriebene amerikanische Originaltext im wesentlichen 
heibehalten worden ist. Es ist zu erwarten, daß auch 
die neue Auflage großen Anklang finden wird. 


Dresden, H. Neuber. 


H. Ziegler (Prof. a. d. ETH. Zürich). Mechanik, 
Bd.III. Dynamik der Systeme (Lehr- und Hand- 
bücher der Ingenieurwissenschaften Nr. 7.) 396 S. mit 
191 Abb. Basel 1952. Verlag Birkhäuser. Preis geb. 
46,80 schw. fr. 

Der vorliegende dritte Band der Sammlung Mecha- 
nik, der von den Professoren E.Meissner und 
H. Ziegler gemeinsam herausgegeben wurde, be- 
handelt die Dynamik der Systeme (bzw. des Punkt- 
haufens). Band behandelte die Statik der starren 
Flüssigkeiten und elastischen Körper und erschien im 
Jahre 1948; Band 2 behandelte die Dynamik der 
starren Körper und erschien im Jahre 1947. Bei der 
Behandlung der Dynamik der Systeme stellt Verf. 
entsprechend der modernen Auffassung die La- 
grangeschen Gleichungen in den Vordergrund 
(Nach dem Vorwort gehören die Lagrangeschen 
Gleichungen an der Eidgenössischen Technischen 
Hochschule Zürich, an der der Verf. als Professor der 
technischen Mechanik tätigist, zum allgemeinen Rüst- 
zeug des heutigen Ingenieurs). Im 1. Kapitel wird die 


Dynamik von Systemen starrer Körper entwickelt, 


wobei auch das d’Alembertsche Prinzip und das 
Prinzip der virtuellen Leistungen im Vordergrund 
stehen. Es enthält ferner neben den allgemeinen 
Sätzen über Schwerpunkt, Impuls, Drehimpuls, 
Energie auch Einiges über freie und erzwungene 
Schwingungen. Eine ausführlichere Behandlung der 
Schwingungsprobleme ist Gegenstand des zweiten 
Abschnittes, der sich in folgende Paragraphen gliedert: 
Schwingungen mit zwei Freiheitsgraden, Normal- 
koordinaten, erzwungene Schwingungen, Schwin- 
gungen mit endlichem Freiheitsgrad, das Rayleigh- 
sche Prinzip, Stabilitätsprobleme, einfache kontinuier- 
liche Schwinger, Schwingungen und Wellen, die Eigen- 
schwingungen, erzwungene Schwingungen, das Ray- 
leighsche Verfahren, Querschwingungen von Stäben, 
erzwungene Querschwingungen. Es wird daher in 
diesem Kapitel eine Reihe technisch wichtiger 
Schwingungsprobleme erörtert, welche sich bei Syste- 
men mit endlichem Freiheitsgrad ergeben. Das 
3. Kapitel behandelt die Systeme mit unendlich 
hohem Freiheitsgrad bzw. die Kontinua, wobei eine 
Reihe von Verfahren der Elastizitätstheorie und 
Hydrodynamik zur Geltung kommen. Es gliedert sich 
in folgende Paragraphen: Der räumliche Spannungs- 
zustand, Verschiebungen und Verzerrungen, Span- 
nungen und Verzerrungen, feste elastische Körper, 
die elastostatischen Grundgleichungen, die Aufgaben 
der Festigkeitslehre, ebene Probleme, Schwi in 
fester Körper, Plastizitätsbedingungen, Flüssigkeiten 
und Gase, die hydrodynamischen Er 
Strom- und Wirbelfelder, einfache hydrodynamische 
Probleme, ebene Potentialströmungen, der unend- 
lich lange Tragflügel, Tragflügel endlicher Länge, 
zähe Flüssigkeiten. Vielfach wird bei der Darstellung 
von der Tensorrechnung Gebrauch gemacht, welche 
fortlaufend mitentwickelt wird und schließlich in 
einem Anhang noch einmal zusammengefaßt er- 
scheint. . 

Beim Studium des Buches fällt auf, daß der/Verf. 
die Mechanik mehr vom Standpunkte des Mathe- 
matikers darlegt und kaum Beispiele bringt. Wie im 
Vorwort erwähnt wird, handelt es sich insbesondere 
beim zweiten und dritten Abschnitt um Wiedergabe 
von Vorlesungen, die für Studierende höherer Se- 
mester bestimmt sind. Der Fachmann wird natürlich 
die konsequente mathematische Durchdri mit 
Befriedigung wahrnehmen. Es erscheint jedoch als 
ein gewisser Mangel, daß nur wenige im Text ver- 
streute Anwendungsbeispiele angegeben sind, .die 


a 


Kr? 
er 


. Halle. - 


cs = 


lichkeiten und 


Bi 


Z. angew. Math. Mech. 
Bd. 32 Nr. 8/9 Aug./Sept. 1952 


*  Buchbesprechungen 


291 


außerdem sehr stark idealisiert erscheinen und sich 
wenig an die technische Praxis anlehnen. Literatur- 
hinweise fehlen fast völlig, lediglich am Ende des 
Buches ist eine Zusammenstellung von Werken bei- 
gefügt, die dem Inhalt nahestehen. Trotz dieser 
Mängel wird das Buch infolge seines hohen theoreti- 
schen Wertes von weiten Kreisen sehr begrüßt werden. 


Dresden. H. Neuber. 


Prof. L. Balser, Einführung in die Karten- 
lehre (Kartennetze). 2.Aufl. 648. mit 
50 Abb. Leipzig 1951. B. G. Teubner Verlagsgesell- 
schaft, Preis 3,— DM. 

Das vorliegende Bändchen bringt die wichtigsten 
Kartenentwürfe und diskutiert deren Eigenschaften. 
Die Darstellung ist elementar und macht nur an 
wenigen Stellen von der Integralrechnung Gebrauch. 
Zur . Veranschaulichung der Längenverzerrungen 
wird für die wichtigsten Entwürfe an einigen Stellen 
der Karte die Indikatrix, das Bild eines kleinen Kreises 
auf der Kugel, angegeben. Nach der Erörterung von 
flächentreuen und winkeltreuen Entwürfen wird im 
letzten Abschnitt auf einige Sonderentwürfe (gno- 
monische Projektion, Entwürfe von Aitoff, Hammer 
und Winkel) und Ergänzungen (Umbezifferung, Be- 
rücksichtigung der Abplattung, Gauß-Krüger- 
sche Koordinaten) eingegangen. Das Büchlein steht 
durch seine anschauliche und elementare Darstellungs- 
art einem weiten Leserkreis offen und kann jedem 
empföhlen werden, der in die Entstehung der Land- 
und Sternkarten einen ersten Einblick gewinnen will. 


Dresden, G. Opitz. 


"Collegue de Geomsötrie dififäsren- 
tielle tenu & Louvain du 11. au 14. avril 1951. 
(Centre belge de recherches mathömatiques.) 235 8. 
Liege et Paris 1951. Georges Thone, Masson & Cie. 
Preis brosch 350 belg. fr., 2450 franz. fr. 

Wie die mathematische Welt es nun schon seit 
einigen Jahren gewohnt ist, veranstalteten die Lüt- 


. - ticher Mathematiker auch in diesem Jahr ein Collo- 
quium, in dem führende Männer aus aller Welt zu 


Worte kamen. — Nachdem in denletzten Jahrzehnten 
eine projektive Differentialgeometrie sich gleichwertig 
neben die klassische metrische gestellt hat, ist es nun 
erstaunlich zu sehen, wie sich in den neueren For- 
schungen mehr und mehr topologische Gesichtspunkte 
in den Vordergrund schieben und bei der Behandlung 
klassischer Probleme ihre Fruchtbarkeit erweisen. 
Andere Vorträge gaben neue Ergebnisse über Rie- 


lität), Steiner (Doppelverhältnis, Fundamentalsatz), 
v. Staudt (harmonische Würfe) vor uns hin. 

Das Büchlein gibt gleichzeitig dem Laien einen ersten 
Einblick in das Gebiet, wie dem Kenner Freude und 
Anregung. 


Halle. Ott-Heinrich Keller. 


D. J. Struik, (Prof. a. Massachusetts Inst. of Tech- 
nology). LecturesonClassical Differen- 
tial Geometry. VIII + 221 S. mit 120 Abb. 
Cambridge (Mass.) 1950. Addison Wesley Press. 
Preis geb. 6,— $. 

Das Buch ist entstanden aus Vorlesungen am Mas- 
sachusetts Institute of Technology. Es bringt die 
klassischen Sätze über Kurven, Flächen, Geometrie 
auf der Fläche, Grundgleichungen und besondere 
Flächen (Minimalflächen, Einhüllende, Regelflächen 
u. &.). Die Figuren sind anschaulich und instruktiv. 
Auf die Auswahl der Übungsbeispiele ist große Sorg- 
falt verwendet; sie sind zahlreich und steigen in jedem 
Kapitel vom Leichteren zum Schwereren auf. Die 
Leistungen der Großen werden in ihrer historischen 
Entwicklung und Bedingtheit dem Leser einprägsam 
vor Augen geführt. Es werden überall die Fragen, die 
über den Rahmen des Buches hinausweisen, aufgezeigt 


und auf die Originalarbeiten verwiesen, so daß auch 


der Student, der sein Studium vertiefen will, mit Ge- 
winn zu dem Buche greifen wird. 


Halle. Ott-HeinrichKeller. 


Dr. P. Zühlke (Honorarprof. a. d. Univ. Marburg), 
KonstruktioneninbegrenzterEbene 
(Mathematisch-physikalischeBibliothek Reihel Nr.11.) 
3. Aufl. 428. mit 65 Abb. Leipzig 1951. B.G. Teubner 
Verlagsgesellschaft. Preis kart. 2,10 DM, 

Verf. will zum Nachdenken darüber anregen, wie 
man sich beim Zeichnen zu helfen habe, wenn wesent- 
liche Punkte außerhalb des Zeichenbrettes fallen. Er 
behandelt: Unzugängliche Punkte von zwei und mehr 
Geraden, Halbierung eines Winkels mit unzugäng- 
lichem Scheitel, unzugängliche Punkte von Kreisen. 
Verf. weist auf den Zusammenhang dieser Aufgaben 


mit denen mit beschränkten Hilfsmitteln, etwa zu 


kleiner Zirkel oder zu kurzer Lineale hin, 

Das Buch wendet sich nicht nur an den, der prak 
tisch zu zeichnen hat, sondern auch der Liebhaber 
und der Lehrer der Geometrie wird aus den vielen 
echt geometrischen Fragen und Gedanken des Büch- 
leins wertvolle Anregung gewinnen, 


Haie, DE Heinz oh elrer 


 mannsche, Einsteinsche, elektromagnetische - 


- Geometrien. Der Bericht wird in erster Linie den For- 
scher auf dem Gebiet der Differentialgeometrie inte- 


ressieren, aber auch der Mathematiker, der nicht un- 
mittelbar auf diesem Gebiet arbeitet, wird sich an 


' Hand dieses Berichtes einen Überblick verschaffen 


können, welche Fragen heute im Vordergrund des In- 
teresses stehen. 
Ott-Heinrich Keller. 


Prof. Dr. M. Zacharias (wissenschaftl. Mitarbeiter 
der Deutschen Akademie der Wissenschaften zu Ber- 


_ lin, EinführungindieprojektiveGeo-. 


_ metrie. (Mathematisch-physikalische Bibliothek 


Reihe 1 Nr. 6.) Vierte Aufl. 54 S. mit 21 Abb. Leipzig 


1951. B. G. Teubner Verlagsgesellschaft. Preis kart. 


‚2,— DM. 


Verf. führt in einer feinen, geistvollen Art den Leser, 
der nur Schulkenntnisse besitzt, in die Schönheiten 
der projektiven Geometrie ein und vermittelt ihm die 


f 
ER Z e 
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Ehrfurcht vorihren Schöpfern, ohne dabei Systematik 
und Vollständigkeit anzustreben. Er stellt die Persön- _ 
'h Leistungen von Desargues (Zentralper- 
 spektive, uneigentliche Elemente, Polarentheorie, In- 
_ volutionen), Pascal (PascalscherSatz), Poncelet (Dua- 
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Dr. A.U. Huggenberger, Talsperren-Meß- 
technik (Meßverfahren, Instrumente und Appa- 
rate für die Prüfung der Bauwerke in Massenbeton) 
VII + 132 S. m. 168 Abb. Berlin/Göttingen/Heidel- 
berg 1951. Springer-Verlag. Preis brosch. 22,50 DM. 

Die im Talsperrenbau auftretenden meßtech- 


nischen Aufgaben bestehen einmal in der Verfolgung 


der Veränderungen, die das Bauwerk während seines 
Entstehens durch die zunehmende Gewichtsbelastung 
und die Temperatur- und Feuchtigkeitsänderungen 


während des Abbindens erleidet, und zweitens in der, 


Überwachung der bei der Inbetriebnahme einerseits 
durch den Wasserdruck, andererseits durch die at- 
'mosphärischen Einflüsse hervorgerufenen Span- 
nungen und der durch sie bewirkten Gestaltsände- 
rungen. Diese Aufgaben stellen besondere Anforde- 
rungen an die zu verwendenden Geräte und ihren 
Einbau, abgesehen von den Messungen der äußeren 
Verformung müssen die Geberteile unzugänglich ein- 
betoniert werden. Eskommen also nur fernanzeigende, 
elektrische Methoden in Frage, und die Ausbildung 
der Geräte muß eine jahrelange Betriebssicherheit 
gewährleisten, da der Ausfall schon eines Gerätes 
die Überwachung gefährden kann. Der Schwerpunkt 
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liegt demnach in der mechanisch, elektrisch und 
korrosionsfesten Gestaltung und Einbettung der Ge- 
räte und Kabel bei gleichzeitig hohen Anforderungen 
an die Empfindlichkeit und Genauigkeit der Messung. 

Im einzelnen handelt es sich innerhalb des Beton- 
körpers um Dehnungs-, Spannungs- und Temperatur- 
messungen, die auf zwei Grundformen basieren, der 
Frequenzänderung einer schwingenden Saite durch 
Längung oder Kürzung, oder der Änderung des 
elektrischen Leitwertes eines gespannten Drahtes 
mit der Spannung desselben. Die letztere Methode 
der Widerstandsmessung wird auch zur gleich- 
zeitigen oder besonderen Temperaturmessung ver- 
wendet, ebenso zur. Feuchtemessung. Eine zweite 
Gruppe umfaßt die Beobachtung des Wasserdruckes 
am Boden der Sperrmauer und des Porendruckes im 
Beton mit Hilfe manometrischer Methoden, Zur 
Messung der äußeren Verformung endlich dienen 
Lotungsgeräte, Nivelliergeräte und geodätische Tri- 
angulierungen. 

Der Verfasser beschreibt in strenger Gliederung 
die verschiedenen Verfahren, die Konstruktion .der 
einzelnen Geräte, ihren Einbau, ihre Bedienung und 
Eichung in aller Ausführlichkeit. Er läßt auch 
Nebendinge nicht unbeachtet, wo sie auf die dauernde 
Betriebssicherheit von Einfluß sein könnten. Durch 
zahlreiche Skizzen und Fotos wird die Konstruktion 
und der Einbau jedes Gerätes klar und überzeugend 
erläutert. 

Nicht nur dem Bauingenieur ist das Buch ein aus- 
gezeichneter und unentbehrlicher Berater. Es ver- 
mag auch für alle anderen Fachgebiete, insbesondere 
da, wo langjährige Betriebssicherheit unzugänglicher 
Meßstellen gefordert wird, eine Fülle wertvoller An- 
regungen zu geben und sei darum jedem Meßtechniker 
wärmstens empfohlen. 


Dresden. H. Faltin 
L, Landau und E. Lifschitz (Akademie der Wissen- 
schaften der UdSSR), The classicaltheory 
of fields. (Aus dem-Russischen ins Englische 
übersetzt von M. Hamermesh). IX -+ 3548. 
mit 12 Abb. Cambridge Mass., USA 1951. Addison- 
i Wesley Press. Preis geb. 7,50 $. 
R} Dieses Buch ist Teil eines Lehrbuches der theo- 
retischen Physik. Es behandelt das elektromagneti- 
Be sche und das Gravitationsfeld unter Ausschluß der 
2: Quantentheorie. Entgegen dem bei uns üblichen 
2 Verfahren, das vom Leichteren zum Schwereren 
fortschreitet, enthält das vorliegende Buch eine 


Dr. W. Krull (o. Prof. a. d. Univ. Bonn), Ele- 
mentare und klassische'fAlgebra 


Göschen Band 930) Bd.I, 2. Aufl. 1368. Berlin 1952. 
Verlag Walter de Gruyter & Co. Preis brosch. 2,40 DM. 


damm 19, Fernsprecher Sammelnummer: 42 55 71. 
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Eingegangene Bücher 


i Z. angew. Math. Mech. 
Bd.32 Nr. 8/9 Aug./Sept. 1952 


systematische Darstellung der Theorie. Ausgangs- 
punktist das Relativitätsprinzip und das Prinzip der 
Konstanz der Lichtgeschwindigkeit, mit denen zu- 
nächst die Kinematik der speziellen Relativitäts- 
theorie abgeleitet wird. Es folgt die Relativitäts- 
mechanik, die vom Prinzip der kleinsten Wirkung 
ausgeht. Dieses Prinzip wird dann erweitert auf den 
Fall, daß sich Ladungen im elektromagnetischen 
Feld bewegen. Auf diesem Wege werden die Feld- 
größen eingeführt. Anschließend werden die Zu- 
sammenhänge zwischen den Feldgrößen, die Max- 


‘wellschen Gleichungen, behandelt und daran an- 


geschlossen das konstante Feld, elektromagnetische 
Wellen, geometrische Optik und Wellenoptik. 
Schließlich werden die Felder und die Strahlung be- 
wegter Ladungen recht ausführlich dargestellt. Den 
Abschluß, fast ein Drittel des Buches, bildet die all- 
gemeine Relativitätstheorie (Gravitationsfelder). Die 
klare, sorgfältige Darstellung macht das Studium des 
Buches zu einem großen Genuß. Für Studenten, die 
zum ersten Mal theoretische Physik lernen, ist es 
aber trotzdem nicht geeignet, weil es zu viel Kennt- 
nisse voraussetzen muß. Wohl aber kann man es 
z. B. unseren Doktoranden und Aspiranten für ihre 
Fortbildung sehr empfehlen. 


Dresden. Recknagel. 


Prof. Dr.F.Hennig, Vontiefenundhohen 
Temperaturen. (Eine allgemein verständ- 
liche Darstellung.) 848. m. 44 Abb. Leipzig 1951. 
B.G. Teubner Verlagsgesellschaft. Preis kart. 
3,40 DM. 

Die vorliegende kleine Schrift bringt n der 
Temperaturmessung dem naturwissenschaftlich 


"bildeten Leser nahe. Mathematische Mittel größeren 


Umfanges werden nicht angewendet, unsachliche 
Popularisierung aber trotzdem sorgfältig vermieden, 

Inhalt: Einfachste Grundlagen der Temperatur- 
skala und wichtigste Temperaturmeßgeräte werden 
eingangs geschildert. Es folgen Zusammenhänge von 
Molekularbewegung und Temperatur und damit die 
genaue Definition der Temperaturskala durch das 
ideale Gas; Zusammenhänge von elektroma, ischer 
Strahlung und Temperatur und damit die optische 
Pyrometrie einschließlich astrophysikalischer Folge- 
rungen; Erzeugung und Messung tiefster Temperatur 
einschließlich thermodynamischer Definition der 
Temperaturskala. 


Dresden. Recknagel, 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Hefte eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Angerer-Ebert, TechnischeKunstgriffe 
bei physikalischen Untersuchungen. 


vom modernen Standpunkt. (Sammlung 8,, vollständig neubearbeitete und erweiterte Auf- 


lage. 350 S. mit 118 Abb. Braunschweig 1952. Ver- 
lag F. Vieweg & Sohn. Preis geb. 14,80 DM. 
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